Multiple parasites mediate balancing selection at two MHC class II genes in the fossorial water vole: insights from multivariate analyses and population genetics

Genome-wide relaxation of selection and the evolution of the island syndrome in Orkney voles

Island populations often experience different ecological and demographic conditions than their counterparts on the continent, resulting in divergent evolutionary forces affecting their genomes. Random genetic drift and selection both may leave their imprints on island populations, although the relative impact depends strongly on the specific conditions. Here we address their contributions to the island syndrome in a rodent with an unusually clear history of isolation. Common voles (Microtus arvalis) were introduced by humans on the Orkney archipelago north of Scotland >5000 years ago and rapidly evolved to exceptionally large size. Our analyses show that the genomes of Orkney voles were dominated by genetic drift, with extremely low diversity, variable Tajima's D, and very high divergence from continental conspecifics. Increased d N/d S ratios over a wide range of genes in Orkney voles indicated genome-wide relaxation of purifying selection. We found evidence of hard sweeps on key genes of the lipid metabolism pathway only in continental voles. The marked increase of body size in Orkney-a typical phenomenon of the island syndrome-may thus be associated to the relaxation of positive selection on genes related to this pathway. On the other hand, a hard sweep on immune genes of Orkney voles likely reflects the divergent ecological conditions and possibly the history of human introduction. The long-term isolated Orkney voles show that adaptive changes may still impact the evolutionary trajectories of such populations despite the pervasive consequences of genetic drift at the genome level.

Infection with a Recently Discovered Gammaherpesvirus Variant in European Badgers, Meles meles, is Associated with Higher Relative Viral Loads in Blood

Herpesviruses are ubiquitous pathogens infecting most animals. Although host immunity continually coevolves to combat virulence, viral variants with enhanced transmissibility or virulence occasionally emerge, resulting in disease burdens in host populations. Mustelid gammaherpesvirus 1 (MusGHV-1) is the only herpesvirus species identified thus far in European badgers, Meles meles. No MusGHV-1 associated pathomorbidity has been reported, but reactivation of MusGHV-1 in genital tracts is linked to impaired female reproductive success. An analysis of a short sequence from the highly conserved DNA polymerase (DNApol) gene previously identified two variants in a single host population. Here we compared genetic variance in blood samples from 66 known individuals of this same free-ranging badger population using a partial sequence comprising 2874 nucleotides of the DNApol gene, among which we identified 15 nucleotide differences resulting in 5 amino acid differences. Prevalence was 86% (59/66) for the common and 17% (11/66) for the novel variant, with 6% (4/66) of badgers presenting with coinfection. MusGHV-1 variants were distributed unevenly across the population, with individuals infected with the novel genotype clustered in 3 of 25 contiguous social groups. Individuals infected with the novel variant had significantly higher MusGHV-1 viral loads in their blood (p = 0.002) after adjusting for age (juveniles > adults, p < 0.001) and season (summer > spring and autumn, p = 0.005; mixed-effect linear regression), likely indicating higher virulence of the novel variant. Further genome-wide analyses of MusGHV-1 host resistance genes and host phenotypic variations are required to clarify the drivers and sequelae of this new MusGHV-1 variant.

Borrelia Infection in Bank Voles Myodes glareolus Is Associated With Specific DQB Haplotypes Which Affect Allelic Divergence Within Individuals

The high polymorphism of Major Histocompatibility Complex (MHC) genes is generally considered to be a result of pathogen-mediated balancing selection. Such selection may operate in the form of heterozygote advantage, and/or through specific MHC allele–pathogen interactions. Specific MHC allele–pathogen interactions may promote polymorphism via negative frequency-dependent selection (NFDS), or selection that varies in time and/or space because of variability in the composition of the pathogen community (fluctuating selection; FS). In addition, divergent allele advantage (DAA) may act on top of these forms of balancing selection, explaining the high sequence divergence between MHC alleles. DAA has primarily been thought of as an extension of heterozygote advantage. However, DAA could also work in concert with NFDS though this is yet to be tested explicitly. To evaluate the importance of DAA in pathogen-mediated balancing selection, we surveyed allelic polymorphism of MHC class II DQB genes in wild bank voles (Myodes glareolus) and tested for associations between DQB haplotypes and infection by Borrelia afzelii, a tick-transmitted bacterium causing Lyme disease in humans. We found two significant associations between DQB haplotypes and infection status: one haplotype was associated with lower risk of infection (resistance), while another was associated with higher risk of infection (susceptibility). Interestingly, allelic divergence within individuals was higher for voles with the resistance haplotype compared to other voles. In contrast, allelic divergence was lower for voles with the susceptibility haplotype than other voles. The pattern of higher allelic divergence in individuals with the resistance haplotype is consistent with NFDS favouring divergent alleles in a natural population, hence selection where DAA works in concert with NFDS.

Pathogen-mediated selection favours the maintenance of innate immunity gene polymorphism in a widespread wild ungulate

Toll-like receptors (TLR) play a central role in recognition and host frontline defence against a wide range of pathogens. A number of recent studies have shown that TLR genes (Tlrs) often exhibit large polymorphism in natural populations. Yet, there is little knowledge on how this polymorphism is maintained and how it influences disease susceptibility in the wild. In previous work, we showed that some Tlrs exhibit similarly high levels of genetic diversity as genes of the Major Histocompatibility Complex (MHC), and signatures of contemporary balancing selection in roe deer (Capreolus capreolus), the most abundant cervid species in Europe. Here, we investigated the evolutionary mechanisms by which pathogen-mediated selection could shape this innate immunity genetic diversity by examining the relationships between Tlr (Tlr2, Tlr4 and Tlr5) genotypes (heterozygosity status and presence of specific alleles) and infections with Toxoplasma and Chlamydia, two widespread intracellular pathogens known to cause reproductive failure in ungulates. We showed that Toxoplasma and Chlamydia exposures vary significantly across years and landscape features with few co-infection events detected and that the two pathogens exert antagonistic selection on Tlr2 polymorphism. By contrast, we found limited support for Tlr heterozygote advantage. Our study confirmed the importance of looking beyond Mhc genes in wildlife immunogenetic studies. It also emphasized the necessity to consider multiple pathogen challenges and their spatiotemporal variation to improve our understanding of vertebrate defence evolution against pathogens.

MHC-associated Baylisascaris schroederi load informs the giant panda reintroduction program

Reintroducing captive giant pandas (Ailuropoda melanoleuca) to the wild is the ultimate goal of their ex situ conservation. Choosing higher fitness candidates to train prior to release is the first step in the giant panda reintroduction program. Disease resistance is one important index of individual fitness and presumed to be related to variation at major histocompatibility complex genes (MHC). Here, we used seven polymorphic functional MHC genes (Aime-C, Aime-I, Aime-L, Aime-DQA1, Aime-DQA2, Aime-DQB1 and Aime-DRB3) and estimate their relationship with Baylisascaris schroederi (Ascarididae) infection in giant panda. We found that DQA1 heterozygous pandas were less frequently infected than homozygotes. The presence of one MHC genotype and one MHC allele were also associated with B. schroederi infection: Aime-C*0203 and Aime-L*08 were both associated with B. schroederi resistance. Our results indicate that both heterozygosity and certain MHC variants are important for panda disease resistance, and should therefore be considered in future reintroduction programs for this species alongside conventional selection criteria (such as physical condition and pedigree-based information).

•

MHC heterozygous pandas were less frequently infected by Baylisascaris schroederi than homozygotes.

•

Presence of Aime-C*0203 and Aime-L*08 are associated with Baylisascaris schroederi resistance.

•

MHC types are important for panda parasite resistance.

Post-glacial phylogeography and variation in innate immunity loci in a sylvatic rodent, bank vole Myodes glareolus

In the northern hemisphere, the spatial structure of many taxa has been shaped by migration patterns after the last glaciation, and phylogeography based on mtDNA variation may reflect the post-glacial demography. The mtDNA lineages are expected to differ in their adaptations to local conditions but little is known about the impact of these conditions on functional genetic variation. Here, we answer this question through an analysis of geographic variation and selection patterns in seven innate immunity genes in free-living bank voles Myodes glareolus from 10 localities across species range assigned to different lineages based on mtDNA. We found clear discrepancies between population structure in mtDNA and each of the studied innate immunity genes. There was no uniform pattern of spatial variation at immunity loci, they differed in the levels of polymorphism, and the results of neutrality tests were not consistent over loci. Each locus comprised a few common haplotypes shared between mitochondrial lineages and studied locations, plus numerous haplotypes unique for each studied site. Our results suggest that the diversity of innate immunity genes cannot be explained solely in terms of demographic processes, and that the observed polymorphism may be attributed to local selection. The strength and direction of selection differed between loci, even within the same gene family, which underlines how crucial it is to take a complex approach while studying the selection patterns acting on immune-related genes.

Both candidate gene and neutral genetic diversity correlate with parasite resistance in female Mediterranean mouflon

BACKGROUND

Parasite infections can have substantial impacts on population dynamics and are accordingly a key challenge for wild population management. Here we studied genetic mechanisms driving parasite resistance in a large herbivore through a comprehensive approach combining measurements of neutral (16 microsatellites) and adaptive (MHC DRB1 exon 2) genetic diversity and two types of gastrointestinal parasites (nematodes and coccidia).

RESULTS

While accounting for other extrinsic and intrinsic predictors known to impact parasite load, we show that both neutral genetic diversity and DRB1 are associated with resistance to gastrointestinal nematodes. Intermediate levels of multi-locus heterozygosity maximized nematodes resistance, suggesting that both in- and outbreeding depression might occur in the population. DRB1 heterozygosity and specific alleles effects were detected, suggesting the occurrence of heterozygote advantage, rare-allele effects and/or fluctuating selection. On the contrary, no association was detected between genetic diversity and resistance to coccidia, indicating that different parasite classes are impacted by different genetic drivers.

CONCLUSIONS

This study provides important insights for large herbivores and wild sheep pathogen management, and in particular suggests that factors likely to impact genetic diversity and allelic frequencies, including global changes, are also expected to impact parasite resistance.

Shared evolutionary origin of major histocompatibility complex polymorphism in sympatric lemurs

Genes of the major histocompatibility complex (MHC) play a central role in adaptive immune responses of vertebrates. They exhibit remarkable polymorphism, often crossing species boundaries with similar alleles or allelic motifs shared across species. This pattern may reflect parallel parasite-mediated selective pressures, either favouring the long maintenance of ancestral MHC allelic lineages across successive speciation events by balancing selection ("trans-species polymorphism"), or alternatively favouring the independent emergence of functionally similar alleles post-speciation via convergent evolution. Here, we investigate the origins of MHC similarity across several species of dwarf and mouse lemurs (Cheirogaleidae). We examined MHC class II variation in two highly polymorphic loci (DRB, DQB) and evaluated the overlap of gut-parasite communities in four sympatric lemurs. We tested for parasite-MHC associations across species to determine whether similar parasite pressures may select for similar MHC alleles in different species. Next, we integrated our MHC data with those previously obtained from other Cheirogaleidae to investigate the relative contribution of convergent evolution and co-ancestry to shared MHC polymorphism by contrasting patterns of codon usage at functional vs. neutral sites. Our results indicate that parasites shared across species may select for functionally similar MHC alleles, implying that the dynamics of MHC-parasite co-evolution should be envisaged at the community level. We further show that balancing selection maintaining trans-species polymorphism, rather than convergent evolution, is the primary mechanism explaining shared MHC sequence motifs between species that diverged up to 30 million years ago.

Genotype Calling from Population-Genomic Sequencing Data

Genotype calling plays important roles in population-genomic studies, which have been greatly accelerated by sequencing technologies. To take full advantage of the resultant information, we have developed maximum-likelihood (ML) methods for calling genotypes from high-throughput sequencing data. As the statistical uncertainties associated with sequencing data depend on depths of coverage, we have developed two types of genotype callers. One approach is appropriate for low-coverage sequencing data, and incorporates population-level information on genotype frequencies and error rates pre-estimated by an ML method. Performance evaluation using computer simulations and human data shows that the proposed framework yields less biased estimates of allele frequencies and more accurate genotype calls than current widely used methods. Another type of genotype caller applies to high-coverage sequencing data, requires no prior genotype-frequency estimates, and makes no assumption on the number of alleles at a polymorphic site. Using computer simulations, we determine the depth of coverage necessary to accurately characterize polymorphisms using this second method. We applied the proposed method to high-coverage (mean 18×) sequencing data of 83 clones from a population of Daphnia pulex. The results show that the proposed method enables conservative and reasonably powerful detection of polymorphisms with arbitrary numbers of alleles. We have extended the proposed method to the analysis of genomic data for polyploid organisms, showing that calling accurate polyploid genotypes requires much higher coverage than diploid genotypes.

Weak link between dispersal and parasite community differentiation or immunogenetic divergence in two sympatric cichlid fishes

Geographical isolation, habitat variation and trophic specialization have contributed to a large extent to the astonishing diversity of cichlid fishes in the Great East African lakes. Because parasite communities often vary across space and environments, parasites can accompany and potentially enhance cichlid species diversification. However, host dispersal may reduce opportunities for parasite-driven evolution by homogenizing parasite communities and allele frequencies of immunity genes. To test for the relationships between parasite community variation, host dispersal and parasite-induced host evolution, we studied two sympatric cichlid species with contrasting dispersal capacities along the shores of southern Lake Tanganyika. Whereas the philopatric Tropheus moorii evolved into several genetically differentiated colour morphs, Simochromis diagramma is phenotypically rather uniform across its distribution range and shows only weak population structure. Populations of both species were infected with divergent parasite communities and harbour differentiated variant pools of an important set of immune genes, the major histocompatibility complex (MHC). The overall extent of geographical variation of parasites and MHC genes was similar between host species. This indicates that immunogenetic divergence among populations of Lake Tanganyika cichlids can occur even in species that are strongly dispersing. However, because this also includes species that are phenotypically uniform, parasite-induced evolution may not represent a key factor underlying species diversification in this system.

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.

Contrasting patterns of selection and drift between two categories of immune genes in prairie-chickens

Immune-receptor genes of the adaptive immune system, such as the major histocompatibility complex (MHC), are involved in recognizing specific pathogens and are known to have high rates of adaptive evolution, presumably as a consequence of rapid co-evolution between hosts and pathogens. In contrast, many 'mediating' genes of the immune system do not interact directly with specific pathogens and are involved in signalling (e.g. cytokines) or controlling immune cell growth. As a consequence, we might expect stronger selection at immune-receptor than mediating genes, but these two types of genes have not been compared directly in wild populations. Here, we tested the hypothesis that selection differs between MHC (class I and II) and mediating genes by comparing levels of population differentiation across the range of greater prairie-chickens (Tympanuchus cupido). As predicted, there was stronger population differentiation and isolation by distance at immune receptor (MHC) than at either mediating genes or neutral microsatellites, suggesting a stronger role of local adaptation at the MHC. In contrast, mediating genes displayed weaker differentiation between populations than neutral microsatellites, consistent with selection favouring similar alleles across populations for mediating genes. In addition to selection, drift also had a stronger effect on immune receptor (MHC) than mediating genes as indicated by the stronger decline of MHC variation in relation to population size. This is the first study in the wild to show that the effects of selection and drift on immune genes vary across populations depending on their functional role.

Emerging horizons for tick-borne pathogens: from the 'one pathogen-one disease' vision to the pathobiome paradigm

Ticks, as vectors of several notorious zoonotic pathogens, represent an important and increasing threat for human and animal health in Europe. Recent applications of new technology revealed the complexity of the tick microbiome, which may affect its vectorial capacity. Appreciation of these complex systems is expanding our understanding of tick-borne pathogens, leading us to evolve a more integrated view that embraces the 'pathobiome'; the pathogenic agent integrated within its abiotic and biotic environments. In this review, we will explore how this new vision will revolutionize our understanding of tick-borne diseases. We will discuss the implications in terms of future research approaches that will enable us to efficiently prevent and control the threat posed by ticks.

The importance of multiparasitism: examining the consequences of co-infections for human and animal health

Most parasites co-occur with other parasites, although the importance of such multiparasitism has only recently been recognised. Co-infections may result when hosts are independently infected by different parasites at the same time or when interactions among parasite species facilitate co-occurrence. Such interactions can have important repercussions on human or animal health because they can alter host susceptibility, infection duration, transmission risks, and clinical symptoms. These interactions may be synergistic or antagonistic and thus produce diverse effects in infected humans and animals. Interactions among parasites strongly influence parasite dynamics and therefore play a major role in structuring parasite populations (both within and among hosts) as well as host populations. However, several methodological challenges remain when it comes to detecting parasite interactions. The goal of this review is to summarise current knowledge on the causes and consequences of multiparasitism and to discuss the different methods and tools that researchers have developed to study the factors that lead to multiparasitism. It also identifies new research directions to pursue.

Among-lake reciprocal transplants induce convergent expression of immune genes in threespine stickleback

Geographic variation in parasite communities can drive evolutionary divergence in host immune genes. However, biotic and abiotic environmental variation can also induce plastic differences in immune function among populations. At present, there is little information concerning the relative magnitudes of heritable vs. induced immune divergence in natural populations. We examined immune gene expression profiles of threespine stickleback (Gasterosteus aculeatus) from six lakes on Vancouver Island, British Columbia. Parasite community composition differs between lake types (large or small, containing limnetic- or benthic-like stickleback) and between watersheds. We observed corresponding differences in immune gene expression profiles among wild-caught stickleback, using a set of seven immune genes representing distinct branches of the immune system. To evaluate the role of environmental effects on this differentiation, we experimentally transplanted wild-caught fish into cages in their native lake, or into a nearby foreign lake. Transplanted individuals' immune gene expression converged on patterns typical of their destination lake, deviating from their native expression profile. Transplant individuals' source population had a much smaller effect, suggesting relatively weak genetic underpinning of population differences in immunity, as viewed through gene expression. This strong environmental regulation of immune gene expression provides a counterpoint to the large emerging literature documenting microevolution and genetic diversification of immune function. Our findings illustrate the value of studying immunity in natural environmental settings where the immune system has evolved and actively functions.

Genetic structure in insular and mainland populations of house sparrows (Passer domesticus) and their hemosporidian parasites

Small and isolated populations usually exhibit low levels of genetic variability, and thus, they are expected to have a lower capacity to adapt to changes in environmental conditions, such as exposure to pathogens and parasites. Comparing the genetic variability of selectively neutral versus functional loci allows one to assess the evolutionary history of populations and their future evolutionary potential. The genes of the major histocompatibility complex (MHC) control immune recognition of parasites, and their unusually high diversity is genes which is likely driven by parasite-mediated balancing selection. Here, we examined diversity and differentiation of neutral microsatellite loci and functional MHC class I genes in house sparrows (Passer domesticus), living in six insular and six mainland populations, and we aimed to determine whether their diversity or differentiation correlates with the diversity and the prevalence of infection of hemosporidian parasites. We found that island bird populations tended to have lower neutral genetic variability, whereas MHC variability gene was similar between island and mainland populations. Similarly, island populations tended to show greater genetic differentiation than mainland populations, especially at microsatellite markers. The maintenance of MHC genetic diversity and its less marked structure in the island populations could be attributed to balancing-selection. The greater MHC differentiation among populations was negatively correlated with similarity in blood parasites (prevalence and diversity of parasite strains) between populations. Even at low prevalence and small geographical scale, haemosporidian parasites might contribute to structure the variability of immune genes among populations of hosts.

No evidence for the effect of MHC on male mating success in the brown bear

Mate choice is thought to contribute to the maintenance of the spectacularly high polymorphism of the Major Histocompatibility Complex (MHC) genes, along with balancing selection from parasites, but the relative contribution of the former mechanism is debated. Here, we investigated the association between male MHC genotype and mating success in the brown bear. We analysed fragments of sequences coding for the peptide-binding region of the highly polymorphic MHC class I and class II DRB genes, while controlling for genome-wide effects using a panel of 18 microsatellite markers. Male mating success did not depend on the number of alleles shared with the female or amino-acid distance between potential mates at either locus. Furthermore, we found no indication of female mating preferences for MHC similarity being contingent on the number of alleles the females carried. Finally, we found no significant association between the number of MHC alleles a male carried and his mating success. Thus, our results provided no support for the role of mate choice in shaping MHC polymorphism in the brown bear.

Major Histocompatibility Complex, demographic, and environmental predictors of antibody presence in a free-ranging mammal

Major Histocompatibility Complex (MHC) variability plays a key role in pathogen resistance, but its relative importance compared to environmental and demographic factors that also influence resistance is unknown. We analyzed the MHC II DRB exon 2 for 165 raccoons (Procyon lotor) in Missouri (USA). For each animal we also determined the presence of immunoglobulin G (IgG) and immunoglobulin M (IgM) antibodies to two highly virulent pathogens, canine distemper virus (CDV) and parvovirus. We investigated the role of MHC polymorphism and other demographic and environmental factors previously associated with predicting seroconversion. In addition, using an experimental approach, we studied the relative importance of resource availability and contact rates. We found important associations between IgG antibody presence and several MHC alleles and supertypes but not between IgM antibody presence and MHC. No effect of individual MHC diversity was found. For CDV, supertype S8, one allele within S8 (Prlo-DRB(∗)222), and a second allele (Prlo-DRB(∗)204) were positively associated with being IgG+, while supertype S4 and one allele within the supertype (Prlo-DRB(∗)210) were negatively associated with being IgG+. Age, year, and increased food availability were also positively associated with being IgG+, but allele Prlo-DRB(∗)222 was a stronger predictor. For parvovirus, only one MHC allele was negatively associated with being IgG+ and age and site were stronger predictors of seroconversion. Our results show that negative-frequency dependent selection is likely acting on the raccoon MHC and that while the role of MHC in relation to other factors depends on the pathogen of interest, it may be one of the most important factors predicting successful immune response.

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre-eminent system for the study of selective pressures that arise from host-pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population-genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.

Immunogenetic factors affecting susceptibility of humans and rodents to hantaviruses and the clinical course of hantaviral disease in humans

We reviewed the associations of immunity-related genes with susceptibility of humans and rodents to hantaviruses, and with severity of hantaviral diseases in humans. Several class I and class II HLA haplotypes were linked with severe or benign hantavirus infections, and these haplotypes varied among localities and hantaviruses. The polymorphism of other immunity-related genes including the C4A gene and a high-producing genotype of TNF gene associated with severe PUUV infection. Additional genes that may contribute to disease or to PUUV infection severity include non-carriage of the interleukin-1 receptor antagonist (IL-1RA) allele 2 and IL-1β (-511) allele 2, polymorphisms of plasminogen activator inhibitor (PAI-1) and platelet GP1a. In addition, immunogenetic studies have been conducted to identify mechanisms that could be linked with the persistence/clearance of hantaviruses in reservoirs. Persistence was associated during experimental infections with an upregulation of anti-inflammatory responses. Using natural rodent population samples, polymorphisms and/or expression levels of several genes have been analyzed. These genes were selected based on the literature of rodent or human/hantavirus interactions (some Mhc class II genes, Tnf promoter, and genes encoding the proteins TLR4, TLR7, Mx2 and β3 integrin). The comparison of genetic differentiation estimated between bank vole populations sampled over Europe, at neutral and candidate genes, has allowed to evidence signatures of selection for Tnf, Mx2 and the Drb Mhc class II genes. Altogether, these results corroborated the hypothesis of an evolution of tolerance strategies in rodents. We finally discuss the importance of these results from the medical and epidemiological perspectives.

To be or not to be associated: power study of four statistical modeling approaches to identify parasite associations in cross-sectional studies

A growing number of studies are reporting simultaneous infections by parasites in many different hosts. The detection of whether these parasites are significantly associated is important in medicine and epidemiology. Numerous approaches to detect associations are available, but only a few provide statistical tests. Furthermore, they generally test for an overall detection of association and do not identify which parasite is associated with which other one. Here, we developed a new approach, the association screening approach, to detect the overall and the detail of multi-parasite associations. We studied the power of this new approach and of three other known ones (i.e., the generalized chi-square, the network and the multinomial GLM approaches) to identify parasite associations either due to parasite interactions or to confounding factors. We applied these four approaches to detect associations within two populations of multi-infected hosts: (1) rodents infected with Bartonella sp., Babesia microti and Anaplasma phagocytophilum and (2) bovine population infected with Theileria sp. and Babesia sp. We found that the best power is obtained with the screening model and the generalized chi-square test. The differentiation between associations, which are due to confounding factors and parasite interactions was not possible. The screening approach significantly identified associations between Bartonella doshiae and B. microti, and between T. parva, T. mutans, and T. velifera. Thus, the screening approach was relevant to test the overall presence of parasite associations and identify the parasite combinations that are significantly over- or under-represented. Unraveling whether the associations are due to real biological interactions or confounding factors should be further investigated. Nevertheless, in the age of genomics and the advent of new technologies, it is a considerable asset to speed up researches focusing on the mechanisms driving interactions between parasites.

Parasite-mediated selection of major histocompatibility complex variability in wild brandt's voles (Lasiopodomys brandtii) from Inner Mongolia, China

Background

Genes of the major histocompatibility complex (MHC) exhibit high levels of variability, which is believed to have arisen through pathogen-mediated selection. We investigated the relationship between parasite load and genetic diversity at selectively neutral, non-coding markers (microsatellites) and adaptive genetic variation at a functionally important part of the MHC in six independent natural populations of Brandt’s voles (Lasiopodomys brandtii) from two regions of the Xilingol Grassland area of Inner Mongolia.

Results

Two-hundred and fifty-two voles were screened for gastrointestinal parasites, and were assessed for genetic variation. Parasite screening was done through non-invasive fecal egg counts, while allelic diversity was determined via single-stranded conformation polymorphism and DNA sequencing. We detected eight distinct helminth egg morphotypes. A total of 10 microsatellite loci were genotyped and 19 unique MHC class II B alleles were isolated. The rate of nonsynonymous substitutions (d_N) exceeded the rate of synonymous substitutions (d_S) at putative antigen binding sites of DRB. Neutral and adaptive genetic diversity differed between the six vole populations. To test the main pathogen-driven selection hypotheses for the maintenance of host MHC diversity and parasite species-specific co-evolutionary effects, multivariate approaches (generalized linear mixed models) were used to test for associations between the MHC class II DRB genotype and infections with nematodes. We found no evidence for heterozygote advantage, and overall heterozygosity was lower than expected in the MHC alleles. We identified an association between the parasite load and specific MHC alleles in the voles, and this pattern varied between geographic regions.

Conclusions

The results suggest that MHC variability in Brandt’s voles is maintained by rare allele advantage and fluctuating selection, but the data failed to show any heterozygote advantage effect. Our results add to a growing body of evidence showing that the mode and relative strength of pathogen-driven selection acting on MHC diversity varies within specific wild populations. In addition, our study contributes to the understanding of what maintains MHC diversity, of host-pathogen coevolution and of how genetic diversity is maintained in voles.

A probabilistic model in cross-sectional studies for identifying interactions between two persistent vector-borne pathogens in reservoir populations

BACKGROUND

In natural populations, individuals are infected more often by several pathogens than by just one. In such a context, pathogens can interact. This interaction could modify the probability of infection by subsequent pathogens. Identifying when pathogen associations correspond to biological interactions is a challenge in cross-sectional studies where the sequence of infection cannot be demonstrated.

METHODOLOGY/PRINCIPAL FINDINGS

Here we modelled the probability of an individual being infected by one and then another pathogen, using a probabilistic model and maximum likelihood statistics. Our model was developed to apply to cross-sectional data, vector-borne and persistent pathogens, and to take into account confounding factors. Our modelling approach was more powerful than the commonly used Chi-square test of independence. Our model was applied to detect potential interaction between Borrelia afzelii and Bartonella spp. that infected a bank vole population at 11% and 57% respectively. No interaction was identified.

CONCLUSIONS/SIGNIFICANCE

The modelling approach we proposed is powerful and can identify the direction of potential interaction. Such an approach can be adapted to other types of pathogens, such as non-persistents. The model can be used to identify when co-occurrence patterns correspond to pathogen interactions, which will contribute to understanding how organism communities are assembled and structured. In the long term, the model's capacity to better identify pathogen interactions will improve understanding of infectious risk.

Inference of selection based on temporal genetic differentiation in the study of highly polymorphic multigene families

The co-evolutionary arms race between host immune genes and parasite virulence genes is known as Red Queen dynamics. Temporal fluctuations in allele frequencies, or the ‘turnover’ of alleles at immune genes, are concordant with predictions of the Red Queen hypothesis. Such observations are often taken as evidence of host-parasite co-evolution. Here, we use computer simulations of the Major Histocompatibility Complex (MHC) of guppies (Poecilia reticulata) to study the turnover rate of alleles (temporal genetic differentiation, G'_ST). Temporal fluctuations in MHC allele frequencies can be ≥≤order of magnitude larger than changes observed at neutral loci. Although such large fluctuations in the MHC are consistent with Red Queen dynamics, simulations show that other demographic and population genetic processes can account for this observation, these include: (1) overdominant selection, (2) fluctuating population size within a metapopulation, and (3) the number of novel MHC alleles introduced by immigrants when there are multiple duplicated genes. Synergy between these forces combined with migration rate and the effective population size can drive the rapid turnover in MHC alleles. We posit that rapid allelic turnover is an inherent property of highly polymorphic multigene families and that it cannot be taken as evidence of Red Queen dynamics. Furthermore, combining temporal samples in spatial F_ST outlier analysis may obscure the signal of selection.

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.

Diversity at the major histocompatibility complex Class II in the platypus, Ornithorhynchus anatinus

The platypus (Ornithorhynchus anatinus) is the sole survivor of a previously widely distributed and diverse lineage of ornithorhynchid monotremes. Its dependence on healthy water systems imposes an inherent sensitivity to habitat degradation and climate change. Here, we compare genetic diversity at the major histocompatibility complex (MHC) Class II-DZB gene and 3 MHC-associated microsatellite markers with diversity at 6 neutral microsatellite markers in 70 platypuses from across their range, including the mainland of Australia and the isolated populations of Tasmania, King Island, and Kangaroo Island. Overall, high DZB diversity was observed in the platypus, with 57 DZB β1 alleles characterized. Significant positive selection was detected within the DZB peptide-binding region, promoting variation in this domain. Low levels of genetic diversity were detected at all markers in the 2 island populations, King Island (endemic) and Kangaroo Island (introduced), with the King Island platypuses monomorphic at the DZB locus. Loss of MHC diversity on King Island is of concern, as the population may have compromised immunological fitness and reduced ability to resist changing environmental conditions.

Insights into the complex associations between MHC class II DRB polymorphism and multiple gastrointestinal parasite infestations in the striped mouse

Differences in host susceptibility to different parasite types are largely based on the degree of matching between immune genes and parasite antigens. Specifically the variable genes of the major histocompatibility complex (MHC) play a major role in the defence of parasites. However, underlying genetic mechanisms in wild populations are still not well understood because there is a lack of studies which deal with multiple parasite infections and their competition within. To gain insights into these complex associations, we implemented the full record of gastrointestinal nematodes from 439 genotyped individuals of the striped mouse, Rhabdomys pumilio. We used two different multivariate approaches to test for associations between MHC class II DRB genotype and multiple nematodes with regard to the main pathogen-driven selection hypotheses maintaining MHC diversity and parasite species-specific co-evolutionary effects. The former includes investigations of a 'heterozygote advantage', or its specific form a 'divergent-allele advantage' caused by highly dissimilar alleles as well as possible effects of specific MHC-alleles selected by a 'rare allele advantage' (= negative 'frequency-dependent selection'). A combination of generalized linear mixed models (GLMMs) and co-inertia (COIA) analyses made it possible to consider multiple parasite species despite the risk of type I errors on the population and on the individual level. We could not find any evidence for a 'heterozygote' advantage but support for 'divergent-allele' advantage and infection intensity. In addition, both approaches demonstrated high concordance of positive as well as negative associations between specific MHC alleles and certain parasite species. Furthermore, certain MHC alleles were associated with more than one parasite species, suggesting a many-to-many gene-parasite co-evolution. The most frequent allele Rhpu-DRB*38 revealed a pleiotropic effect, involving three nematode species. Our study demonstrates the co-existence of specialist and generalist MHC alleles in terms of parasite detection which may be an important feature in the maintenance of MHC polymorphism.

The impact of multiple infections on wild animal hosts: a review

Field parasitological studies consistently demonstrate the reality of polyparasitism in natural systems. However, only recently, studies from ecological and evolutionary fields have emphasised a broad spectrum of potential multiple infections-related impacts. The main goal of our review is to reunify the different approaches on the impacts of polyparasitism, not only from laboratory or human medical studies but also from field or theoretical studies. We put forward that ecological and epidemiological determinants to explain the level of polyparasitism, which regularly affects not only host body condition, survival or reproduction but also host metabolism, genetics or immune investment. Despite inherent limitations of all these studies, multiple infections should be considered more systematically in wildlife to better appreciate the importance of parasite diversity in wildlife, cumulative effects of parasitism on the ecology and evolution of their hosts.

Varying disease-mediated selection at different life-history stages of Atlantic salmon in fresh water

Laboratory studies on associations between disease resistance and susceptibility and major histocompatibility (MH) genes in Atlantic salmon Salmo salar have shown the importance of immunogenetics in understanding the capacity of populations to fight specific diseases. However, the occurrence and virulence of pathogens may vary spatially and temporally in the wild, making it more complicated to predict the overall effect that MH genes exert on fitness of natural populations and over several life-history stages. Here we show that MH variability is a significant determinant of salmon survival in fresh water, by comparing observed and expected genotype frequencies at MH and control microsatellite loci at parr and migrant stages in the wild. We found that additive allelic effects at immunogenetic loci were more likely to determine survival than dominance deviation, and that selection on certain MH alleles varied with life stage, possibly owing to varying pathogen prevalence and/or virulence over time. Our results highlight the importance of preserving genetic diversity (particularly at MH loci) in wild populations, so that they have the best chance of adapting to new and increased disease challenges as a result of projected climate warming and increasing aquaculture.

MHC variation, multiple simultaneous infections and physiological condition in the subterranean rodent Ctenomys talarum

Parasites and pathogens can play a significant role in shaping the genetic diversity of host populations, particularly at genes associated with host immune response. To explore this relationship in a natural population of vertebrates, we characterized Major Histocompatibility Complex (MHC) variation in the subterranean rodent Ctenomys talarum (the talas tuco-tuco) as a function of parasite load and ability to mount an adaptive immune response against a novel antigen. Specifically, we quantified genotypic diversity at the MHC class II DRB locus in relation to (1) natural variation in infection by multiple genera of parasites (potential agents of selection on MHC genes) and (2) antibody production in response to injection with sheep red blood cells (a measure of immunocompetence). Data were analyzed using co-inertia multivariate statistics, with epidemiological proxies for individual condition (hematocrit, leukocyte profile, body weight) and risk of parasite exposure (season of capture, sex). A significant excess of DRB heterozygotes was evident in the study population. Co-inertia analyses revealed significant associations between specific DRB alleles and both parasite load and intensity of humoral immune response against sheep red blood cells. The presence of specific DRB aminoacid sequences appeared to be more strongly associated with parasite load and response to a novel antigen than was heterozygosity at the DRB locus. These data suggest a role for parasite-driven balancing selection in maintaining MHC variation in natural populations of C. talarum. At the same time, these findings underscore the importance of using diverse parameters to study interactions among physiological conditions, immunocompetence, and MHC diversity in free-living animals that are confronted with multiple simultaneous immune challenges.

Major histocompatibility complex polymorphism: dynamics and consequences of parasite-mediated local adaptation in fishes

Parasitism is a common form of life and represents a strong selective pressure for host organisms. In response to this evolutionary pressure, vertebrates have developed genetically coded defences such as the major histocompatibility complex (MHC). Mechanisms of parasite-mediated selection not only maintain outstanding polymorphism in these genes but have also been proposed to further promote host population divergence and ultimately speciation because it can drive evolution of local adaptation in which MHC genes play a crucial role. This review first highlights the dynamics and complexity of parasite-mediated selection in natural systems, which not only depends on dominating parasite strategies and on the taxonomic diversity of the parasite community but also includes the differences in parasite communities between habitats and niches, creating divergent selection on locally adapted populations. Then the different ways in which MHC genes potentially allow vertebrates to respond to these dynamics and to adapt locally are outlined. Finally, it is proposed that varying selection strength in time and space may lead to variation in the strength of precopulatory reproductive isolation which has evolved to maintain local adaptation.

Parasite-mediated evolution of the functional part of the MHC in primates

The major histocompatibility complex (MHC) is a key model of genetic polymorphism, but the mechanisms underlying its extreme variability are debated. Most hypotheses for MHC diversity focus on pathogen-driven selection and predict that MHC polymorphism evolves under the pressure of a diverse parasite fauna. Several studies reported that certain alleles offer protection against certain parasites, yet it remains unclear whether variation in parasite pressure more generally covaries with allelic diversity and rates of molecular evolution of MHC across species. We tested this prediction in a comparative study of 41 primate species. We characterized polymorphism of the exon 2 of DRB region of the MHC class II. Our phylogenetic analyses controlled for the potential effects of neutral mutation rate, population size, geographic origin and body mass and revealed that nematode species richness associates positively with nonsynonymous nucleotide substitution rate at the functional part of the molecule. We failed to find evidence for allelic diversity being strongly related to parasite species richness. Continental distribution was a strong predictor of both allelic diversity and substitution rate, with higher values in Malagasy and Neotropical primates. These results indicate that parasite pressure can influence the different estimates of MHC polymorphism, whereas geography plays an independent role in the natural history of MHC.

Genetic structure and contrasting selection pattern at two major histocompatibility complex genes in wild house mouse populations

The mammalian major histocompatibility complex (MHC) is a tightly linked cluster of immune genes, and is often thought of as inherited as a unit. This has led to the hope that studying a single MHC gene will reveal patterns of evolution representative of the MHC as a whole. In this study we analyse a 1000-km transect of MHC variation traversing the European house mouse hybrid zone to compare signals of selection and patterns of diversification at two closely linked MHC class II genes, H-2Aa and H-2Eb. We show that although they are 0.01 cM apart (that is, recombination is expected only once in 10 000 meioses), disparate evolutionary patterns were detected. H-2Aa shows higher allelic polymorphism, faster allelic turnover due to higher mutation rates, stronger positive selection at antigen-binding sites and higher population structuring than H-2Eb. H-2Eb alleles are maintained in the gene pool for longer, including over separation of the subspecies, some H-2Eb alleles are positively and others negatively selected and some of the alleles are not expressed. We conclude that studies on MHC genes in wild-living vertebrates can give substantially different results depending on the MHC gene examined and that the level of polymorphism in a related species is a poor criterion for gene choice.

Major histocompatibility complex variation and age-specific endoparasite load in subadult European rabbits

Genes of the major histocompatibility complex (MHC) play a fundamental role in the vertebrate immune response and are amongst the most polymorphic genes in vertebrate genomes. It is generally agreed that the highly polymorphic nature of the MHC is maintained through host-parasite co-evolution. Two nonexclusive mechanisms of selection are supposed to act on MHC genes: superiority of MHC heterozygous individuals (overdominance) and an advantage for rare MHC alleles. However, the precise mechanisms and their relative importance are still unknown. Here, we examined MHC dependent parasite load in European rabbits (Oryctolagus cuniculus) from a distinct population with low MHC diversity (three alleles, six genotypes). Using a multivariate approach, we tested for associations of individual MHC class II DRB constitution and the rabbits' intestinal burden with nematodes and coccidia. Rabbits having a particular allele showed lower infestations with hepatic coccidia (E. stiedai). However, a comparison of all six genotypes in the population revealed that carriers of this allele only benefit when they are heterozygous, and furthermore, MHC heterozygosity in general did not affect individual parasite load. In conclusion, this study suggests an immunogenetic basis of European rabbit resistance to hepatic coccidiosis, which can strongly limit survival to maturity in this species. Our study gives a complex picture of MHC-parasite correlations, unveiling the limits of the classical hypotheses of how MHC polymorphism is maintained in natural systems.

Functional associations of similar MHC alleles and shared parasite species in two sympatric lemurs

Several recent studies of animals in their natural surroundings found evidence for effects of certain major histocompatibility complex (MHC) immune gene alleles on the parasite load. However, in multi-infected individuals the particular selection pressure exerted by specific parasites has rarely been explored. In this study we took advantage of the parasitological and genetic data of two previously investigated Malagasy lemur species (Cheirogaleus medius and Microcebus murinus). We investigated whether the two sympatric and ecologically similar primates are infected by similar parasite species and explored if certain parasites are associated with particular MHC alleles. Our study revealed that most of the parasite egg morphotypes were found in both hosts. In each lemur species we identified one MHC allele which was positively associated with Ascaris-infection. Interestingly, these MHC alleles were very similar to each other but differed from all other investigated MHC alleles in an amino acid substitution in a putative functional important antigen binding site. Thus, our study gives first intriguing evidence for a direct connection between certain antigen binding sites of MHC molecules with a particular parasite in two wild primate populations. This may indicate that indeed certain parasites exert direct selective pressure on the MHC of wild living hosts.

Negative relationships between cellular immune response, Mhc class II heterozygosity and secondary sexual trait in the montane water vole

Heterogeneities in immune responsiveness may affect key epidemiological parameters and the dynamics of pathogens. The roles of immunogenetics in these variations remain poorly explored. We analysed the influence of Major histocompatibility complex (Mhc) genes and epigamic traits on the response to phytohaemagglutinin in males from cyclic populations of the montane water vole (Arvicola scherman). Besides, we tested the relevance of lateral scent glands as honest signals of male quality. Our results did not corroborate neither the hypotheses of genome-wide heterozygosity-fitness correlation nor the Mhc heterozygote advantage. We found a negative relationship between Mhc hetetozygosity and response to phytohaemagglutinin, mediated by a specific Mhc homozygous genotype. Our results therefore support the hypothesis of the Arte-Dqa-05 homozygous genotype being a ‘good’ Mhc variant in terms of immunogenetic quality. The development of the scent glands seems to be an honest signal for mate choice as it is negatively correlated with helminth load. The ‘good gene’ hypothesis was not validated as Arte-Dqa-05 homozygous males did not exhibit larger glands. Besides, the negative relationship observed between the size of these glands and the response to phytohaemagglutinin, mainly for Mhc homozygotes, corroborates the immunocompetence handicap hypothesis. The Mhc variants associated with larger glands remain yet to be determined.

How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings

Major histocompatibility complex (MHC) genes have been put forward as a model for studying how genetic diversity is maintained in wild populations. Pathogen-mediated selection (PMS) is believed to generate the extraordinary levels of MHC diversity observed. However, establishing the relative importance of the three proposed mechanisms of PMS (heterozygote advantage, rare-allele advantage and fluctuating selection) has proved extremely difficult. Studies have attempted to differentiate between mechanisms of PMS using two approaches: (i) comparing MHC diversity with that expected under neutrality and (ii) relating MHC diversity to pathogen regime. Here, we show that in many cases the same predictions arise from the different mechanisms under these approaches, and that most studies that have inferred one mechanism of selection have not fully considered the alternative explanations. We argue that, while it may be possible to demonstrate that particular mechanisms of PMS are occurring, resolving their relative importance within a system is probably impossible. A more realistic target is to continue to demonstrate when and where the different mechanisms of PMS occur, with the aim of determining their relative importance across systems. We put forward what we believe to be the most promising approaches that will allow us to progress towards achieving this.

Major histocompatibility complex heterozygote advantage and widespread bacterial infections in populations of Chinook salmon (Oncorhynchus tshawytscha)

Despite growing evidence for parasite-mediated selection on the vertebrate major histocompatibility complex (MHC), little is known about variation in the bacterial parasite community within and among host populations or its influence on MHC evolution. In this study, we characterize variation in the parasitic bacterial community associated with Chinook salmon (Oncorhynchus tshawytscha) fry in five populations in British Columbia (BC), Canada across 2 years, and examine whether bacterial infections are a potential source of selection on the MHC. We found an unprecedented diversity of bacteria infecting fry with a total of 55 unique bacteria identified. Bacterial infection rates varied from 9% to 29% among populations and there was a significant isolation by distance relationship in bacterial community phylogenetic similarity across the populations. Spatial variation in the frequency of infections and in the phylogenetic similarity of bacterial communities may result in differential parasite-mediated selection at the MHC across populations. Across all populations, we found evidence of a heterozygote advantage at the MHC class II, which may be a source of balancing selection on this locus. Interestingly, a co-inertia analysis indicated only susceptibility associations between a few of the MHC class I and II alleles and specific bacterial parasites; there was no evidence that any of the alleles provided resistance to the bacteria. Our results reveal a complex bacterial community infecting populations of a fish and underscore the importance of considering the role of multiple pathogens in the evolution of host adaptations.