Major histocompatibility complex (MHC) heterozygote superiority to natural multi-parasite infections in the water vole (Arvicola terrestris)

Heterozygote advantage can explain the extraordinary diversity of immune genes

The majority of highly polymorphic genes are related to immune functions and with over 100 alleles within a population, genes of the major histocompatibility complex (MHC) are the most polymorphic loci in vertebrates. How such extraordinary polymorphism arose and is maintained is controversial. One possibility is heterozygote advantage (HA), which can in principle maintain any number of alleles, but biologically explicit models based on this mechanism have so far failed to reliably predict the coexistence of significantly more than 10 alleles. We here present an eco-evolutionary model showing that evolution can result in the emergence and maintenance of more than 100 alleles under HA if the following two assumptions are fulfilled: first, pathogens are lethal in the absence of an appropriate immune defence; second, the effect of pathogens depends on host condition, with hosts in poorer condition being affected more strongly. Thus, our results show that HA can be a more potent force in explaining the extraordinary polymorphism found at MHC loci than currently recognised.

Climate heterogeneity, season variation, and sexual dimorphism modulate the association between MHC II diversity and parasite variation in striped hamster

Parasite-mediated selection is widely believed to play a crucial role in maintaining the diversity of the major histocompatibility complex (MHC) genes, which is thought to be maintained through heterozygote advantage, rare-allele advantage, and fluctuating selection. However, the relationship between parasite pressure and MHC diversity has yielded inconsistent findings. These inconsistencies may arise from the influence of environmental factors and individual variations in traits on host-parasite interactions. To address these issues, our study extensively investigated populations of striped hamsters inhabiting regions characterized by environmental heterogeneity. The primary objective was to examine the universality of parasite-mediated selection mechanisms. Our observations revealed the presence of multiple parasite infections, accompanied by spatial and temporal variations in parasite communities and infection patterns among individual hamsters. Specifically, the temperature was found to influence all four parasite indices, while the presence of gamasid mites and parasite richness decreased with increasing precipitation. We also noted significant seasonal variation in parasite dynamics. Moreover, a significant sexual dimorphism was observed with males exhibiting a considerably higher parasite burden compared to their female counterparts. Lastly, we identified the maintenance of MHC polymorphism in striped hamsters as being driven by the heterozygote advantage and fluctuating selection mechanisms. This study underscores the significance of ecological processes in comprehending host-parasite systems and highlights the necessity of considering environmental factors and individual traits when elucidating the mechanisms underlying MHC diversity mediated by parasites.

MHC Class II Supertypes Affect Survival and Lifetime Reproductive Success in a Migratory Songbird

The major histocompatibility complex (MHC) plays a critical role in the immune response against pathogens. Its high polymorphism is thought to be mainly the consequence of host-pathogen co-evolution, but elucidating the mechanism(s) driving MHC evolution remains challenging for natural populations. We investigated the diversity of MHC class II genes in a wild population of pied flycatchers Ficedula hypoleuca and tested its associations with two key components of individual fitness: lifetime reproductive success and survival. Among 180 breeding adults in our study population, we found 182 unique MHC class II exon 2 alleles. The alleles showed a strong signal of positive selection and grouped into nine functional supertypes based on physicochemical properties at the inferred antigen-binding sites. Three supertypes were found in > 98% of the sampled individuals, indicating that they are nearly fixed in the population. We found no rare supertypes in the population, as all supertypes were present in > 70% of individuals. Three supertypes were related to different components of individual fitness: two were associated with lower offspring production over time, while the third was positively associated with survival. Overall, the substantial allelic and functional diversity and the relationship between specific supertypes and fitness are in accordance with the notion that balancing selection maintains MHC class II diversity in the study population, possibly with fluctuating selection as the underlying mechanism. The absence of rare supertypes in the population suggests that the balancing selection is not driven by rare-allele advantage.

Putative SNPs in Ovar-DRB1 and GALNTL6 Genes Conferring Susceptibility to Natural Infection of Haemonchus Contortus in Southern Indian Sheep

AIM

To explore associations between phenotypic traits and polymorphisms in the DRB1 and GALNT6 gene in Nellore, Deccani and Kenguri sheep naturally infected with Haemonchus contortus.

MATERIALS AND METHODS

Blood and faecal samples were collected to evaluate fecal worm egg counts (FEC), packed cell volume (PCV), hemoglobin (Hb), eosinophilia and for DNA isolation.

RESULTS

Animals were grouped into susceptible and resistant groups based on EPG counts. FEC and circulating eosinophilia were higher in a susceptible group. Log FEC was negatively correlated (P < 0.01) with PCV, and Hb estimates. The second exon of DRB1 and intron variant of GALNTL6 genes were amplified from DNA samples of resistant and susceptible sheep. Characterization of Ovar-DRB1 amplicon by RFLP revealed two genotypes ('bb' and 'ab'). The genotype frequencies differed significantly between both groups (P < 0.05). The 'bb' genotypes had higher (P < 0.05) log FEC value than 'ab' genotypes and 'b' allele was linked with susceptibility to haemonchosis in sheep. The mean FEC of Nellore sheep was high indicating susceptibility of the breed and also in which the frequency of 'b' allele was more compared to the other two breeds. OVAR-DRB1 genotypes associated with FEC did not affect PCV and Hb. PCR-RFLP assay developed to determine the genotypes with respect to SNP rs424521894 of GALNTL6 revealed monomorphic nature at the locus in the breeds studied.

CONCLUSION

MHC polymorphism could be used as a genetic marker for the selection of sheep resistant to H. contortus. However, a more intensive study, involving controlled infections and other GALNTL6 SNPs may be enforced to make any decisive assertion.

Immunogenetic-pathogen networks shrink in Tome's spiny rat, a generalist rodent inhabiting disturbed landscapes

Anthropogenic disturbance may increase the emergence of zoonoses. Especially generalists that cope with disturbance and live in close contact with humans and livestock may become reservoirs of zoonotic pathogens. Yet, whether anthropogenic disturbance modifies host-pathogen co-evolutionary relationships in generalists is unknown. We assessed pathogen diversity, neutral genome-wide diversity (SNPs) and adaptive MHC class II diversity in a rodent generalist inhabiting three lowland rainforest landscapes with varying anthropogenic disturbance, and determined which MHC alleles co-occurred more frequently with 13 gastrointestinal nematodes, blood trypanosomes, and four viruses. Pathogen-specific selection pressures varied between landscapes. Genome-wide diversity declined with the degree of disturbance, while MHC diversity was only reduced in the most disturbed landscape. Furthermore, pristine forest landscapes had more functional important MHC-pathogen associations when compared to disturbed forests. We show co-evolutionary links between host and pathogens impoverished in human-disturbed landscapes. This underscores that parasite-mediated selection might change even in generalist species following human disturbance which in turn may facilitate host switching and the emergence of zoonoses.

Immunogenetics, sylvatic plague and its vectors: insights from the pathogen reservoir Mastomys natalensis in Tanzania

Yersinia pestis is a historically important vector-borne pathogen causing plague in humans and other mammals. Contemporary zoonotic infections with Y. pestis still occur in sub-Saharan Africa, including Tanzania and Madagascar, but receive relatively little attention. Thus, the role of wildlife reservoirs in maintaining sylvatic plague and spillover risks to humans is largely unknown. The multimammate rodent Mastomys natalensis is the most abundant and widespread rodent in peri-domestic areas in Tanzania, where it plays a major role as a Y. pestis reservoir in endemic foci. Yet, how M. natalensis' immunogenetics contributes to the maintenance of plague has not been investigated to date. Here, we surveyed wild M. natalensis for Y. pestis vectors, i.e., fleas, and tested for the presence of antibodies against Y. pestis using enzyme-linked immunosorbent assays (ELISA) in areas known to be endemic or without previous records of Y. pestis in Tanzania. We characterized the allelic and functional (i.e., supertype) diversity of the major histocompatibility complex (MHC class II) of M. natalensis and investigated links to Y. pestis vectors and infections. We detected antibodies against Y. pestis in rodents inhabiting both endemic areas and areas considered non-endemic. Of the 111 nucleotide MHC alleles, only DRB*016 was associated with an increased infestation with the flea Xenopsylla. Surprisingly, we found no link between MHC alleles or supertypes and antibodies of Y. pestis. Our findings hint, however, at local adaptations towards Y. pestis vectors, an observation that more exhaustive sampling could unwind in the future.

Population-wide genetic analysis of Ovar-DQA1 and DQA2 loci across sheep breeds in India revealed their evolutionary importance and fitness of sheep in a tropical climate

Genetic variability at the major histocompatibility complex (MHC) is important in any species due to significant role played by MHC for antigen presentation. DQA locus has not been studied for its genetic variability across sheep population in India. In the present study, MHC of sheep at DQA1 and DQA2 loci were evaluated across 17 Indian sheep breeds. Results revealed high degree of heterozygosity (10.34% to 100% for DQA1 and 37.39 to 100% for DQA2). 18 DQA1 alleles and 22 DQA2 alleles were isolated in different breeds. Nucleotide content for DQA region revealed richness of AT content (54.85% for DQA1 and 53.89% for DQA2). DQA1 and DQA2 sequences clustered independently. We could see evidence of divergence of DQA as DQA1 and DQA2 across sheep breeds. Wu-Kabat variability index revealed vast genetic variation across DQA1 and DQA2, specifically at peptide binding sites (PBS) that consisted 21 residues for DQA1 and 17 residues for DQA2. Evolutionary analysis revealed the presence of positive and balancing selection for DQA1 locus, however DQA2 was under purifying selection across sheep breeds. Higher heterozygosity and large diversity at both loci especially at PBS indicated the fitness of the sheep population for evading pathogens and adapt to the harsh tropical climate.

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.

Genomic consequences of a century of inbreeding and isolation in the Danish wild boar population

Demographic events such as series of bottlenecks impact the genetic variation and adaptive potential of populations. European megafauna, such as wild boars (Sus scrofa), have experienced severe climatic and size fluctuations that have shaped their genetic variation. Habitat fragmentation and human-mediated translocations have further contributed to the complex demographic history of European wild boar. Danish wild boars represent an extreme case of a small and isolated population founded by four wild boars from Germany. Here, we explore the genetic composition of the Danish wild boar population in Klelund. We genotyped all 21 Danish wild boars that were recently transferred from the source population in Lille Vildmose into the Klelund Plantation to establish a novel wild boar population. We compared the Danish wild boars with high-density single-nucleotide polymorphism genotypes from a comprehensive reference set of 1263 wild and domesticated pigs, including 11 individuals from Ulm, one of two presumed founder locations in Germany. Our findings support the European wild background of the Danish population, and no traces of gene flow with wild or domesticated pigs were found. The narrow genetic origin of the Danish wild boars is illustrated by extremely long and frequent runs of homozygous stretches in their genomes, indicative of recent inbreeding. This study provides the first insights into one of the most inbred wild boar populations globally established a century ago from a narrow base of only four founders.

Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics

Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes - arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.

Selection and demography drive range-wide patterns of MHC-DRB variation in mule deer

BACKGROUND

Standing genetic variation is important especially in immune response-related genes because of threats to wild populations like the emergence of novel pathogens. Genetic variation at the major histocompatibility complex (MHC), which is crucial in activating the adaptive immune response, is influenced by both natural selection and historical population demography, and their relative roles can be difficult to disentangle. To provide insight into the influences of natural selection and demography on MHC evolution in large populations, we analyzed geographic patterns of variation at the MHC class II DRB exon 2 locus in mule deer (Odocoileus hemionus) using sequence data collected across their entire broad range.

RESULTS

We identified 31 new MHC-DRB alleles which were phylogenetically similar to other cervid MHC alleles, and one allele that was shared with white-tailed deer (Odocoileus virginianus). We found evidence for selection on the MHC including high dN/dS ratios, positive neutrality tests, deviations from Hardy-Weinberg Equilibrium (HWE) and a stronger pattern of isolation-by-distance (IBD) than expected under neutrality. Historical demography also shaped variation at the MHC, as indicated by similar spatial patterns of variation between MHC and microsatellite loci and a lack of association between genetic variation at either locus type and environmental variables.

CONCLUSIONS

Our results show that both natural selection and historical demography are important drivers in the evolution of the MHC in mule deer and work together to shape functional variation and the evolution of the adaptive immune response in large, well-connected populations.

Neutral and Selective Processes Shape MHC Diversity in Roe Deer in Slovenia

Simple Summary

Disease prevention and appropriate wildlife management are among the major challenges in wildlife conservation. In the present study, we made a first assessment of the variability of major histocompatibility complex (MHC) genes in roe deer in Slovenia and evaluated local population adaptation by comparing MHC variability with neutral microsatellites. We discovered three new MHC DRB exon 2 alleles in addition to seven previously described in the literature. Moreover, we found evidence of historical positive selection, as selection analysis indicated that approx. 10% of the encoded amino acids were subjected to episodic positive selection. This study provides the basis for further research on immunogenetic variation in roe deer and highlights opportunities to incorporate genetic data into science-based population management.

Abstract

Disease control and containment in free-ranging populations is one of the greatest challenges in wildlife management. Despite the importance of major histocompatibility complex (MHC) genes for immune response, an assessment of the diversity and occurrence of these genes is still rare in European roe deer, the most abundant and widespread large mammal in Europe. Therefore, we examined immunogenetic variation in roe deer in Slovenia to identify species adaptation by comparing the genetic diversity of the MHC genes with the data on neutral microsatellites. We found ten MHC DRB alleles, three of which are novel. Evidence for historical positive selection on the MHC was found using the maximum likelihood codon method. Patterns of MHC allelic distribution were not congruent with neutral population genetic findings. The lack of population genetic differentiation in MHC genes compared to existing structure in neutral markers suggests that MHC polymorphism was influenced primarily by balancing selection and, to a lesser extent, by neutral processes such as genetic drift, with no clear evidence of local adaptation. Selection analyses indicated that approx. 10% of amino acids encoded under episodic positive selection. This study represents one of the first steps towards establishing an immunogenetic map of roe deer populations across Europe, aiming to better support science-based management of this important game species.

Genetic diversity and disease: The past, present, and future of an old idea

Why do infectious diseases erupt in some host populations and not others? This question has spawned independent fields of research in evolution, ecology, public health, agriculture, and conservation. In the search for environmental and genetic factors that predict variation in parasitism, one hypothesis stands out for its generality and longevity: genetically homogeneous host populations are more likely to experience severe parasitism than genetically diverse populations. In this perspective piece, I draw on overlapping ideas from evolutionary biology, agriculture, and conservation to capture the far-reaching implications of the link between genetic diversity and disease. I first summarize the development of this hypothesis and the results of experimental tests. Given the convincing support for the protective effect of genetic diversity, I then address the following questions: (1) Where has this idea been put to use, in a basic and applied sense, and how can we better use genetic diversity to limit disease spread? (2) What new hypotheses does the established disease-diversity relationship compel us to test? I conclude that monitoring, preserving, and augmenting genetic diversity is one of our most promising evolutionarily informed strategies for buffering wild, domesticated, and human populations against future outbreaks.

Contemporary selection on MHC genes in a free-living ruminant population

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate contemporary natural selection is acting on MHC class II genes in Soay sheep and the mode of selection on specific fitness components can be different mode from selection on total fitness.

Heterozygosity of the major histocompatibility complex predicts later self-reported pubertal maturation in men

Individual variation in the age of pubertal onset is linked to physical and mental health, yet the factors underlying this variation are poorly understood. Life history theory predicts that individuals at higher risk of mortality due to extrinsic causes such as infectious disease should sexually mature and reproduce earlier, whereas those at lower risk can delay puberty and continue to invest resources in somatic growth. We examined relationships between a genetic predictor of infectious disease resistance, heterozygosity of the major histocompatibility complex (MHC), referred to as the human leukocyte antigen (HLA) gene in humans, and self-reported pubertal timing. In a combined sample of men from Canada (n = 137) and the United States (n = 43), MHC heterozygosity predicted later self-reported pubertal development. These findings suggest a genetic trade-off between immunocompetence and sexual maturation in human males.

Predicting insect outbreaks using machine learning: A mountain pine beetle case study

Planning forest management relies on predicting insect outbreaks such as mountain pine beetle, particularly in the intermediate-term future, e.g., 5-year. Machine-learning algorithms are potential solutions to this challenging problem due to their many successes across a variety of prediction tasks. However, there are many subtle challenges in applying them: identifying the best learning models and the best subset of available covariates (including time lags) and properly evaluating the models to avoid misleading performance-measures. We systematically address these issues in predicting the chance of a mountain pine beetle outbreak in the Cypress Hills area and seek models with the best performance at predicting future 1-, 3-, 5- and 7-year infestations. We train nine machine-learning models, including two generalized boosted regression trees (GBM) that predict future 1- and 3-year infestations with 92% and 88% AUC, and two novel mixed models that predict future 5- and 7-year infestations with 86% and 84% AUC, respectively. We also consider forming the train and test datasets by splitting the original dataset randomly rather than using the appropriate year-based approach and show that this may obtain models that score high on the test dataset but low in practice, resulting in inaccurate performance evaluations. For example, a k-nearest neighbor model with the actual performance of 68% AUC, scores the misleadingly high 78% on a test dataset obtained from a random split, but the more accurate 66% on a year-based split. We then investigate how the prediction accuracy varies with respect to the provided history length of the covariates and find that neural network and naive Bayes, predict more accurately as history-length increases, particularly for future 1- and 3-year predictions, and roughly the same holds with GBM. Our approach is applicable to other invasive species. The resulting predictors can be used in planning forest and pest management and planning sampling locations in field studies.

Functional immunogenetic variation, rather than local adaptation, predicts ectoparasite infection intensity in a model fish species

Natural host populations differ in their susceptibility to infection by parasites, and these intrapopulation differences are still an incompletely understood component of host‐parasite dynamics. In this study, we used controlled infection experiments with wild‐caught guppies (Poecilia reticulata) and their ectoparasite Gyrodactylus turnbulli to investigate the roles of local adaptation and host genetic composition (immunogenetic and neutral) in explaining differences in susceptibility to infection. We found differences between our four study host populations that were consistent between two parasite source populations, with no indication of local adaptation by either host or parasite at two tested spatial scales. Greater values of host population genetic variability metrics broadly aligned with lower population mean infection intensity, with the best alignments associated with major histocompatibility complex (MHC) “supertypes”. Controlling for intrapopulation differences and potential inbreeding variance, we found a significant negative relationship between individual‐level functional MHC variability and infection: fish carrying more MHC supertypes experienced infections of lower severity, with limited evidence for supertype‐specific effects. We conclude that population‐level differences in host infection susceptibility probably reflect variation in parasite selective pressure and/or host evolutionary potential, underpinned by functional immunogenetic variation.

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.

Positive selection on the MHC class II DLA-DQA1 gene in golden jackals (Canis aureus) from their recent expansion range in Europe and its effect on their body mass index

Background

In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data.

Results

Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy–Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001.

Conclusions

Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01856-z.

Trade-offs in expressed major histocompatibility complex diversity seen on a macroevolutionary scale among songbirds

To survive organisms must defend themselves against pathogens. Classical Major Histocompatibility Complex (MHC) genes play a key role in pathogen defense by encoding molecules involved in pathogen recognition. MHC gene diversity influences the variety of pathogens individuals can recognize and respond to and has consequently been a popular genetic marker for disease resistance in ecology and evolution. However, MHC diversity is predominantly estimated using genomic DNA (gDNA) with little knowledge of expressed diversity. This limits our ability to interpret the adaptive significance of variation in MHC diversity, especially in species with very many MHC genes such as songbirds. Here, we address this issue using phylogenetic comparative analyses of the number of MHC class I alleles (MHC-I diversity) in gDNA and complementary DNA (cDNA), that is, expressed alleles, across 13 songbird species. We propose three theoretical relationships that could be expected between genomic and expressed MHC-I diversity on a macroevolutionary scale and test which of these are best supported. In doing so, we show that significantly fewer MHC-I alleles than the number available are expressed, suggesting that optimal MHC-I diversity could be achieved by modulating gene expression. Understanding the relationship between genomic and expressed MHC diversity is essential for interpreting variation in MHC diversity in an evolutionary context.

Spatiotemporal adaptive evolution of an MHC immune gene in a frog-fungus disease system

Genetic diversity of major histocompatibility complex (MHC) genes is linked to reduced pathogen susceptibility in amphibians, but few studies also examine broad spatial and temporal patterns of MHC and neutral genetic diversity. Here, we characterized range-wide MHC diversity in the Northern leopard frog, Rana pipiens, a species found throughout North America that is experiencing disease-related declines. We used previously sequenced neutral markers (mitochondrial DNA and microsatellites), sequenced an expressed MHC class IIß gene fragment, and measured infection prevalence and intensity of the global fungal pathogen Batrachochytrium dendrobatidis (Bd) across 14 populations. Four populations were sampled across two decades, enabling temporal comparisons of selection and demography. We recovered 37 unique MHC alleles, including 17 that were shared across populations. Phylogenetic and population genetic patterns between MHC and neutral markers were incongruent, and five MHC codon positions associated with peptide binding were under positive selection. MHC heterozygosity, but not neutral marker heterozygosity, was a significant factor explaining spatial patterns of Bd prevalence, whereas only environmental variables predicted Bd intensity. MHC allelic richness (AR) decreased significantly over time but microsatellite-based AR did not, highlighting a loss of functional immunogenetic diversity that may be associated with Bd selective pressures. MHC supertype 4 was significantly associated with an elevated risk of Bd infection, whereas one supertype 2 allele was associated with a nearly significant reduced risk of Bd. Taken together, these results provide evidence that positive selection contributes to MHC class IIß evolution in R. pipiens and suggest that functional MHC differences across populations may contribute to disease adaptation.

Balancing selection versus allele and supertype turnover in MHC class II genes in guppies

Selection pressure from parasites is thought to be a major force shaping the extreme polymorphism of the major histocompatibility complex (MHC) genes, but the modes and consequences of selection remain unclear. Here, we analyse MHC class II and microsatellite diversity in 16 guppy (Poecilia reticulata) populations from two islands (Trinidad and Tobago) that have been separated for at least 10 ky. Within-population MHC diversity was high, but allele sharing was limited within islands and even lower between islands, suggesting relatively fast turnover of alleles. Allelic lineages strongly supported in phylogenetic analyses tended to be island-specific, suggesting rapid lineage sorting, and an expansion of an allelic lineage private to Tobago was observed. New alleles appear to be generated locally at a detectably high frequency. We did not detect a consistent signature of local adaptation, but F_ST outlier analysis suggested that balancing selection may be the more general process behind spatial variation in MHC allele frequencies in this system, particularly within Trinidad. We found no evidence for divergent allele advantage within populations, or for decreased genetic structuring of MHC supertypes compared to MHC alleles. The dynamic and complex nature of MHC evolution we observed in guppies, coupled with some evidence for balancing selection shaping MHC allele frequencies, are consistent with Red Queen processes of host-parasite coevolution.

Sex and hatching order modulate the association between MHC-II diversity and fitness in early-life stages of a wild seabird

Genes of the major histocompatibility complex (MHC) play a pivotal role in parasite resistance, and their allelic diversity has been associated with fitness variations in several taxa. However, studies report inconsistencies in the direction of this association, with either positive, quadratic or no association being described. These discrepancies may arise because the fitness costs and benefits of MHC diversity differ among individuals depending on their exposure and immune responses to parasites. Here, we investigated in black-legged kittiwake (Rissa tridactyla) chicks whether associations between MHC class-II diversity and fitness vary with sex and hatching order. MHC-II diversity was positively associated with growth and tick clearance in female chicks, but not in male chicks. Our data also revealed a positive association between MHC-II diversity and survival in second-hatched female chicks (two eggs being the typical clutch size). These findings may result from condition-dependent parasite infections differentially impacting sexes in relation to hatching order. We thus suggest that it may be important to account for individual heterogeneities in traits that potentially exert selective pressures on MHC diversity in order to properly predict MHC-fitness associations.

No evidence for a relationship between MHC heterozygosity and life history strategy in a sample of North American undergraduates

Although allelic diversity at the major histocompatibility complex (MHC) has implications for adaptive immunity, mate choice, and social signalling, how diversity at the MHC influences the calibration of life history strategies remains largely uninvestigated. The current study investigated whether greater MHC heterozygosity was associated with markers of slower life history strategies in a sample of 789 North American undergraduates. Contrary to preregistered predictions and to previously published findings, MHC heterozygosity was not related to any of the psychological life history-relevant variables measured (including short- vs. long-term sexual strategy, temporal discounting, the Arizona life history battery, past and current health, disgust sensitivity, and Big Five personality traits). Further, no meaningful effects emerged when analysing women and men separately. Possible reasons for why the current results are inconsistent with previous work are discussed.

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.

Spatio-temporal trends in richness and persistence of bacterial communities in decline-phase water vole populations

Understanding the driving forces that control vole population dynamics requires identifying bacterial parasites hosted by the voles and describing their dynamics at the community level. To this end, we used high-throughput DNA sequencing to identify bacterial parasites in cyclic populations of montane water voles that exhibited a population outbreak and decline in 2014-2018. An unexpectedly large number of 155 Operational Taxonomic Units (OTUs) representing at least 13 genera in 11 families was detected. Individual bacterial richness was higher during declines, and vole body condition was lower. Richness as estimated by Chao2 at the local population scale did not exhibit clear seasonal or cycle phase-related patterns, but at the vole meta-population scale, exhibited seasonal and phase-related patterns. Moreover, bacterial OTUs that were detected in the low density phase were geographically widespread and detected earlier in the outbreak; some were associated with each other. Our results demonstrate the complexity of bacterial community patterns with regard to host density variations, and indicate that investigations about how parasites interact with host populations must be conducted at several temporal and spatial scales: multiple times per year over multiple years, and at both local and long-distance dispersal scales for the host(s) under consideration.

Rats as potential reservoirs for neglected zoonotic Bartonella species in Flanders, Belgium

BACKGROUND

Bartonella spp. are vector-borne pathogens transmitted to humans via blood-sucking arthropods. Rodents such as the black rat (Rattus rattus) and Norway rat (R. norvegicus) are thought to be the main reservoirs. An infection with rodent-associated Bartonella spp. may cause severe symptoms in humans such as endocarditis and neuroretinitis. The current knowledge of Bartonella prevalence in rats from western Europe is scarce.

METHODS

Rats and a few other rodent by-catches were trapped in the context of a rodenticide resistance study at different sites in Flanders, Belgium. During dissection, biometric data were collected, and spleen tissues were taken. DNA was extracted from spleen samples and tested for Bartonella spp. by conventional generic polymerase chain reaction (PCR). To determine the Bartonella species, a selected number of amplicons were sequenced and compared with GenBank entries.

RESULTS

In total, 1123 rodents were trapped. The predominate species was R. norvegicus (99.64%). Other rodents trapped included: two water voles (Arvicola amphibius, 0.18%); one colour rat (R. norvegicus forma domestica, 0.09%); and one muskrat (Ondatra zibethicus, 0.09%). PCR analysis of 1097 rodents resulted in 410 (37.37%, 95% CI: 34.50-40.31%) Bartonella spp. DNA-positive samples. Bartonella tribocorum (94.68%, 95% CI: 88.02-98.25%) was the most frequently detected Bartonella species, followed by B. grahamii (3.19%, 95% CI: 0.66-9.04%) and B. doshiae (1.06%, 95% CI: 0.03-5.79%). An uncultured Bartonella species occurred in one water vole (1.06%, 95% CI: 0.03-5.79%). There was a significantly higher Bartonella prevalence in older rats compared to juveniles and a significant difference in Bartonella prevalence concerning the localisation of trapping sites. In contrast, there was no statistically significant difference in Bartonella prevalence regarding sex, degree of urbanisation and season.

CONCLUSIONS

Based on the high prevalence found, we conclude that the Norway rat seems to be a key reservoir host for zoonotic B. tribocorum in Belgium.

Association between louse abundance and MHC II supertypes in Galápagos mockingbirds

Major histocompatibility complex class II (MHC II) is an essential molecule triggering the adaptive immune response by the presentation of pathogens to helper T cells. The association between individual MHC II variants and various parasites has become a frequent finding in studies of vertebrate populations. However, although bird ectoparasites have a significant effect on their host's fitness, and the host's immune system can regulate ectoparasitic infections, no study has yet investigated the association between MHC II polymorphism and ectoparasite infection in the populations of free-living birds. Here, we test whether an association exists between the abundance of a chewing louse (Myrsidea nesomimi) and MHC II polymorphism of its hosts, the Galápagos mockingbirds (Mimus). We have found that the presence of two MHC II supertypes (functionally differentiated clusters) was significantly associated with louse abundance. This pattern supports the theory that a co-evolutionary interaction stands behind the maintenance of MHC polymorphism. Moreover, we have found a positive correlation between louse abundance and heterophil/lymphocyte ratio (an indicator of immunological stress) that serves as an additional piece of evidence that ectoparasite burden is affected by immunological state of Galápagos mockingbirds.

Wetter climates select for higher immune gene diversity in resident, but not migratory, songbirds

Pathogen communities can vary substantially between geographical regions due to different environmental conditions. However, little is known about how host immune systems respond to environmental variation across macro-ecological and evolutionary scales. Here, we select 37 species of songbird that inhabit diverse environments, including African and Palaearctic residents and Afro-Palaearctic migrants, to address how climate and habitat have influenced the evolution of key immune genes, the major histocompatibility complex class I (MHC-I). Resident species living in wetter regions, especially in Africa, had higher MHC-I diversity than species living in drier regions, irrespective of the habitats they occupy. By contrast, no relationship was found between MHC-I diversity and precipitation in migrants. Our results suggest that the immune system of birds has evolved greater pathogen recognition in wetter tropical regions. Furthermore, evolving transcontinental migration appears to have enabled species to escape wet, pathogen-rich areas at key periods of the year, relaxing selection for diversity in immune genes and potentially reducing immune system costs.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

Does MHC heterozygosity influence microbiota form and function?

MHC molecules are essential for the adaptive immune response, and they are the most polymorphic genetic loci in vertebrates. Extreme genetic variation at these loci is paradoxical given their central importance to host health. Classic models of MHC gene evolution center on antagonistic host-pathogen interactions to promote gene diversification and allelic diversity in host populations. However, all multicellular organisms are persistently colonized by their microbiota that perform essential metabolic functions for their host and protect from infection. Here, we provide data to support the hypothesis that MHC heterozygote advantage (a main force of selection thought to drive MHC gene evolution), may operate by enhancing fitness advantages conferred by the host’s microbiome. We utilized fecal 16S rRNA gene sequences and their predicted metagenome datasets collected from multiple MHC congenic homozygote and heterozygote mouse strains to describe the influence of MHC heterozygosity on microbiome form and function. We find that in contrast to homozygosity at MHC loci, MHC heterozygosity promotes functional diversification of the microbiome, enhances microbial network connectivity, and results in enrichment for a variety of microbial functions that are positively associated with host fitness. We demonstrate that taxonomic and functional diversity of the microbiome is positively correlated in MHC heterozygote but not homozygote animals, suggesting that heterozygote microbiomes are more functionally adaptive under similar environmental conditions than homozygote microbiomes. Our data complement previous observations on the role of MHC polymorphism in sculpting microbiota composition, but also provide functional insights into how MHC heterozygosity may enhance host health by modulating microbiome form and function. We also provide evidence to support that MHC heterozygosity limits functional redundancy among commensal microbes and may enhance the metabolic versatility of their microbiome. Results from our analyses yield multiple testable predictions regarding the role of MHC heterozygosity on the microbiome that will help guide future research in the area of MHC-microbiome interactions.

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.

Major histocompatibility complex class I diversity limits the repertoire of T cell receptors

Major histocompatibility complex (MHC) genes encode proteins that initiate adaptive immune responses through the presentation of foreign antigens to T cells. The high polymorphism found at these genes, thought to be promoted and maintained by pathogen-mediated selection, contrasts with the limited number of MHC loci found in most vertebrates. Although expressing many diverse MHC genes should broaden the range of detectable pathogens, it has been hypothesized to also cause deletion of larger fractions of self-reactive T cells, leading to a detrimental reduction of the T cell receptor (TCR) repertoire. However, a key prediction of this TCR depletion hypothesis, that the TCR repertoire should be inversely related to the individual MHC diversity, has never been tested. Here, using high-throughput sequencing and advanced sequencing error correction, we provide evidence of such an association in a rodent species with high interindividual variation in the number of expressed MHC molecules, the bank vole (Myodes glareolus). Higher individual diversity of MHC class I, but not class II, was associated with smaller TCR repertoires. Our results thus provide partial support for the TCR depletion model, while also highlighting the complex, potentially MHC class-specific mechanisms by which autoreactivity may trade off against evolutionary expansion of the MHC gene family.

Plumage microbiota covaries with the major histocompatibility complex in blue petrels

To increase fitness, a wide range of vertebrates preferentially mate with partners that are dissimilar at the major histocompatibility complex (MHC) or that have high MHC diversity. Although MHC often can be assessed through olfactory cues, the mechanism by which MHC genes influence odour remains largely unclear. MHC class IIB molecules, which enable recognition and elimination of extracellular bacteria, have been suggested to influence odour indirectly by shaping odour-producing microbiota, i.e. bacterial communities. However, there is little evidence of the predicted covariation between an animal's MHC genotype and its bacterial communities in scent-producing body surfaces. Here, using high-throughput sequencing, we tested the covariation between MHC class IIB genotypes and feather microbiota in the blue petrel (Halobaena caerulea), a seabird with highly developed olfaction that has been suggested to rely on oduor cues during an MHC-based mate choice. First, we show that individuals with similar MHC class IIB profiles also have similar bacterial assemblages in their feathers. Then, we show that individuals with high MHC diversity have less diverse feather microbiota and also a reduced abundance of a bacterium of the genus Arsenophonus, a genus in which some species are symbionts of avian ectoparasites. Our results, showing that feather microbiota covary with MHC, are consistent with the hypothesis that individual MHC genotype may shape the semiochemical-producing microbiota in birds.

Specific MHC class I supertype associated with parasite infection and color morph in a wild lizard population

The major histocompatibility complex (MHC) is a large gene family that plays a central role in the immune system of all jawed vertebrates. Nonavian reptiles are underrepresented within the MHC literature and little is understood regarding the mechanisms maintaining MHC diversity in this vertebrate group. Here, we examined the relative roles of parasite-mediated selection and sexual selection in maintaining MHC class I diversity of a color polymorphic lizard. We discovered evidence for parasite-mediated selection acting via rare-allele advantage or fluctuating selection as ectoparasite load was significantly lower in the presence of a specific MHC supertype (functional clustering of alleles): supertype four. Based on comparisons between ectoparasite prevalence and load, and assessment of the impact of ectoparasite load on host fitness, we suggest that supertype four confers quantitative resistance to ticks or an intracellular tickborne parasite. We found no evidence for MHC-associated mating in terms of pair genetic distance, number of alleles, or specific supertypes. An association was uncovered between supertype four and male throat color morph. However, it is unlikely that male throat coloration acts as a signal of MHC genotype to conspecifics because we found no evidence to suggest that male throat coloration predicts male mating status. Overall, our results suggest that parasite-mediated selection plays a role in maintaining MHC diversity in this population via rare-allele advantage and/or fluctuating selection. Further work is required to determine whether sexual selection also plays a role in maintaining MHC diversity in agamid lizards.

MHC-mediated sexual selection on birdsong: Generic polymorphism, particular alleles and acoustic signals

Several hypotheses predict that the major histocompatibility complex (MHC) drives mating preference in females. Olfactory, colour or morphological traits are often found as reliable signals of the MHC profile, but the role of avian song mediating MHC-based female choice remains largely unexplored. We investigated the relationship between several MHC and acoustic features in the collared flycatcher (Ficedula albicollis), a European passerine with complex songs. We screened a fragment of the class IIB second exon of the MHC molecule, of which individuals harbour 4-15 alleles, while considerable sequence diversity is maintained at the population level. To make statistical inferences from a large number of comparisons, we adopted both null-hypothesis testing and effect size framework in combination with randomization procedures. After controlling for potential confounding factors, neither MHC allelic diversity nor the presence of particular alleles was associated remarkably with the investigated qualitative and quantitative song traits. Furthermore, genetic similarity among males based on MHC sequences was not reflected by the similarity in their song based on syllable content. Overall, these results suggest that the relationship between features of song and the allelic composition and diversity of MHC is not strong in the studied species. However, a biologically motivated analysis revealed that individuals that harbour an MHC allele that impairs survival perform songs with broader frequency range. This finding suggests that certain aspects of the song may bear reliable information concerning the MHC profile of the individuals, which can be used by females to optimize mate choice.

Lack of evidence for selection favouring MHC haplotypes that combine high functional diversity

High rates of gene duplication and the highest levels of functional allelic diversity in vertebrate genomes are the main hallmarks of the major histocompatibility complex (MHC), a multigene family with a primordial role in pathogen recognition. The usual tight linkage among MHC gene duplicates may provide an opportunity for the evolution of haplotypes that associate functionally divergent alleles and thus grant the transmission of optimal levels of diversity to coming generations. Even though such associations may be a crucial component of disease resistance, this hypothesis has been given little attention in wild populations. Here, we leveraged pedigree data from a barn owl (Tyto alba) population to characterize MHC haplotype structure across two MHC class I (MHC-I) and two MHC class IIB (MHC-IIB) duplicates, in order to test the hypothesis that haplotypes' genetic diversity is higher than expected from randomly associated alleles. After showing that MHC loci are tightly linked within classes, we found limited evidence for shifts towards MHC haplotypes combining high diversity. Neither amino acid nor functional within-haplotype diversity were significantly higher than in random sets of haplotypes, regardless of MHC class. Our results therefore provide no evidence for selection towards high-diversity MHC haplotypes in barn owls. Rather, high rates of concerted evolution may constrain the evolution of high-diversity haplotypes at MHC-I, while, in contrast, for MHC-IIB, fixed differences among loci may provide barn owls with already optimized functional diversity. This suggests that at the MHC-I and MHC-IIB respectively, different evolutionary dynamics may govern the evolution of within-haplotype diversity.

MHC class II DRB diversity predicts antigen recognition and is associated with disease severity in California sea lions naturally infected with Leptospira interrogans

We examined the associations between California sea lion MHC class II DRB (Zaca-DRB) configuration and diversity, and leptospirosis. As Zaca-DRB gene sequences are involved with antigen presentation of bacteria and other extracellular pathogens, we predicted that they would play a role in determining responses to these pathogenic spirochaetes. Specifically, we investigated whether Zaca-DRB diversity (number of genes) and configuration (presence of specific genes) explained differences in disease severity, and whether higher levels of Zaca-DRB diversity predicted the number of specific Leptospira interrogans serovars that a sea lion's serum would react against. We found that serum from diseased sea lions with more Zaca-DRB loci reacted against a wider array of serovars. Specific Zaca-DRB loci were linked to reactions with particular serovars. Interestingly, sea lions with clinical manifestation of leptospirosis that had higher numbers of Zaca-DRB loci were less likely to recover from disease than those with lower diversity, and those that harboured Zaca-DRB.C or -G were 4.5 to 5.3 times more likely to die from leptospirosis, regardless of the infective serovars. We propose that for leptospirosis, a disadvantage of having a wider range of antigen presentation might be increased disease severity due to immunopathology. Ours is the first study to examine the importance of Zaca-DRB diversity for antigen detection and disease severity following natural exposure to infective leptospires.

Genetic wealth, population health: Major histocompatibility complex variation in captive and wild ring-tailed lemurs (Lemur catta)

Across species, diversity at the major histocompatibility complex (MHC) is critical to individual disease resistance and, hence, to population health; however, MHC diversity can be reduced in small, fragmented, or isolated populations. Given the need for comparative studies of functional genetic diversity, we investigated whether MHC diversity differs between populations which are open, that is experiencing gene flow, versus populations which are closed, that is isolated from other populations. Using the endangered ring-tailed lemur (Lemur catta) as a model, we compared two populations under long-term study: a relatively "open," wild population (n = 180) derived from Bezà Mahafaly Special Reserve, Madagascar (2003-2013) and a "closed," captive population (n = 121) derived from the Duke Lemur Center (DLC, 1980-2013) and from the Indianapolis and Cincinnati Zoos (2012). For all animals, we assessed MHC-DRB diversity and, across populations, we compared the number of unique MHC-DRB alleles and their distributions. Wild individuals possessed more MHC-DRB alleles than did captive individuals, and overall, the wild population had more unique MHC-DRB alleles that were more evenly distributed than did the captive population. Despite management efforts to maintain or increase genetic diversity in the DLC population, MHC diversity remained static from 1980 to 2010. Since 2010, however, captive-breeding efforts resulted in the MHC diversity of offspring increasing to a level commensurate with that found in wild individuals. Therefore, loss of genetic diversity in lemurs, owing to small founder populations or reduced gene flow, can be mitigated by managed breeding efforts. Quantifying MHC diversity within individuals and between populations is the necessary first step to identifying potential improvements to captive management and conservation plans.

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.

Natural selection on MHC IIβ in parapatric lake and stream stickleback: Balancing, divergent, both or neither?

Major histocompatibility complex (MHC) genes encode proteins that play a central role in vertebrates' adaptive immunity to parasites. MHC loci are among the most polymorphic in vertebrates' genomes, inspiring many studies to identify evolutionary processes driving MHC polymorphism within populations and divergence between populations. Leading hypotheses include balancing selection favouring rare alleles within populations, and spatially divergent selection. These hypotheses do not always produce diagnosably distinct predictions, causing many studies of MHC to yield inconsistent or ambiguous results. We suggest a novel strategy to distinguish balancing vs. divergent selection on MHC, taking advantage of natural admixture between parapatric populations. With divergent selection, individuals with immigrant alleles will be more infected and less fit because they are susceptible to novel parasites in their new habitat. With balancing selection, individuals with locally rare immigrant alleles will be more fit (less infected). We tested these contrasting predictions using three-spine stickleback from three replicate pairs of parapatric lake and stream habitats. We found numerous positive and negative associations between particular MHC IIβ alleles and particular parasite taxa. A few allele-parasite comparisons supported balancing selection, and others supported divergent selection between habitats. But, there was no overall tendency for fish with immigrant MHC alleles to be more or less heavily infected. Instead, locally rare MHC alleles (not necessarily immigrants) were associated with heavier infections. Our results illustrate the complex relationship between MHC IIβ allelic variation and spatially varying multispecies parasite communities: different hypotheses may be concurrently true for different allele-parasite combinations.

Testing mate choice and overdominance at MH in natural families of Atlantic salmon Salmo salar

This study aimed to test mate choice and selection during early life stages on major histocompatibility (MH) genotype in natural families of Atlantic salmon Salmo salar spawners and juveniles, using nine microsatellites to reconstruct families, one microsatellite linked to an MH class I gene and one minisatellite linked to an MH class II gene. MH-based mate choice was only detected for the class I locus on the first year, with lower expected heterozygosity in the offspring of actually mated pairs than predicted under random mating. The genotype frequencies of MH-linked loci observed in the juveniles were compared with frequencies expected from Mendelian inheritance of parental alleles to detect selection during early life stages. No selection was detected on the locus linked to class I gene. For the locus linked to class II gene, observed heterozygosity was higher than expected in the first year and lower in the second year, suggesting overdominance and underdominance, respectively. Within family, juveniles' body size was linked to heterozygosity at the same locus, with longer heterozygotes in the first year and longer homozygotes in the second year. Selection therefore seems to differ from one locus to the other and from year to year.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

A quick and robust MHC typing method for free-ranging and captive primate species

Gene products of the major histocompatibility complex (MHC) of human and non-human primates play a crucial role in adaptive immunity, and most of the relevant genes not only show a high degree of variability (polymorphism) but also copy number variation (CNV) is observed. Due to this diversity, MHC proteins influence the capability of individuals to cope with various pathogens. MHC and/or MHC-linked gene products such as odorant receptor genes are thought to influence mate choice and reproductive success. Therefore, MHC typing of wild and captive primate populations is considered to be useful in conservation biology, which is, however, often hampered by the need of invasive and time-consuming methods. All intact Mhc-DRB genes in primates appear to possess a complex and highly divergent microsatellite, DRB-STR. A panel of 154 pedigreed olive baboons (Papio anubis) was examined for their DRB content by DRB-STR analysis of genomic DNA. Using the same methodology on DNA of feces samples, DRB variability of a silvery gibbon population (Hylobates moloch) (N = 24), an endangered species, could successfully be studied. In both species, length determination of the DRB-STR resulted in the definition of unique genotyping patterns that appeared to be specific for a certain chromosome. Moreover, the different STR lengths were shown to segregate with the allelic variation of the respective gene. The results obtained expand data gained previously on DRB-STR typing in macaques, great apes, and humans and strengthen the conclusion that this protocol is applicable in molecular ecology, conservation biology, and colony management, especially of endangered primate species.