Evolutionary genetics of immunological supertypes reveals two faces of the Red Queen

In silico functional analysis of the human, chimpanzee, and gorilla MHC-A repertoires

T cells recognize peptides displayed on the surface of cells on MHC molecules. Genetic variation in MHC genes alters their peptide-binding repertoire and thus influences the potential immune response generated against pathogens. Both gorillas and chimpanzees show reduced diversity at their MHC class I A (MHC-A) locus compared to humans, which has been suggested to be the result of a pathogen-mediated selective sweep. More specifically, gorillas lack A3 lineage alleles while chimpanzees seem to have lost the A2 lineage. While previous studies showed this using phylogenetic analysis, here, we take an in silico functional approach and use the peptide-MHC binding prediction software NetMHCpan to examine the peptide-binding repertoires of common human, chimpanzee, and gorilla MHC-A molecules. We find that both gorillas and chimpanzees lack the A02 peptide binding specificity (supertype) despite gorillas being expected to have this specificity since they kept the A2 lineage. Additionally, we show that human MHC molecules with the A02 specificity bind fewer virus-derived peptides than other MHC molecules. We also do not find differential presentation of self-peptides by the A02 supertype, making the purpose of maintaining this specificity in high frequencies in the human population unclear. Taken together, we hypothesize that poor presentation of viral peptides by A02 supertype MHC molecules could have resulted in a selective sweep in chimpanzees and/or gorillas, though we could not identify a specific virus that may have caused this sweep.

The 2024 IPD-MHC database update: a comprehensive resource for major histocompatibility complex studies

The IPD-MHC Database project (http://www.ebi.ac.uk/ipd/mhc/) serves as a comprehensive and expertly curated repository for major histocompatibility complex (MHC) sequences from non-human species, providing the necessary infrastructure and tools to study the function and evolution of this highly polymorphic genomic region. In its latest version, the IPD-MHC database has expanded both in content and in the tools for data visualization and comparison. The database now hosts over 18 000 MHC alleles from 125 species, organized into eleven taxonomic groups, all manually curated and named by the Comparative MHC Nomenclature Committee. A cetacean section has recently been included, offering researchers valuable data to study the immune system of whales, dolphins, and porpoises, as well establishing the official nomenclature platform for the Cetacea Leukocyte Antigens (CeLA). In response to user demand and reflecting broader trends in bioinformatics and immunogenetics, IPD-MHC now includes the predicted tertiary structure of over 8000 alleles and allows comparison and visualisation of allele variation within and between species at single residue resolution. These latest developments maintain the critically important link between official nomenclature of curated alleles and the ability to analyse this complex polymorphism using the most up to date methods within a single repository.

What shapes a microbiome? Differences in bacterial communities associated with helminth-amphipod interactions

The fast technological advances of molecular tools have enabled us to uncover a new dimension hidden within parasites and their hosts: their microbiomes. Increasingly, parasitologists characterise host microbiome changes in the face of parasitic infections, revealing the potential of these microscopic fast-evolving entities to influence host-parasite interactions. However, most of the changes in host microbiomes seem to depend on the host and parasite species in question. Furthermore, we should understand the relative role of parasitic infections as microbiome modulators when compared with other microbiome-impacting factors (e.g., host size, age, sex). Here, we characterised the microbiome of a single intermediate host species infected by two parasites belonging to different phyla: the acanthocephalan Plagiorhynchus allisonae and a dilepidid cestode, both infecting Transorchestia serrulata amphipods collected simultaneously from the same locality. We used the v4 hypervariable region of the 16S rRNA prokaryotic gene to identify the hemolymph bacterial community of uninfected, acanthocephalan-infected, and cestode-infected amphipods, as well as the bacteria in the amphipods' immediate environment and in the parasites infecting them. Our results show that parasitic infections were more strongly associated with differences in host bacterial community richness than amphipod size, presence of amphipod eggs in female amphipods, and even parasite load. Amphipods infected by acanthocephalans had the most divergent bacterial community, with a marked decrease in alpha diversity compared with cestode-infected and uninfected hosts. In accordance with the species-specific nature of microbiome changes in parasitic infections, we found unique microbial taxa associating with hosts infected by each parasite species, as well as taxa only shared between a parasite species and their infected hosts. However, there were some bacterial taxa detected in all parasitised amphipods (regardless of the parasite species), but not in uninfected amphipods or the environment. We propose that shared bacteria associated with all hosts parasitised by distantly related helminths may be important either in helping host defences or parasites' success, and could thus interact with host-parasite evolution.

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.

Immunopeptidomics of Salmonella enterica Serovar Typhimurium-Infected Pig Macrophages Genotyped for Class II Molecules

Salmonellosis is a zoonotic infection that has a major impact on human health; consuming contaminated pork products is the main source of such infection. Vaccination responses to classic vaccines have been unsatisfactory; that is why peptide subunit-based vaccines represent an excellent alternative. Immunopeptidomics was used in this study as a novel approach for identifying antigens coupled to major histocompatibility complex class II molecules. Three homozygous individuals having three different haplotypes (Lr-0.23, Lr-0.12, and Lr-0.21) were thus selected as donors; peripheral blood macrophages were then obtained and stimulated with Salmonella typhimurium (MOI 1:40). Although similarities were observed regarding peptide length distribution, elution patterns varied between individuals; in total, 1990 unique peptides were identified as follows: 372 for Pig 1 (Lr-0.23), 438 for Pig 2 (Lr.0.12) and 1180 for Pig 3 (Lr.0.21). Thirty-one S. typhimurium unique peptides were identified; most of the identified peptides belonged to outer membrane protein A and chaperonin GroEL. Notably, 87% of the identified bacterial peptides were predicted in silico to be elution ligands. These results encourage further in vivo studies to assess the immunogenicity of the identified peptides, as well as their usefulness as possible protective vaccine candidates.

Extensive immune receptor repertoire diversity in disease-resistant rice landraces

Plants have powerful defense mechanisms and extensive immune receptor repertoires, yet crop monocultures are prone to epidemic diseases. Rice (Oryza sativa) is susceptible to many diseases, such as rice blast caused by Magnaporthe oryzae. Varietal resistance of rice to blast relies on intracellular nucleotide binding, leucine-rich repeat (NLR) receptors that recognize specific pathogen molecules and trigger immune responses. In the Yuanyang terraces in southwest China, rice landraces rarely show severe losses to disease whereas commercial inbred lines show pronounced field susceptibility. Here, we investigate within-landrace NLR sequence diversity of nine rice landraces and eleven modern varieties using complexity reduction techniques. We find that NLRs display high sequence diversity in landraces, consistent with balancing selection, and that balancing selection at NLRs is more pervasive in landraces than modern varieties. Notably, modern varieties lack many ancient NLR haplotypes that are retained in some landraces. Our study emphasizes the value of standing genetic variation that is maintained in farmer landraces as a resource to make modern crops and agroecosystems less prone to disease. The conservation of landraces is, therefore, crucial for ensuring food security in the face of dynamic biotic and abiotic threats.

MHCtools 1.5: Analysis of MHC Sequencing Data in R

The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.

Eggshell colour differences in a classic example of coevolved eggshell mimicry

Avian brood parasitism is a model system for understanding coevolutionary arms races, and the great reed warbler (Acrocephalus arundinaceus, hereafter 'warbler') and its parasite the common cuckoo (Cuculus canorus, hereafter 'cuckoo') are prime examples of this coevolutionary struggle. Here, warblers select for egg colour mimicry by rejecting poorly matched cuckoo eggs. Contrary to long-held assumptions, recent work showed that warblers tend to reject lighter and browner eggs but tended to accept darker and bluer eggs rather than basing rejection decisions solely on perceived colour differences (i.e. the degree of mimicry). This counterintuitive, colour-biased rejection behaviour would select for bluer and darker cuckoo eggs, but would only be adaptive if cuckoos were consistently lighter and browner than warbler eggs. Therefore, we tested whether warbler eggs were consistently bluer and darker than cuckoo eggs. To do so, we re-analysed eggshell reflectance spectra of warblers and the cuckoos that parasitized them in the Czech Republic. As expected, we found that warbler eggs were significantly bluer and darker than the cuckoo eggs at the population level. Thus, we demonstrate imperfect mimicry in a long-coevolved cuckoo host-race and provide insights for exploring the coevolutionary interactions among hosts and their brood parasites.

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

The impact of HLA polymorphism on herpesvirus infection and disease

Human Leukocyte Antigens (HLA) are cell surface molecules, central in coordinating innate and adaptive immune responses, that are targets of strong diversifying natural selection by pathogens. Of these pathogens, human herpesviruses have a uniquely ancient relationship with our species, where coevolution likely has reciprocating impact on HLA and viral genomic diversity. Consistent with this notion, genetic variation at multiple HLA loci is strongly associated with modulating immunity to herpesvirus infection. Here, we synthesize published genetic associations of HLA with herpesvirus infection and disease, both from case/control and genome-wide association studies. We analyze genetic associations across the eight human herpesviruses and identify HLA alleles that are associated with diverse herpesvirus-related phenotypes. We find that whereas most HLA genetic associations are virus- or disease-specific, HLA-A*01 and HLA-A*02 allotypes may be more generally associated with immune susceptibility and control, respectively, across multiple herpesviruses. Connecting genetic association data with functional corroboration, we discuss mechanisms by which diverse HLA and cognate receptor allotypes direct variable immune responses during herpesvirus infection and pathogenesis. Together, this review examines the complexity of HLA-herpesvirus interactions driven by differential T cell and Natural Killer cell immune responses.

Both selection and drift drive the spatial pattern of adaptive genetic variation in a wild mammal

The major histocompatibility complex (MHC) has been intensively studied for the relative effects of different evolutionary forces in recent decades. Pathogen-mediated balancing selection is generally thought to explain the high polymorphism observed in MHC genes, but it is still unclear to what extent MHC diversity is shaped by selection relative to neutral drift. In this study, we genotyped MHC class II DRB genes and 15 neutral microsatellite loci across 26 geographic populations of European badgers (Meles meles) covering most of their geographic range. By comparing variation of microsatellite and diversity of MHC at different levels, we demonstrate that both balancing selection and drift have shaped the evolution of MHC genes. When only MHC allelic identity was investigated, the spatial pattern of MHC variation was similar to that of microsatellites. By contrast, when functional aspects of the MHC diversity (e.g., immunological supertypes) were considered, balancing selection appears to decrease genetic structuring across populations. Our comprehensive sampling and analytical approach enable us to conclude that the likely mechanisms of selection are heterozygote advantage and/or rare-allele advantage. This study is a clear demonstration of how both balancing selection and genetic drift simultaneously affect the evolution of MHC genes in a widely distributed wild mammal.

Incipient Parallel Evolution of SARS-CoV-2 Deltacron Variant in South Brazil

With the coexistence of multiple lineages and increased international travel, recombination and gene flow are likely to become increasingly important in the adaptive evolution of SARS-CoV-2. These processes could result in genetic introgression and the incipient parallel evolution of multiple recombinant lineages. However, identifying recombinant lineages is challenging, and the true extent of recombinant evolution in SARS-CoV-2 may be underestimated. This study describes the first SARS-CoV-2 Deltacron recombinant case identified in Brazil. We demonstrate that the recombination breakpoint is at the beginning of the Spike gene. The 5' genome portion (circa 22 kb) resembles the AY.101 (Delta), and the 3' genome portion (circa 8 kb nucleotides) is most similar to the BA.1.1 (Omicron). Furthermore, evolutionary genomic analyses indicate that the new strain emerged after a single recombination event between lineages of diverse geographical locations in December 2021 in South Brazil. This Deltacron, AYBA-RS, is one of the dozens of recombinants described in 2022. The submission of only four sequences in the GISAID database suggests that this lineage had a minor epidemiological impact. However, the recent emergence of this and other Deltacron recombinant lineages (XD, XF, and XS) suggests that gene flow and recombination may play an increasingly important role in the COVID-19 pandemic. We explain the evolutionary and population genetic theory that supports this assertion, concluding that this stresses the need for continued genomic surveillance. This monitoring is vital for countries where multiple variants are present, as well as for countries that receive significant inbound international travel.

Meta-analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Pathogen-mediated selection and sexual selection are important drivers of evolution. Both processes are known to target genes of the major histocompatibility complex (MHC), a gene family encoding cell-surface proteins that display pathogen peptides to the immune system. The MHC is also a model for understanding processes such as gene duplication and trans-species allele sharing. The class II MHC protein is a heterodimer whose peptide-binding groove is encoded by an MHC-IIA gene and an MHC-IIB gene. However, our literature review found that class II MHC papers on infectious disease or sexual selection included IIA data only 18% and 9% of the time, respectively. To assess whether greater emphasis on MHC-IIA is warranted, we analysed MHC-IIA sequence data from 50 species of vertebrates (fish, amphibians, birds, mammals) to test for polymorphism and positive selection. We found that the number of MHC-IIA alleles within a species was often high, and covaried with sample size and number of MHC-IIA genes assayed. While MHC-IIA variability tended to be lower than that of MHC-IIB, the difference was only ~25%, with ~3 fewer IIA alleles than IIB. Furthermore, the unexpectedly high MHC-IIA variability showed clear signatures of positive selection in most species, and positive selection on MHC-IIA was stronger in fish than in other surveyed vertebrate groups. Our findings underscore that MHC-IIA can be an important target of selection. Future studies should therefore expand the characterization of MHC-IIA at both allelic and genomic scales, and incorporate MHC-IIA into models of fitness consequences of MHC variation.

Sequence diversity and differences at the highly duplicated MHC-I gene reflect viral susceptibility in sympatric pinniped species

Differences in disease susceptibility among species can result from rapid host-pathogen coevolution and differences in host species ecology that affect the strength and direction of natural selection. Among two sympatric pinniped species that differ in sociality and putative disease exposure, we investigate observed differences in susceptibility through an analysis of a highly variable, duplicated gene family involved in the vertebrate immune response. Using high-throughput amplicon sequencing, we characterize diversity at the two exons that encode the peptide binding region of the major histocompatibility complex class I (MHC-I) gene in harbor (N = 60) and gray (N = 90) seal populations from the Northwest Atlantic. Across species, we identified 106 full-length exon 2 and 103 exon 3 sequence variants and a minimum of 11 duplicated MHC-I loci. The sequence variants clustered in 15 supertypes defined by the physiochemical properties of the peptide binding region, including a putatively novel Northwest Atlantic MHC-I diversity sublineage. Trans-species polymorphisms, dN/dS ratios, and evidence of gene conversion among supertypes are consistent with balancing selection acting on this gene. High functional redundancy suggests particularly strong selection among gray seals at the novel Northwest Atlantic MHC-I diversity sublineage. At exon 2, harbor seals had a significantly greater number of variants per individual than gray seals, but fewer supertypes. Supertype richness and private supertypes are hypothesized to contribute to observed differences in disease resistance between species, as consistently, across the North Atlantic and many disease outbreaks, gray seals appear to be more resistant to respiratory viruses than harbor seals.

Joint estimation of selection intensity and mutation rate under balancing selection with applications to HLA

A composite likelihood method is introduced for jointly estimating the intensity of selection and the rate of mutation, both scaled by the effective population size, when there is balancing selection at a single multi-allelic locus in an isolated population at demographic equilibrium. The performance of the method is tested using simulated data. Average estimated mutation rates and selection intensities are close to the true values but there is considerable variation about the averages. Allowing for both population growth and population subdivision does not result in qualitative differences but the estimated mutation rates and selection intensities do not in general reflect the current effective population size. The method is applied to three class I (HLA-A, HLA-B and HLA-C) and two class II loci (HLA-DRB1 and HLA-DQA1) in the 1000 Genomes populations. Allowing for asymmetric balancing selection has only a slight effect on the results from the symmetric model. Mutations that restore symmetry of the selection model are preferentially retained because of the tendency of natural selection to maximize average fitness. However, slight differences in selective effects result in much longer persistence time of some alleles. Trans-species polymorphism (TSP), which is characteristic of major-histocompatibility loci in vertebrates, is more likely when there are small differences in allelic fitness than when complete symmetry is assumed. Therefore, variation in allelic fitness expands the range of parameter values consistent with observations of TSP.

Divergent and non-parallel evolution of MHC IIB in the Neotropical Midas cichlid species complex

Background

Ecological diversification is the result of divergent natural selection by contrasting habitat characteristics that favours the evolution of distinct phenotypes. This process can happen in sympatry and in allopatry. Habitat-specific parasite communities have the potential to drive diversification among host populations by imposing selective pressures on their host's immune system. In particular, the hyperdiverse genes of the major histocompatibility complex (MHC) are implicated in parasite-mediated host divergence. Here, we studied the extent of divergence at MHC, and discuss how it may have contributed to the Nicaraguan Midas cichlid species complex diversification, one of the most convincing examples of rapid sympatric parallel speciation.

Results

We genotyped the MHC IIB for individuals from six sympatric Midas cichlid assemblages, each containing species that have adapted to exploit similar habitats. We recovered large allelic and functional diversity within the species complex. While most alleles were rare, functional groups of alleles (supertypes) were common, suggesting that they are key to survival and that they were maintained during colonization and subsequent radiations. We identified lake-specific and habitat-specific signatures for both allelic and functional diversity, but no clear pattern of parallel divergence among ecomorphologically similar phenotypes.

Conclusions

Colonization and demographic effects of the fish could have contributed to MHC evolution in the Midas cichlid in conjunction with habitat-specific selective pressures, such as parasites associated to alternative preys or environmental features. Additional ecological data will help evaluating the role of host–parasite interactions in the Midas cichlid radiations and aid in elucidating the potential role of non-parallel features differentiating crater lake species assemblages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01997-9.

Variation in volatile organic compounds in Atlantic salmon mucus is associated with resistance to salmon lice infection

Salmon lice are ectoparasites that threaten wild and farmed salmonids. Artificial selection of salmon for resistance to the infectious copepodid lice stage currently relies on in vivo challenge trials on thousands of salmon a year. We challenged 5750 salmon with salmon lice (Lepeophtheirus salmonis) from two distinct farmed strains of salmon in two separate trials. We found that volatile organic compounds (VOC), 1-penten-3-ol, 1-octen-3-ol and 6-methyl-5-hepten-2-one in the mucus of the salmon host after salmon lice infection, were significantly associated with lice infection numbers across a range of water temperatures (5 °C, 10 °C, 17 °C). Some VOCs (benzene, 1-octen-3-ol and 3,5,5-trimethyl-2-hexene) were significantly different between lines divergently selected for salmon lice resistance. In a combined population assessment, selected VOCs varied between families in the range of 47- 59% indicating a genetic component and were positively correlated to the salmon hosts estimated breeding values 0.59–0.74. Mucosal VOC phenotypes could supplement current breeding practices and have the potential to be a more direct and ethical proxy for salmon lice resistance provided they can be measured prior to lice infestation.

Global population genomics of two subspecies of Cryptosporidium hominis during 500 years of evolution

Cryptosporidiosis is a major global health problem and a primary cause of diarrhea, particularly in young children in low- and middle-income countries (LMICs). The zoonotic Cryptosporidium parvum and anthroponotic Cryptosporidium hominis cause most human infections. Here, we present a comprehensive whole-genome study of C. hominis, comprising 114 isolates from 16 countries within five continents. We detect two lineages with distinct biology and demography, which diverged circa 500 years ago. We consider these lineages two subspecies and propose the names C. hominis hominis and C. hominis aquapotentis (gp60 subtype IbA10G2). In our study, C. h. hominis is almost exclusively represented by isolates from LMICs in Africa and Asia and appears to have undergone recent population contraction. In contrast, C. h. aquapotentis was found in high-income countries, mainly in Europe, North America, and Oceania, and appears to be expanding. Notably, C. h. aquapotentis is associated with high rates of direct human-to-human transmission, which may explain its success in countries with well-developed environmental sanitation infrastructure. Intriguingly, we detected genomic regions of introgression following secondary contact between the subspecies. This resulted in high diversity and divergence in genomic islands of putative virulence genes, including muc5 (CHUDEA2_430) and a hypothetical protein (CHUDEA6_5270). This diversity is maintained by balancing selection, suggesting a co-evolutionary arms race with the host. Finally, we find that recent gene flow from C. h. aquapotentis to C. h. hominis, likely associated with increased human migration, maybe driving the evolution of more virulent C. hominis variants.

MHC class I-α population differentiation in a commercial fish, the European sea bass (Dicentrarchus labrax)

Identifying population structuring in highly fecund marine species with high dispersal rates is challenging, but critical for conservation and stock delimitation for fisheries management. European sea bass (Dicentrarchus labrax) is a commercial species of fisheries and aquaculture relevance whose stocks are declining in the North Atlantic, despite management measures to protect them and identifying their fine population structure is needed for managing their exploitation. As for other marine fishes, neutral genetic markers indicate that eastern Atlantic sea bass form a panmictic population and is currently managed as arbitrarily divided stocks. The genes of the major histocompatibility complex (MHC) are key components of the adaptive immune system and ideal candidates to assess fine structuring arising from local selective pressures. We used Illumina sequencing to characterise allelic composition and signatures of selection at the MHC class I-α region of six D. labrax populations across the Atlantic range. We found high allelic diversity driven by positive selection, corresponding to moderate supertype diversity, with 131 alleles clustering into four to eight supertypes, depending on the Bayesian information criterion threshold applied, and a mean number of 13 alleles per individual. Alleles could not be assigned to particular loci, but private alleles allowed us to detect regional genetic structuring not found previously using neutral markers. Our results suggest that MHC markers can be used to detect cryptic population structuring in marine species where neutral markers fail to identify differentiation. This is particularly critical for fisheries management, and of potential use for selective breeding or identifying escapes from sea farms.

Coevolutionary theory of hosts and parasites

Host and parasite evolution are closely intertwined, with selection for adaptations and counter-adaptations forming a coevolutionary feedback loop. Coevolutionary dynamics are often difficult to intuit due to these feedbacks and are hard to demonstrate empirically in most systems. Theoretical models have therefore played a crucial role in shaping our understanding of host-parasite coevolution. Theoretical models vary widely in their assumptions, approaches and aims, and such variety makes it difficult, especially for non-theoreticians and those new to the field, to: (1) understand how model approaches relate to one another; (2) identify key modelling assumptions; (3) determine how model assumptions relate to biological systems; and (4) reconcile the results of different models with contrasting assumptions. In this review, we identify important model features, highlight key results and predictions and describe how these pertain to model assumptions. We carry out a literature survey of theoretical studies published since the 1950s (n = 219 papers) to support our analysis. We identify two particularly important features of models that tend to have a significant qualitative impact on the outcome of host-parasite coevolution: population dynamics and the genetic basis of infection. We also highlight the importance of other modelling features, such as stochasticity and whether time proceeds continuously or in discrete steps, that have received less attention but can drastically alter coevolutionary dynamics. We finish by summarizing recent developments in the field, specifically the trend towards greater model complexity, and discuss likely future directions for research.

Trans-specific polymorphism and the convergent evolution of supertypes in major histocompatibility complex class II genes in darters (Etheostoma)

Major Histocompatibility Complex (MHC) genes are one of the most polymorphic gene groups known in vertebrates. MHC genes also exhibit allelic variants that are shared among taxa, referred to as trans-specific polymorphism (TSP). The role that selection plays in maintaining such high diversity within species, as well as TSP, is an ongoing discussion in biology. In this study, we used deep-sequencing techniques to characterize MHC class IIb gene diversity in three sympatric species of darters. We found at least 5 copies of the MHC gene in darters, with 126 genetic variants encoding 122 unique amino acid sequences. We identified four supertypes based on the binding properties of proteins encoded by the sequences. Although each species had a unique pool of variants, many variants were shared between species pairs and across all three species. Phylogenetic analysis showed that the variants did not group together monophyletically based on species identity or on supertype. An expanded phylogenetic analysis showed that some darter alleles grouped together with alleles from other percid fishes. Our findings show that TSP occurs in darters, which suggests that balancing selection is acting at the genotype level. Supertypes, however, are most likely evolving convergently, as evidenced by the fact that alleles do not form monophyletic groups based on supertype. Our research demonstrates that selection may be acting differently on MHC genes at the genotype and supertype levels, selecting for the maintenance of high genotypic diversity while driving the convergent evolution of similar MHC phenotypes across different species.

Ecological basis and genetic architecture of crypsis polymorphism in the desert clicker grasshopper (Ligurotettix coquilletti)

Color polymorphic species can offer exceptional insight into the ecology and genetics of adaptation. Although the genetic architecture of animal coloration is diverse, many color polymorphisms are associated with large structural variants and maintained by biotic interactions. Grasshoppers are notably polymorphic in both color and karyotype, which makes them excellent models for understanding the ecological drivers and genetic underpinnings of color variation. Banded and uniform morphs of the desert clicker grasshopper (Ligurotettix coquilletti) are found across the western deserts of North America. To address the hypothesis that predation maintains local color polymorphism and shapes regional crypsis variation, we surveyed morph frequencies and tested for covariation with two predation environments. Morphs coexisted at intermediate frequencies at most sites, consistent with local balancing selection. Morph frequencies covaried with the appearance of desert substrate-an environment used only by females-suggesting that ground-foraging predators are major agents of selection on crypsis. We next addressed the hypothesized link between morph variation and genome structure. To do so, we designed an approach for detecting inversions and indels using only RADseq data. The banded morph was perfectly correlated with a large putative indel. Remarkably, indel dominance differed among populations, a rare example of dominance evolution in nature.

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.

Differential patterns of diversity at neutral and adaptive loci in endangered Rhodeus pseudosericeus populations

Given the fact that threatened species are often composed of isolated small populations, spatial continuity or demography of the populations may be major factors that have shaped the species' genetic diversity. Thus, neutral loci have been the most commonly-used markers in conservation genetics. However, the populations under the influence of different environmental factors may have evolved in response to different selective pressures, which cannot be fully reflected in neutral genetic variation. Rhodeus pseudosericeus, a bitterling species (Acheilognathidae; Cypriniformes) endemic to the Korean Peninsula, are only found in some limited areas of three rivers, Daecheon, Han and Muhan, that flow into the west coast. Here, we genotyped 24 microsatellite loci and two loci (DAB1 and DAB3) of MHC class II peptide-binding β1 domain for 222 individuals collected from seven populations. Our microsatellite analysis revealed distinctive differentiation between the populations of Daecheon and Muhan Rivers and the Han River populations, and populations were structured into two subgroups within the Han River. Apparent positive selection signatures were found in the peptide-binding residues (PBRs) of the MHC loci. The allelic distribution of MHC showed a degree of differentiation between the populations of Daecheon and Muhan Rivers and the Han River populations, partially similar to the results obtained for microsatellites, however showed rather complex patterns among populations in the Han River. Considering the apparent differences in the distribution of supertypes obtained based on the physicochemical differences induced by the polymorphisms of these PBRs, the differentiation in DAB1 between the two regional groups may result in the differences in immune function. No differentiation between these two regions was observed in the supertyping of DAB3, probably indicating that only DAB1 was associated with the response to locally specialized antigenic peptides.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen-mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC-DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human-mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

Depletion of MHC supertype during domestication can compromise immunocompetence

The major histocompatibility complex (MHC) plays an important role in infectious disease resistance. The presence of certain MHC alleles and functionally similar groups of MHC alleles (i.e., supertypes) has been associated with resistance to particular parasite species. Farmed and domesticated fish stocks are often depleted in their MHC alleles and supertype diversity, possibly as a consequence of artificial selection for desirable traits, inbreeding (loss of heterozygosity), genetic drift (loss of allelic diversity) and/or reduced parasite biodiversity. Here we quantify the effects of depletion of MHC class II genotype and supertype variation on resistance to the parasite Gyrodactylus turnbulli in guppies (Poecilia reticulata). Compared to the descendants of wild‐caught guppies, ornamental fish had a significantly reduced MHC variation (i.e., the numbers of MHC alleles and supertypes per individual, and per population). In addition, ornamental fish were significantly more susceptible to G. turnbulli infections, accumulating peak intensity 10 times higher than that of their wildtype counterparts. Four out of 13 supertypes were associated with a significantly reduced parasite load, and the presence of some supertypes had a dramatic effect on the intensity of infection. Remarkably, the ornamental and wildtype fish differed in the supertypes that were associated with parasite resistance. Analysis with a genetic algorithm showed that resistance‐conferring supertypes of the wildtype and ornamental fish shared two unique amino acids in the peptide‐binding region of the MHC that were not found in any other alleles. These data show that the supertype demarcation captures some, but not all, of the variation in the immune function of the alleles. This study highlights the importance of managing functional MHC diversity in livestock, and suggests there might be some immunological redundancy among MHC supertypes.

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.

Analysis of structural variants in four African cichlids highlights an association with developmental and immune related genes

BACKGROUND

East African lake cichlids are one of the most impressive examples of an adaptive radiation. Independently in Lake Victoria, Tanganyika, and Malawi, several hundreds of species arose within the last 10 million to 100,000 years. Whereas most analyses in cichlids focused on nucleotide substitutions across species to investigate the genetic bases of this explosive radiation, to date, no study has investigated the contribution of structural variants (SVs) in the evolution of adaptive traits across the three Great Lakes of East Africa.

RESULTS

Here, we annotate and characterize the repertoires and evolutionary potential of different SV classes (deletion, duplication, inversion, insertions and translocations) in four cichlid species: Haplochromis burtoni, Metriaclima zebra, Neolamprologus brichardi and Pundamilia nyererei. We investigate the patterns of gain and loss evolution for each SV type, enabling the identification of lineage specific events. Both deletions and inversions show a significant overlap with SINE elements, while inversions additionally show a limited, but significant association with DNA transposons. Inverted regions are enriched for genes regulating behaviour, or involved in skeletal and visual system development. We also find that duplicated regions show enrichment for genes associated with "antigen processing and presentation" and other immune related categories. Our pipeline and results were further tested by PCR validation of selected deletions and inversions, which confirmed respectively 7 out of 10 and 6 out of 9 events.

CONCLUSIONS

Altogether, we provide the first comprehensive overview of rearrangement evolution in East African cichlids, and some important insights into their likely contribution to adaptation.

Divergent Allele Advantage at Human MHC Genes: Signatures of Past and Ongoing Selection

The highly polymorphic genes of the major histocompatibility complex (MHC) play a key role in adaptive immunity. Divergent allele advantage, a mechanism of balancing selection, is proposed to contribute to their exceptional polymorphism. It assumes that MHC genotypes with more divergent alleles allow for broader antigen-presentation to immune effector cells, by that increasing immunocompetence. However, the direct correlation between pairwise sequence divergence and the corresponding repertoire of bound peptides has not been studied systematically across different MHC genes. Here, we investigated this relationship for five key classical human MHC genes (human leukocyte antigen; HLA-A, -B, -C, -DRB1, and -DQB1), using allele-specific computational binding prediction to 118,097 peptides derived from a broad range of human pathogens. For all five human MHC genes, the genetic distance between two alleles of a heterozygous genotype was positively correlated with the total number of peptides bound by these two alleles. In accordance with the major antigen-presentation pathway of MHC class I molecules, HLA-B and HLA-C alleles showed particularly strong correlations for peptides derived from intracellular pathogens. Intriguingly, this bias coincides with distinct protein compositions between intra- and extracellular pathogens, possibly suggesting adaptation of MHC I molecules to present specifically intracellular peptides. Eventually, we observed significant positive correlations between an allele’s average divergence and its population frequency. Overall, our results support the divergent allele advantage as a meaningful quantitative mechanism through which pathogen-mediated selection leads to the evolution of MHC diversity.

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.

Parasite diversity and ecology in a model species, the guppy (Poecilia reticulata) in Trinidad

The guppy (Poecilia reticulata) is a model species in ecology and evolution. Many studies have examined effects of predators on guppy behaviour, reproduction, survival strategies, feeding and other life-history traits, but few have studied variation in their parasite diversity. We surveyed parasites of 18 Trinidadian populations of guppy, to provide insight on the geographical mosaic of parasite variability, which may act as a source of natural selection acting on guppies. We found 21 parasite species, including five new records for Trinidad. Spatial variation in parasite diversity was significantly higher than that of piscine predators, and significant variation in parasite richness among individuals and populations was correlated with: (i) host size, (ii) snail species richness, and (iii) the distance between populations. Differences in parasite species richness are likely to play an important, yet underestimated role in the biology of this model species of vertebrate ecology and evolution.

Lack of statistical power as a major limitation in understanding MHC-mediated immunocompetence in wild vertebrate populations

Disentangling the sources of variation in developing an effective immune response against pathogens is of major interest to immunoecology and evolutionary biology. To date, the link between immunocompetence and genetic variation at the major histocompatibility complex (MHC) has received little attention in wild animals, despite the key role of MHC genes in activating the adaptive immune system. Although several studies point to a link between MHC and immunocompetence, negative findings have also been reported. Such disparate findings suggest that limited statistical power might be affecting studies on this topic, owing to insufficient sample sizes and/or a generally small effect of MHC on the immunocompetence of wild vertebrates. To clarify this issue, we investigated the link between MHC variation and seven immunocompetence proxies in a large sample of barn owls and estimated the effect sizes and statistical power of this and published studies on this topic. We found that MHC poorly explained variation in immunocompetence of barn owls, with small-to-moderate associations between MHC and immunocompetence in owls (effect size: .1 ≥ r ≤ .3) similar to other vertebrates studied to date. Such small-to-moderate effects were largely associated with insufficient power, which was only sufficient (>0.8) to detect moderate-to-large effect sizes (r ≥ .3). Thus, studies linking MHC variation with immunocompetence in wild populations are underpowered to detect MHC effects, which are likely to be of generally small magnitude. Larger sample sizes (>200) will be required to achieve sufficient power in future studies aiming to robustly test for a link between MHC variation and immunocompetence.

Red Queen revisited: Immune gene diversity and parasite load in the asexual Poecilia formosa versus its sexual host species P. mexicana

In accordance with the Red Queen hypothesis, the lower genotypic diversity in clonally reproducing species should make them easier targets for pathogen infection, especially when closely related sexually reproducing species occur in close proximity. We analyzed two populations of clonal P. formosa and their sexual parental species P. mexicana by correlating individual parasite infection with overall and immune genotype. Our study revealed lower levels of overall genotypic diversity and marginally fewer MHC class I alleles in P. formosa individuals compared to sexually reproducing P. mexicana. Parasite load, however, differed only between field sites but not between species. We hypothesize that this might be due to slightly higher genotypic diversity in P. formosa at the innate immune system (toll like receptor 8) which is likely due to the species’ hybrid origin. In consequence, it appears that clonal individuals do not necessarily suffer a disadvantage compared to sexual individuals when fighting parasite infection.

Arms races with mitochondrial genome soft sweeps in a gynodioecious plant, Plantago lanceolata

Biological situations involving conflict can create arms race situations with repeated fixations of different functional variants, producing selective sweeps and lowering neutral diversity in genome regions linked to the functional locus. However, they can sometimes lead to balancing selection, potentially creating long coalescent times for sites with functionally different variants, and, if recombination occurs rarely, for extended haplotypes carrying such variants. We tested between these possibilities in a gynodioecious plant, Plantago lanceolata, in which cytoplasmic male-sterility factors conflict with nuclear restorers of male fertility. We find low mitochondrial diversity, which does not support very long-term coexistence of highly diverged mitochondrial haplotypes. Interestingly, however, we found a derived haplotype that is associated with male fertility in a restricted geographic region, and that has fixed differences from the ancestral sequence in several genes, suggesting that it did not arise very recently. Taken together, the results suggest arms race events that involved "soft" selective sweeps involving a moderately old-established haplotype, consistent with the frequency fluctuations predicted by theoretical models of gynodioecy.

Reduced bonobo MHC class I diversity predicts a reduced viral peptide binding ability compared to chimpanzees

BACKGROUND

The highly polymorphic genes of the major histocompatibility complex (MHC) class I are involved in defense against viruses and other intracellular pathogens. Although several studies found reduced MHC class I diversity in bonobos in comparison to the closely related chimpanzee, it is unclear if this lower diversity also influences the functional ability of MHC class I molecules in bonobos. Here, we use a bioinformatic approach to analyze the viral peptide binding ability of all published bonobo MHC class I molecules (n = 58) in comparison to all published chimpanzee MHC class I molecules (n = 161) for the class I loci A, B, C and A-like.

RESULTS

We examined the peptide binding ability of all 219 different MHC class I molecules to 5,788,712 peptides derived from 1432 different primate viruses and analyzed the percentage of bound peptides and the overlap of the peptide binding repertoires of the two species. We conducted multiple levels of analysis on the "species"-, "population"- and "individual"-level to account for the characterization of MHC variation in a larger number of chimpanzees and their broader geographic distribution. We found a lower percentage of bound peptides in bonobos at the B locus in the "population"-level comparison and at the B and C loci in the "individual"-level comparison. Furthermore, we found evidence of a limited peptide binding repertoire in bonobos by tree-based visualization of functional clustering of MHC molecules, as well as an analysis of peptides bound by both species.

CONCLUSION

Our results suggest a reduced MHC class I viral peptide binding ability at the B and C loci in bonobos compared to chimpanzees. The effects of this finding on the immune defense against viruses in wild living bonobos are unclear. However, special caution is needed to prevent introduction and spread of new viruses to bonobos, as their defensive ability to cope with new viruses could be limited compared to chimpanzees.

Antibiotic-driven escape of host in a parasite-induced Red Queen dynamics

Winnerless coevolution of hosts and parasites could exhibit Red Queen dynamics, which is characterized by parasite-driven cyclic switching of expressed host phenotypes. We hypothesize that the application of antibiotics to suppress the reproduction of parasites can provide an opportunity for the hosts to escape such winnerless coevolution. Here, we formulate a minimal mathematical model of host-parasite interaction involving multiple host phenotypes that are targeted by adapting parasites. Our model predicts the levels of antibiotic effectiveness that can steer the parasite-driven cyclic switching of host phenotypes (oscillations) to a stable equilibrium of host survival. Our simulations show that uninterrupted application of antibiotic with high-level effectiveness (greater than 85%) is needed to escape the Red Queen dynamics. Interrupted and low level of antibiotic effectiveness are indeed useless to stop host-parasite coevolution. This study can be a guide in designing good practices and protocols to minimize the risk of further progression of parasitic infections.

The Effect of Balancing Selection on Population Differentiation: A Study with HLA Genes

Balancing selection is defined as a class of selective regimes that maintain polymorphism above what is expected under neutrality. Theory predicts that balancing selection reduces population differentiation, as measured by FST. However, balancing selection regimes in which different sets of alleles are maintained in different populations could increase population differentiation. To tackle the connection between balancing selection and population differentiation, we investigated population differentiation at the HLA genes, which constitute the most striking example of balancing selection in humans. We found that population differentiation of single nucleotide polymorphisms (SNPs) at the HLA genes is on average lower than that of SNPs in other genomic regions. We show that these results require using a computation that accounts for the dependence of FST on allele frequencies. However, in pairs of closely related populations, where genome-wide differentiation is low, differentiation at HLA is higher than in other genomic regions. Such increased population differentiation at HLA genes for recently diverged population pairs was reproduced in simulations of overdominant selection, as long as the fitness of the homozygotes differs between the diverging populations. The results give insight into a possible "divergent overdominance" mechanism for the nature of balancing selection on HLA genes across human populations.

Immunogenetic novelty confers a selective advantage in host-pathogen coevolution

The major histocompatibility complex (MHC) is crucial to the adaptive immune response of vertebrates and is among the most polymorphic gene families known. Its high diversity is usually attributed to selection imposed by fast-evolving pathogens. Pathogens are thought to evolve to escape recognition by common immune alleles, and, hence, novel MHC alleles, introduced through mutation, recombination, or gene flow, are predicted to give hosts superior resistance. Although this theoretical prediction underpins host-pathogen "Red Queen" coevolution, it has not been demonstrated in the context of natural MHC diversity. Here, we experimentally tested whether novel MHC variants (both alleles and functional "supertypes") increased resistance of guppies (Poecilia reticulata) to a common ectoparasite (Gyrodactylus turnbulli). We used exposure-controlled infection trials with wild-sourced parasites, and Gyrodactylus-naïve host fish that were F2 descendants of crossed wild populations. Hosts carrying MHC variants (alleles or supertypes) that were new to a given parasite population experienced a 35-37% reduction in infection intensity, but the number of MHC variants carried by an individual, analogous to heterozygosity in single-locus systems, was not a significant predictor. Our results provide direct evidence of novel MHC variant advantage, confirming a fundamental mechanism underpinning the exceptional polymorphism of this gene family and highlighting the role of immunogenetic novelty in host-pathogen coevolution.