MHC evolution in three salmonid species: a comparison between class II alpha and beta genes

The genes of the major histocompatibility complex (MHC) are amongst the most variable in vertebrates and represent some of the best candidates to study processes of adaptive evolution. However, despite the number of studies available, most of the information on the structure and function of these genes come from studies in mammals and birds in which the MHC class I and II genes are tightly linked and class II alpha exhibits low variability in many cases. Teleost fishes are among the most primitive vertebrates with MHC and represent good organisms for the study of MHC evolution because their class I and class II loci are not physically linked, allowing for independent evolution of both classes of genes. We have compared the diversity and molecular mechanisms of evolution of classical MH class II α and class II β loci in farm populations of three salmonid species: Oncorhynchus kisutch, Oncorhynchus mykiss and Salmo salar. We found single classical class II loci and high polymorphism at both class II α and β genes in the three species. Mechanisms of evolution were common for both class II genes, with recombination and point mutation involved in generating diversity and positive selection acting on the peptide-binding residues. These results suggest that the maintenance of variability at the class IIα gene could be a mechanism to increase diversity in the MHC class II in salmonids in order to compensate for the expression of one single classical locus and to respond to a wider array of parasites.     

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

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

Evaluation of two approaches to genotyping major histocompatibility complex class I in a passerine-CE-SSCP and 454 pyrosequencing

Genes of the highly dynamic major histocompatibility complex (MHC) are directly linked to individual fitness and are of high interest in evolutionary ecology and conservation genetics. Gene duplication and positive selection usually lead to high levels of polymorphism in the MHC region, making genotyping of MHC a challenging task. Here, we compare the performance of two methods for MHC class I genotyping in a passerine with highly duplicated MHC class I genes: capillary electrophoresis-single-strand conformation polymorphism (CE-SSCP) analysis and 454 GS FLX Titanium pyrosequencing. According to our findings, the number of MHC variants (called alleles for simplicity) detected by CE-SSCP is significantly lower than detected by 454. To resolve discrepancies between the two methods, we cloned and Sanger sequenced a MHC class I amplicon for an individual with high number of alleles. We found a perfect congruence between cloning/Sanger sequencing results and 454. Thus, in case of multi-locus amplification, CE-SSCP considerably underestimates individual MHC diversity. However, numbers of alleles detected by both methods are significantly correlated, although the correlation is weak (r = 0.32). Thus, in systems with highly duplicated MHC, 454 provides more reliable information on individual diversity than CE-SSCP.     

Avian blood parasite richness decreases with major histocompatibility complex class I loci number

Major histocompatibility complex (MHC) genes are among the most polymorphic in the vertebrate genome. The high allele diversity is believed to be maintained primarily by sexual and pathogen-mediated balancing selection. The number of MHC loci also varies greatly across vertebrates, most notably across birds. MHC proteins play key roles in presenting antigens on the cell surface for recognition by T cells, with class I proteins specifically targeting intracellular pathogens. Here, we explore the hypothesis that MHC class I diversity (measured as loci number) coevolves with haemosporidian parasite burden of the host. Using data on 54 bird species, we demonstrate that high-MHC class I diversity is associated with significantly lower richness of Plasmodium, Haemoproteus as well as overall haemosporidian parasite lineages, the former thus indicating more efficient protection against intracellular pathogens. Nonetheless, the latter associations were only detected when MHC diversity was assessed using cloning and not 454 pyrosequencing-based studies, nor across all genotyping methods combined. Our results indicate that high-MHC class I diversity might play a key role in providing qualitative resistance against diverse haemosporidian parasites in birds, but further clarification is needed for the origin of contrasting results when using different genotyping methods for MHC loci quantification.     

Family‐assisted inference of the genetic architecture of major histocompatibility complex variation

With their direct link to individual fitness, genes of the major histocompatibility complex (MHC) are a popular system to study the evolution of adaptive genetic diversity. However, owing to the highly dynamic evolution of the MHC region, the isolation, characterization and genotyping of MHC genes remain a major challenge. While high‐throughput sequencing technologies now provide unprecedented resolution of the high allelic diversity observed at the MHC, in many species, it remains unclear (i) how alleles are distributed among MHC loci, (ii) whether MHC loci are linked or segregate independently and (iii) how much copy number variation (CNV) can be observed for MHC genes in natural populations. Here, we show that the study of allele segregation patterns within families can provide significant insights in this context. We sequenced two MHC class I (MHC‐I) loci in 1267 European barn owls (Tyto alba), including 590 offspring from 130 families using Illumina MiSeq technology. Coupled with a high per‐individual sequencing coverage (~3000×), the study of allele segregation patterns within families provided information on three aspects of the architecture of MHC‐I variation in barn owls: (i) extensive sharing of alleles among loci, (ii) strong linkage of MHC‐I loci indicating tandem architecture and (iii) the presence of CNV in the barn owl MHC‐I. We conclude that the additional information that can be gained from high‐coverage amplicon sequencing by investigating allele segregation patterns in families not only helps improving the accuracy of MHC genotyping, but also contributes towards enhanced analyses in the context of MHC evolutionary ecology.     

Individual MHC class I and MHC class IIB diversities are associated with male and female reproductive traits in the three-spined stickleback

Genes of the major histocompatibility complex (MHC) are indispensable for pathogen defence in vertebrates. With wild-caught three-spined sticklebacks (Gasterosteus aculeatus) we conducted the first study to relate individual reproductive parameters to both MHC class I and II diversities. An optimal MHC class IIB diversity was found for male nest quality. However, male breeding colouration was most intense at a maximal MHC class I diversity. One MHC class I allele was associated with male redness. Similarly, one MHC class IIB allele was associated with continuous rather than early female reproduction, possibly extending the reproductive period. Both alleles occurred more frequently with increasing individual allele diversity. We suggest that if an allele is currently not part of the optimum, it had not been propagated by choosy females. The parasite against which this allele provides resistance is therefore unlikely to have been predominant the previous year - a step to negative frequency-dependent selection.     

Evolution of MHC class I genes in the European badger (Meles meles)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune system and provides a good model with which to understand the evolutionary processes underlying functional genes. Trans-species polymorphism and orthology are both commonly found in MHC genes; however, mammalian MHC class I genes tend to cluster by species. Concerted evolution has the potential to homogenize different loci, whereas birth-and-death evolution can lead to the loss of orthologs; both processes result in monophyletic groups within species. Studies investigating the evolution of MHC class I genes have been biased toward a few particular taxa and model species. We present the first study of MHC class I genes in a species from the superfamily Musteloidea. The European badger (Meles meles) exhibits moderate variation in MHC class I sequences when compared to other carnivores. We identified seven putatively functional sequences and nine pseudogenes from genomic (gDNA) and complementary (cDNA) DNA, signifying at least two functional class I loci. We found evidence for separate evolutionary histories of the α1 and α2/α3 domains. In the α1 domain, several sequences from different species were more closely related to each other than to sequences from the same species, resembling orthology or trans-species polymorphism. Balancing selection and probable recombination maintain genetic diversity in the α1 domain, evidenced by the detection of positive selection and a recombination event. By comparison, two recombination breakpoints indicate that the α2/α3 domains have most likely undergone concerted evolution, where recombination has homogenized the α2/α3 domains between genes, leading to species-specific clusters of sequences. Our findings highlight the importance of analyzing MHC domains separately.     

Pathogen burden,co‐infection and major histocompatibility complex variability in the European badger (Meles meles)

Pathogen‐mediated selection is thought to maintain the extreme diversity in the major histocompatibility complex (MHC) genes, operating through the heterozygote advantage, rare‐allele advantage and fluctuating selection mechanisms. Heterozygote advantage (i.e. recognizing and binding a wider range of antigens than homozygotes) is expected to be more detectable when multiple pathogens are considered simultaneously. Here, we test whether MHC diversity in a wild population of European badgers (Meles meles) is driven by pathogen‐mediated selection. We examined individual prevalence (infected or not), infection intensity and co‐infection of 13 pathogens from a range of taxa and examined their relationships with MHC class I and class II variability. This population has a variable, but relatively low, number of MHC alleles and is infected by a variety of naturally occurring pathogens, making it very suitable for the investigation of MHC–pathogen relationships. We found associations between pathogen infections and specific MHC haplotypes and alleles. Co‐infection status was not correlated with MHC heterozygosity, but there was evidence of heterozygote advantage against individual pathogen infections. This suggests that rare‐allele advantages and/or fluctuating selection, and heterozygote advantage are probably the selective forces shaping MHC diversity in this species. We show stronger evidence for MHC associations with infection intensity than for prevalence and conclude that examining both pathogen prevalence and infection intensity is important. Moreover, examination of a large number and diversity of pathogens, and both MHC class I and II genes (which have different functions), provide an improved understanding of the mechanisms driving MHC diversity.     

Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next‐generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans‐species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long‐term survival of the species.     

Diversity and evolutionary patterns of immune genes in free-ranging Namibian leopards (Panthera pardus pardus)

The genes of the major histocompatibility complex (MHC) are a key component of the mammalian immune system and have become important molecular markers for fitness-related genetic variation in wildlife populations. Currently, no information about the MHC sequence variation and constitution in African leopards exists. In this study, we isolated and characterized genetic variation at the adaptively most important region of MHC class I and MHC class II-DRB genes in 25 free-ranging African leopards from Namibia and investigated the mechanisms that generate and maintain MHC polymorphism in the species. Using single-stranded conformation polymorphism analysis and direct sequencing, we detected 6 MHC class I and 6 MHC class II-DRB sequences, which likely correspond to at least 3 MHC class I and 3 MHC class II-DRB loci. Amino acid sequence variation in both MHC classes was higher or similar in comparison to other reported felids. We found signatures of positive selection shaping the diversity of MHC class I and MHC class II-DRB loci during the evolutionary history of the species. A comparison of MHC class I and MHC class II-DRB sequences of the leopard to those of other felids revealed a trans-species mode of evolution. In addition, the evolutionary relationships of MHC class II-DRB sequences between African and Asian leopard subspecies are discussed.     

Gene duplication,allelic diversity,selection processes and adaptive value of MHC class II <Emphasis Type="Italic">DRB</Emphasis> genes of the bank vole, <Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.     

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.Subject terms: Population genetics, Evolutionary genetics     

Regulation of MHC class I and II gene transcription: differences and similarities

Major histocompatibility complex (MHC) molecules serve as peptide receptors. These peptides are derived from processed cellular or extra-cellular antigens. The MHC gene complex encodes two major classes of molecules, MHC class I and class II, whose function is to present peptides to CD8⁺ (cytotoxic) and CD4⁺ (helper) T cells, respectively. The genes encoding both classes of MHC molecules seem to originate from a common ancestral gene. One of the hallmarks of the MHC is its extensive polymorphism which displays locus and allele-specific characteristics among the various MHC class I and class II genes. Because of its central role in immunosurveillance and in various disease states, the MHC is one of the best studied genetic systems. This review addresses several aspects of MHC class I and class II gene regulation in human and in particular, the contribution to the constitutive and cytokine-induced expression of MHC class I and II genes of MHC class-specific regulatory elements and regulatory elements which apparently are shared by the promoters of MHC class I and class II genes. Received: 12 January 1998     

Two patterns of variation among MHC class I loci in Tuatara (Sphenodon punctatus)

The genes of the major histocompatibility complex (MHC) are a central component of the immune system in vertebrates and have become important markers of functional, fitness-related genetic variation. We have investigated the evolutionary processes that generate diversity at MHC class I genes in a large population of an archaic reptile species, the tuatara (Sphenodon punctatus), found on Stephens Island, Cook Strait, New Zealand. We identified at least 2 highly polymorphic (UA type) loci and one locus (UZ) exhibiting low polymorphism. The UZ locus is characterized by low nucleotide diversity and weak balancing selection and may be either a nonclassical class I gene or a pseudogene. In contrast, the UA-type alleles have high nucleotide diversity and show evidence of balancing selection at putative peptide-binding sites. Twenty-one different UA-type genotypes were identified among 26 individuals, suggesting that the Stephens Island population has high levels of MHC class I variation. UA-type allelic diversity is generated by a mixture of point mutation and gene conversion. As has been found in birds and fish, gene conversion obscures the genealogical relationships among alleles and prevents the assignment of alleles to loci. Our results suggest that the molecular mechanisms that underpin MHC evolution in nonmammals make locus-specific amplification impossible in some species.     

Modes of salmonid MHC class I and II evolution differ from the primate paradigm

Rainbow trout (Oncorhynchus mykiss) and brown trout (Salmo trutta) represent two salmonid genera separated for 15--20 million years. cDNA sequences were determined for the classical MHC class I heavy chain gene UBA and the MHC class II beta-chain gene DAB from 15 rainbow and 10 brown trout. Both genes are highly polymorphic in both species and diploid in expression. The MHC class I alleles comprise several highly divergent lineages that are represented in both species and predate genera separation. The class II alleles are less divergent, highly species specific, and probably arose after genera separation. The striking difference in salmonid MHC class I and class II evolution contrasts with the situation in primates, where lineages of class II alleles have been sustained over longer periods of time relative to class I lineages. The difference may arise because salmonid MHC class I and II genes are not linked, whereas in mammals they are closely linked. A prevalent mechanism for evolving new MHC class I alleles in salmonids is recombination in intron II that shuffles alpha 1 and alpha 2 domains into different combinations.     

Characterization of MHC class I and II genes in a subantarctic seabird,the blue petrel, <Emphasis Type="Italic">Halobaena caerulea</Emphasis> (Procellariiformes)

The great polymorphism observed in the major histocompatibility complex (MHC) genes is thought to be maintained by pathogen-mediated selection possibly combined with MHC-disassortative mating, guided by MHC-determined olfactory cues. Here, we partly characterize the MHC class I and II B of the blue petrel, Halobaena caerulea (Procellariiformes), a bird with significant olfactory abilities that lives under presumably low pathogen burdens in Subantarctica. Blue petrels are long-lived, monogamous birds which suggest the necessity of an accurate mate choice process. The species is ancestral to songbirds (Passeriformes; many MHC loci), although not to gamefowls (Galliformes; few MHC loci). Considering the phylogenetic relationships and the low subantarctic pathogen burden, we expected few rather than many MHC loci in the blue petrel. However, when we analysed partial MHC class I and class II B cDNA and gDNA sequences we found evidence for as many as at least eight MHC class I loci and at least two class II B loci. These class I and II B sequences showed classical MHC characteristics, e.g. high nucleotide diversity, especially in putative peptide-binding regions where signatures of positive selection was detected. Trans-species polymorphism was found between MHC class II B sequences of the blue petrel and those of thin-billed prion, Pachyptila belcheri, two species that diverged ∼25 MYA. The observed MHC allele richness in the blue petrel may well serve as a basis for mate choice, especially since olfactory discrimination of MHC types may be possible in this species.     

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.     

Characterization of antisera raised against stickleback (Gasterosteus aculeatus) MHC class I and class II molecules

The three-spined stickleback (Gasterosteus aculeatus) is an important model organism for investigations on the maintenance of polymorphism of the major histocompatibility complex (MHC) of vertebrates. Analysis of functional aspects of MHC diversity in stickleback would benefit from the availability of MHC specific reagents. Here we characterize antisera raised against recombinant fusion proteins of stickleback MHC class I alpha and class II alpha and beta. Western blot analysis using recombinant proteins confirmed the specificity of the antisera. In brain and muscle preparations, neither of the MHC types was detectable. High levels of each MHC receptor type were observed in gills and spleen and lower levels in head kidneys. In histological sections of gills, epithelial cells of primary and secondary lamellae stained positive with MHC class I antiserum, while single, scattered cells stained positive for MHC class II. In sections of spleen and head kidney, considerable numbers of cells positive for either MHC type were detected. Molecular weight shift in SDS-PAGE after deglycosylation of MHC class I alpha and class II beta confirmed the predicted glyco-protein character of the molecules. The majority of MHC II alpha was not glycosylated; only a small fraction of MHC II alpha was susceptible to deglycosylation. This suggests differential expression of the two stickleback MHC II alpha genes (Gaac-DAA, Gaac-DBA) only one of which (Gaac-DBA) has a site for N-linked glycosylation.     

Blood parasites shape extreme major histocompatibility complex diversity in a migratory passerine

Pathogens are one of the main forces driving the evolution and maintenance of the highly polymorphic genes of the vertebrate major histocompatibility complex (MHC). Although MHC proteins are crucial in pathogen recognition, it is still poorly understood how pathogen‐mediated selection promotes and maintains MHC diversity, and especially so in host species with highly duplicated MHC genes. Sedge warblers (Acrocephalus schoenobaenus) have highly duplicated MHC genes, and using data from high‐throughput MHC genotyping, we were able to investigate to what extent avian malaria parasites explain temporal MHC class I supertype fluctuations in a long‐term study population. We investigated infection status and infection intensities of two different strains of Haemoproteus, that is avian malaria parasites that are known to have significant fitness consequences in sedge warblers. We found that prevalence of avian malaria in carriers of specific MHC class I supertypes was a significant predictor of their frequency changes between years. This finding suggests that avian malaria infections partly drive the temporal fluctuations of the MHC class I supertypes. Furthermore, we found that individuals with a large number of different supertypes had higher resistance to avian malaria, but there was no evidence for an optimal MHC class I diversity. Thus, the two studied malaria parasite strains appear to select for a high MHC class I supertype diversity. Such selection may explain the maintenance of the extremely high number of MHC class I gene copies in sedge warblers and possibly also in other passerines where avian malaria is a common disease.