Signals of major histocompatibility complex overdominance in a wild salmonid population

The major histocompatibility complex (MHC) contains the most variable genes in vertebrates, but despite extensive research, the mechanisms maintaining this polymorphism are still unresolved. One hypothesis is that MHC polymorphism is a result of balancing selection operating by overdominance, but convincing evidence for overdominant selection in natural populations has been lacking. We present strong evidence consistent with MHC-specific overdominance in a free-living population of Arctic charr (Salvelinus alpinus) in northernmost Europe. In this population, where just two MHC alleles were observed, MHC heterozygous fish had a lower parasite load, were in better condition (as estimated by a fatness indicator) and had higher survival under stress than either of the homozygotes. Conversely, there was no consistent association between these fitness measures and assumedly neutral microsatellite variability, indicating an MHC-specific effect. Our results provide convincing empirical evidence consistent with the notion that overdominance can be an important evolutionary mechanism contributing to MHC polymorphism in wild animal populations. They also support a recent simulation study indicating that the number of alleles expected to be maintained at an MHC loci can be low, even under strong heterozygote advantage.     

MHC class IIB alleles contribute to both additive and nonadditive genetic effects on survival in Chinook salmon

The genes of the major histocompatibility complex (MHC) are found in all vertebrates and are an important component of individual fitness through their role in disease and pathogen resistance. These genes are among the most polymorphic in genomes and the mechanism that maintains the diversity has been actively debated with arguments for natural selection centering on either additive or nonadditive genetic effects. Here, we use a quantitative genetics breeding design to examine the genetic effects of MHC class IIB alleles on offspring survivorship in Chinook salmon (Oncorhynchus tshawytscha). We develop a novel genetic algorithm that can be used to assign values to specific alleles or genotypes. We use this genetic algorithm to show simultaneous additive and nonadditive effects of specific MHC class IIB alleles and genotypes on offspring survivorship. The additive effect supports the rare-allele hypothesis as a potential mechanism for maintaining genetic diversity at the MHC. However, contrary to the overdominance hypothesis, the nonadditive effect led to underdominance at one heterozygous genotype, which could instead reduce variability at the MHC. Our algorithm is an advancement over traditional animal models that only partition variance in fitness to additive and nonadditive genetic effects, but do not allocate these effects to specific alleles and genotypes. Additionally, we found evidence of nonrandom segregation during meiosis in females that promotes an MHC allele that is associated with higher survivorship. Such nonrandom segregation could further reduce variability at the MHC and may explain why Chinook salmon has one of the lowest levels of MHC diversity of all vertebrates.     

Pathogen resistance and genetic variation at MHC loci

Abstract.— Balancing selection in the form of heterozygote advantage, frequency-dependent selection, or selection that varies in time and/or space, has been proposed to explain the high variation at major histocompatibility complex (MHC) genes. Here the effect of variation of the presence and absence of pathogens over time on genetic variation at multiallelic loci is examined. In the basic model, resistance to each pathogen is conferred by a given allele, and this allele is assumed to be dominant. Given that s is the selective disadvantage for homozygotes (and heterozygotes) without the resistance allele and the proportion of generations, which a pathogen is present, is e , fitnesses for homozygotes become (1 — s )^(n-1)e and the fitnesses for heterozygotes become (1 — s )^(n-2)e, where n is the number of alleles. In this situation, the conditions for a stable, multiallelic polymorphism are met even though there is no intrinsic heterozygote advantage. The distribution of allele frequencies and consequently heterozygosity are a function of the autocorrelation of the presence of the pathogen in subsequent generations. When there is a positive autocorrelation over generations, the observed heterozygosity is reduced. In addition, the effects of lower levels of selection and dominance and the influence of genetic drift were examined. These effects were compared to the observed heterozygosity for two MHC genes in several South American Indian samples. Overall, resistance conferred by specific alleles to temporally variable pathogens may contribute to the observed polymorphism at MHC genes and other similar host defense loci.     

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.     

Cryptic preference for MHC-dissimilar females in male red junglefowl,Gallus gallus

An increasing number of studies test the idea that females increase offspring fitness by biasing fertilization in favour of genetically compatible partners; however, few have investigated or controlled for corresponding preferences in males. Here, we experimentally test whether male red junglefowl, Gallus gallus, prefer genetically compatible females, measured by similarity at the major histocompatibility complex (MHC), a key gene complex in vertebrate immune function. Theory predicts that because some degree of MHC heterozygosity favours viability, individuals should prefer partners that carry MHC alleles different from their own. While male fowl showed no preference when simultaneously presented with an MHC-similar and an MHC-dissimilar female, they showed a ‘cryptic’ preference, by allocating more sperm to the most MHC-dissimilar of two sequentially presented females. These results provide the first experimental evidence that males might respond to the MHC similarity of a female through differential ejaculate expenditure. By revealing that cryptic male behaviours may bias fertilization success in favour of genetically compatible partners, this study demonstrates the need to experimentally disentangle male and female effects when studying preferences for genetically compatible partners.     

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

The fundamental role of the major histocompatibility complex (MHC) in immune recognition has led to a general consensus that the characteristically high levels of functional polymorphism at MHC genes is maintained by balancing selection operating through host–parasite coevolution. However, the actual mechanism by which selection operates is unclear. Two hypotheses have been proposed: overdominance (or heterozygote superiority) and negative frequency-dependent selection. Evidence for these hypotheses was evaluated by examining MHC–parasite relationships in an island population of water voles (Arvicola terrestris). Generalized linear mixed models were used to examine whether individual variation at an MHC class II DRB locus explained variation in the individual burdens of five different parasites. MHC genotype explained a significant amount of variation in the burden of gamasid mites, fleas (Megabothris walkeri) and nymphs of sheep ticks (Ixodes ricinus). Additionally, MHC heterozygotes were simultaneously co-infected by fewer parasite types than homozygotes. In each case where an MHC-dependent effect on parasite burden was resolved, the heterozygote genotype was associated with fewer parasites, and the heterozygote outperformed each homozygote in two of three cases, suggesting an overall superiority against parasitism for MHC heterozygote genotypes. This is the first demonstration of MHC heterozygote superiority against multiple parasites in a natural population, a mechanism that could help maintain high levels of functional MHC genetic diversity in natural populations.     

Methods for MHC genotyping in non-model vertebrates

Genes of the major histocompatibility complex (MHC) are considered a paradigm of adaptive evolution at the molecular level and as such are frequently investigated by evolutionary biologists and ecologists. Accurate genotyping is essential for understanding of the role that MHC variation plays in natural populations, but may be extremely challenging. Here, I discuss the DNA-based methods currently used for genotyping MHC in non-model vertebrates, as well as techniques likely to find widespread use in the future. I also highlight the aspects of MHC structure that are relevant for genotyping, and detail the challenges posed by the complex genomic organization and high sequence variation of MHC loci. Special emphasis is placed on designing appropriate PCR primers, accounting for artefacts and the problem of genotyping alleles from multiple, co-amplifying loci, a strategy which is frequently necessary due to the structure of the MHC. The suitability of typing techniques is compared in various research situations, strategies for efficient genotyping are discussed and areas of likely progress in future are identified. This review addresses the well established typing methods such as the Single Strand Conformation Polymorphism (SSCP), Denaturing Gradient Gel Electrophoresis (DGGE), Reference Strand Conformational Analysis (RSCA) and cloning of PCR products. In addition, it includes the intriguing possibility of direct amplicon sequencing followed by the computational inference of alleles and also next generation sequencing (NGS) technologies; the latter technique may, in the future, find widespread use in typing complex multilocus MHC systems.     

MHC diversity and the association to nematode parasitism in the yellow-necked mouse (Apodemus flavicollis)

In vertebrates, the genes of the major histocompatibility complex (MHC) are among the most debated candidates accounting for co-evolutionary processes of host-parasite interaction at the molecular level. The exceptionally high allelic polymorphism found in MHC loci is believed to be maintained by pathogen-driven selection, mediated either through heterozygous advantage or rare allele advantage (= frequency dependent selection). While investigations under natural conditions are still very rare, studies on humans or mice under laboratory conditions revealed support for both hypotheses. We investigated nematode burden and allelic diversity of a functional important MHC class II gene (DRB exon2) in free-ranging yellow-necked mice (Apodemus flavicollis). Twenty-seven distinct Apfl-DRB alleles were detected in 146 individuals with high levels of amino acid sequence divergence, especially at the antigen binding sites (ABS), indicating selection processes acting on this locus. Heterozygosity had no influence on the infection status (being infected or not), the number of different nematode infections (NNI) or the intensity of infection, measured as the individual faecal egg count (FEC). However, significant associations of specific Apfl-DRB alleles to both nematode susceptibility and resistance were found, for all nematodes as well as in separate analyses of the two most common nematodes. Apodemus flavicollis individuals carrying the alleles Apfl-DRB*5 or Apfl-DRB*15 revealed significantly higher FEC than individuals with other alleles. In contrast, the allele Apfl-DRB*23 showed a significant association to low FEC of the most common nematode. Thus, our results provide evidence for pathogen-driven selection acting through rare allele advantage under natural conditions.     

Conserved haplotype blocks within the sheep MHC and low SNP heterozygosity in the Class IIa subregion

This report describes single-nucleotide polymorphisms (SNPs) in the sheep major histocompatibility complex (MHC) class II and class III regions and provides insights into the internal structure of this important genomic complex. MHC haplotypes were deduced from sheep family trios based on genotypes from 20 novel SNPs representative of the class II region and 10 previously described SNPs spanning the class III region. All 30 SNPs exhibited Hardy-Weinberg proportions in the sheep population studied. Recombination within an extended sire haplotype was observed within the class II region for 4 of 20 sheep chromosomes, thereby supporting the presence of separated IIa and IIb subregions similar to those present in cattle. SNP heterozygosity varied across the class II and III regions. One segment of the class IIa subregion manifested very low heterozygosity for several SNPs spanning approximately 120 Kbp. This feature corresponds to a subregion within the human MHC class II region previously described as a 'SNP desert' because of its paucity of SNPs. Linkage disequilibrium (LD) was reduced at the junction separating the putative class IIb and IIa subregions and also between the class IIa and the class III subregions. The latter observation is consistent with either an unmapped physical separation at this location or more likely a boundary characterized by more frequent recombination between two conserved subregions, each manifesting high within-block LD. These results identify internal blocks of loci in the sheep MHC, within which recombination is relatively rare.     

Microsatellite variation in the equine MHC

Genes within the major histocompatibility complex (MHC) encode proteins involved in innate and adaptive immune responses. Genetic variation in this region can influence the immune response of an individual animal to challenges from a variety of pathogens; however, a complete documentation of genetic variation in the MHC is lacking for most domestic animals, including horses. To provide additional genetic markers for study of the horse MHC, or ELA (equine lymphocyte antigen), we identified 37 polymorphic microsatellite repeats in ELA and used these variations separately and together with published SNPs to investigate linkage disequilibrium (LD) and haplotype structure in a sample of Thoroughbred horses. ELA SNPs alone detected little LD, but microsatellites, either separately or combined with SNPs, revealed substantially more LD. A subset of markers in very high LD across the breadth of ELA may be predictive of structural polymorphisms or linked epistases that are important drivers of haplotype structure in Thoroughbreds.     

MHC variation in birds and reptiles

The major histocompatibility complex (MHC) has been studied in a multitude of mammals by now, but much less is known about its organisation and variation in other vertebrate species. The mammalian MHC is organised as a single gene cluster, but recent studies on birds suggest that this paradigm of MHC organisation has to be supplemented. The domestic chicken thus possesses two separate gene clusters which both contain MHC class I and class II B genes, and we have shown that the ring‐necked pheasant Phasianus colchicus also has two unlinked clusters of class II B genes. We are studying the effect of the MHC on mate choice, survival and reproductive success in natural populations of birds and reptiles. For this reason, we are developing DNA techniques to determine the animals' MHC genotype. The amplification of the hypervariable exon 3 of the class I gene from songbirds and reptiles has provided us with species specific probes that can be used in Southern blot analysis. The first results indicate very extensive variation in all studied species, that is starlings Sturnus vulgaris, great reed warblers Acrocephalus arundinaceus and water pythons Liasis fuscus. The restriction fragment length polymorphism (RFLP) analysis also suggests that the number of MHC genes is significantly larger in these species than in pheasants and domestic chickens. This revised version was published online in July 2006 with corrections to the Cover Date.     

Historical and contemporary selection of teleost MHC genes: did we leave the past behind?

The extreme polymorphism of antigen‐presenting genes of the major histocompatibility complex (MHC) has spurred intense research unparalleled for any other gene family. This applies also to teleosts where sequence information is available for 3559 MHC class I and class II allelic variants from 137 species. This review summarizes current knowledge on the origin and maintenance of diversity at classical MHC loci. Most studies identified positive selection (i.e. elevated rates of non‐synonymous over synonymous substitutions, dN/dS) as a sign of balancing selection. A meta‐analysis on nine species with sufficient numbers of class I and class II sequences revealed that recombination rate and intensity of positive selection were positively correlated, suggesting that recombination and gene conversion played a significant role in shaping the allelic repertoire. Processes that create diversity over long timescales need to be complemented by contemporary balancing selection, either through overdominance or frequency‐dependent selection, in order to explain the high allelic diversity observed today. While some evidence for overdominance exists for a few taxa (mainly salmonids) by correlating parasite infection data or survival to MHC genotypes, field or experimental data on negative frequency‐dependent selection are lacking altogether, even though some fish species are particularly suitable as model systems. Theoretical predictions suggest that negative frequency‐dependent selection is necessary to maintain the existing polymorphism. Hence, future empirical studies should focus on detecting signals that differentiate between mechanisms of contemporary selection rather than repeatedly showing historical selection events.     

Overdominance in the human genome and olfactory receptor activity

We investigate the contribution of overdominance to the maintenanceof polymorphism in the human genome during the recent evolutionof our species. Using the HapMap genotypic information, we havedetected that the Gene Ontology term "olfactory receptor activity"is a molecular function overrepresented in genes that have SNPs(Single Nucleotide Polymorphisms) showing higher than expectednumber of heterozygotes in the HapMap populations. Our resultssuggest that the diversity of a subset of human olfactory receptors(ORs) may have been maintained by balancing selection, in theform of overdominance. This observation may suggest that theloss of OR genes during the evolution of the human lineage mayhave been accompanied by an increased capability to discriminateodorants with closely similar structures.     

Disentangling the roles of natural selection and genetic drift in shaping variation at MHC immunity genes

The major histocompatibility complex (MHC) forms an integral component of the vertebrate immune response and, due to strong selection pressures, is one of the most polymorphic regions of the entire genome. Despite over 15 years of research, empirical studies offer highly contradictory explanations of the relative roles of different evolutionary forces, selection and genetic drift, acting on MHC genes during population bottlenecks. Here, we take a meta-analytical approach to quantify the results of studies into the effects of bottlenecks on MHC polymorphism. We show that the consequences of selection acting on MHC loci prior to a bottleneck event, combined with drift during the bottleneck, will result in overall loss of MHC polymorphism that is ～15% greater than loss of neutral genetic diversity. These results are counter to general expectations that selection should maintain MHC polymorphism, but do agree with the results of recent simulation models and at least two empirical studies. Notably, our results suggest that negative frequency-dependent selection could be more important than overdominance for maintaining high MHC polymorphism in pre-bottlenecked populations.     

濒危动物主要组织相容性复合体基因的研究现状及前景

主要组织相容性复合体（MHC）基因是脊椎动物体内与免疫应答调节密切相关的一个基因家族,由紧密连锁的高度多态性基因座位组成。MHC基因具有高变异性,在机体免疫系统中发挥着非常重要的作用,而且与物种的抗病性和易感性,以及种群繁殖力和生存力密切相关。对MHC基因进行研究,在种群遗传学,特别是濒危动物的保护遗传学中具有独特的优势和应用前景。     

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.     

Playing chicken with red junglefowl: identifying phenotypic markers of genetic purity in Gallus gallus

We report the results of a novel experiment, in which genetically pure male red junglefowl Gallus gallus (Richardson strain) were deliberately crossed with domestic female chickens to create contaminated lines of known purity, reaching as high as 93.75%. Phenotypic characters generally used as indicators of purity (reduced or absent female comb, male eclipse plumage, etc.) all appeared to at least some extent in domestically contaminated progeny and moreso in successively more pure generations of the experiment, suggesting that such phenotypic characters may have little, if any, utility in characterizing red junglefowl stocks as to their genetic purity.     

MHC class I variation associates with parasite resistance and longevity in tropical pythons

Using restriction fragment length polymorphism (RFLP) we identified 26 unique major histocompatibility complex (MHC) genotypes in 104 water pythons. We observed a significant independent association between reduced blood parasite load (Hepatozoon sp.) and python body length/age, presence of a specific RFLP fragment (C-fragment) and the overall number of fragments. The parasite has a negative impact on several python life-history traits such as growth, nutritional status and longevity. Thus, the C-fragment could be considered a 'good gene' (a fitness-enhancing genetic element). However, while the number of fragments affected parasite load, the association between level of parasitaemia and fragment number was not linear, and, hence, minimum parasite infection level was achieved at an intermediate number of fragments. Intermediate MHC fragment numbers were also observed among the largest/oldest pythons, suggesting that both a specific fragment and intermediate levels of MHC polymorphism enhanced python longevity. Thus, our results suggest python MHC is subject to both frequency-dependent and balancing selection.     

Resistance to malignant catarrhal fever in American bison (Bison bison) is associated with MHC class IIa polymorphisms

The Rhadinovirus ovine herpesvirus-2 (OvHV-2) is the most common causative agent of malignant catarrhal fever (MCF) in clinically susceptible ruminants including cattle and bison. American bison (Bison bison) are highly susceptible to clinical MCF. Nevertheless, approximately 20% of bison on ranches or in feedlots become infected with the virus without developing clinical disease. Defining the genetic basis for differences in susceptibility between bison could facilitate development of improved control strategies for MCF. One genetic region that influences susceptibility to infectious diseases is the major histocompatibility complex (MHC). In this study, a Bison bison (Bibi) DRB3 oligonucleotide microarray was used to type 189 bison from 10 herds where MCF outbreaks had occurred. Binary logistic regression was used to classify DRB3 alleles as resistant (R), susceptible (S) or neutral (N). Animals were reclassified using six DRB3 genotype categories: N/N, N/R, N/S, R/S, R/R and S/S. Analysis of homogeneity across herds showed that there was a herd effect. Consequently, a penalized logistic regression model was run with herd and genotype categories as the explanatory variables. The R/R genotype was associated with resistance to MCF (P = 0.0327), while the S/S genotype was associated with clinical MCF (P = 0.0069). This is the first evidence that MHC class IIa polymorphism is associated with resistance or susceptibility to OvHV-2-induced MCF.