Trans-species polymorphism and evidence of selection on class II MHC loci in tuco-tucos (Rodentia: Ctenomyidae)

Balancing selection acting over the evolutionary history of a lineage can result in the retention of alleles among species for longer than expected under neutral evolution. The associated pattern of trans-species polymorphism, in which similar or even identical alleles are shared among species, is often used to infer that balancing selection has occurred. The genes of the major histocompatibility complex (MHC) are thought to be subject to balancing selection that maintains alleles associated with response to specific pathogens. To explore the role of balancing selection in shaping MHC diversity in ctenomyid rodents, we examined allelic variability at the class II DRB and DQA loci in 18 species in the genus Ctenomys. Previous studies of four of these species had revealed significant within-population evidence of positive selection on MHC loci. The current study expands upon these analyses to (1) evaluate among-species evidence of positive selection and (2) explore the potential for balancing selection on MHC genes. Interspecific nucleotide sequence variation revealed significant evidence of positive selection on the DRB and DQA loci. At the same time, comparisons of phylogenetic trees for these MHC loci with a putative species tree based on mitochondrial sequence data revealed multiple examples of trans-specific polymorphism, including sharing of identical DRB and DQA alleles among distantly related species of Ctenomys. These findings suggest that MHC genes in these animals have historically been subject to balancing selection and yield new insights into the complex suite of forces shaping MHC diversity in free-living vertebrates.     

Microrecombinations generate sequence diversity in the murine major histocompatibility complex: analysis of the Kbm3, Kbm4, Kbm10, and Kbm11 mutants.

The mechanism that generates spontaneous mutants of the Kb histocompatibility gene was analyzed. Nucleotide sequence analysis of four mutant genes (Kbm3, Kbm4, Kbm10, and Kbm11) revealed that each mutant K gene contains clustered, multiple nucleotide substitutions. Hybridization analyses of parental B6 genomic DNA and cloned class I genes with mutant-specific oligonucleotide probes, followed by sequence analyses, have identified major histocompatibility complex class I genes in the K, D, and Tla regions (K1, Db, and T5, respectively) that contain the exact sequences as substituted into mutant Kb genes. These data provide evidence for the hypothesis that the mutant Kb genes are generated by a microrecombination (gene conversion) mechanism that results in the transfer of small DNA segments from class I genes of all four regions of the major histocompatibility complex (K, D, Qa, and Tla) to Kb. Many of the nucleotides substituted into the mutant Kb genes were identical to those found in other naturally occurring K alleles such as Kd. Thus, we propose that the accumulation of microrecombination products within the K genes of a mouse population is responsible for the high sequence diversity among H-2 alleles.     

Distinct evolutionary trajectories of MHC class I and class II genes in Old World finches and buntings

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.Subject terms: Molecular evolution, Immunogenetics     

Evidence of gene orthology and trans‐species polymorphism,but not of parallel evolution,despite high levels of concerted evolution in the major histocompatibility complex of flamingo species

The major histocompatibility complex (MHC) is a cornerstone in the study of adaptive genetic diversity. Intriguingly, highly polymorphic MHC sequences are often not more similar within species than between closely related species. Divergent selection of gene duplicates, balancing selection maintaining trans‐species polymorphism (TSP) that predate speciation and parallel evolution of species sharing similar selection pressures can all lead to higher sequence similarity between species. In contrast, high rates of concerted evolution increase sequence similarity of duplicated loci within species. Assessing these evolutionary models remains difficult as relatedness and ecological similarities are often confounded. As sympatric species of flamingos are more distantly related than allopatric species, flamingos represent an ideal model to disentangle these evolutionary models. We characterized MHC Class I exon 3, Class IIB exon 2 and exon 3 of the six extant flamingo species. We found up to six MHC Class I loci and two MHC Class IIB loci. As all six species shared the same number of MHC Class IIB loci, duplication appears to predate flamingo speciation. However, the high rate of concerted evolution has prevented the divergence of duplicated loci. We found high sequence similarity between all species regardless of codon position. The latter is consistent with balancing selection maintaining TSP, as under this mechanism amino acid sites under pathogen‐mediated selection should be characterized by fewer synonymous codons (due to their common ancestry) than under parallel evolution. Overall, balancing selection maintaining TSP appears to result in high MHC similarity between species regardless of species relatedness and geographical distribution.     

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.     

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.     

Greater prairie chickens have a compact MHC-B with a single class IA locus

The major histocompatibility complex (MHC) plays a central role in innate and adaptive immunity, but relatively little is known about the evolution of the number and arrangement of MHC genes in birds. Insights into the evolution of the MHC in birds can be gained by comparing the genetic architecture of the MHC between closely related species. We used a fosmid DNA library to sequence a 60.9-kb region of the MHC of the greater prairie chicken (Tympanuchus cupido), one of five species of Galliformes with a physically mapped MHC. Greater prairie chickens have the smallest core MHC yet observed in any bird species, and major changes are observed in the number and arrangement of MHC loci. In particular, the greater prairie chicken differs from other Galliformes in the deletion of an important class I antigen binding gene. Analysis of the remaining class IA gene in a population of greater prairie chickens in Wisconsin, USA revealed little evidence for selection at the region responsible for antigen binding.     

Differential modes of MHC class IIB gene evolution in cichlid fishes

Cichlid fishes are emblematic models for the study of adaptive radiation, driven by natural and sexual selection. Parasite mediated selection is an important component in these processes, and the evolution of their immune system therefore merits special attention. In this study, light is shed on the phylogeny of the b family of cichlid major histocompatibility complex (MHC) class IIB genes. Full-length coding sequences were used to reconstruct phylogenies using criteria of maximum parsimony, maximum likelihood and Bayesian inference. All analyses suggest monophyly of the b family of cichlid MHC class IIB genes, although sequences of the cichlid sister taxa are currently not available. Two evolutionary lineages of these genes, respectively encompassing the recently defined genomic regions DBB-DEB-DFB and DCB-DDB, show highly contrasting levels of differentiation. To explore putative causes for these differences, exon 2 sequences were screened for variation in recombination rate and strength of selection. The more diversified lineage of cichlid MHC class IIB b genes was found to have higher levels of both recombination and selection. This is consistent with the observation in other taxa that recombination facilitates the horizontal spread of positively selected sites across MHC loci and hence contributes to fast sequence evolution. In contrast, the lineage that showed low diversification might either be under stabilizing selection or is evolutionary constrained by its low recombination rate. We speculate whether this lineage might include MHC genes with non-classical functions.     

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.     

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     

Balancing selection, random genetic drift, and genetic variation at the major histocompatibility complex in two wild populations of guppies (Poecilia reticulata)

Our understanding of the evolution of genes of the major histocompatibility complex (MHC) is rapidly increasing, but there are still enigmatic questions remaining, particularly regarding the maintenance of high levels of MHC polymorphisms in small, isolated populations. Here, we analyze the genetic variation at eight microsatellite loci and sequence variation at exon 2 of the MHC class IIB (DAB) genes in two wild populations of the Trinidadian guppy, Poecilia reticulata. We compare the genetic variation of a small (Ne, 100) and relatively isolated upland population to that of its much larger (Ne approximately 2400) downstream counterpart. As predicted, microsatellite diversity in the upland population is significantly lower and highly differentiated from the population further downstream. Surprisingly, however, these guppy populations are not differentiated by MHC genetic variation and show very similar levels of allelic richness. Computer simulations indicate that the observed level of genetic variation can be maintained with overdominant selection acting at three DAB loci. The selection coefficients differ dramatically between the upland (s > or = 0.2) and lowland (s < or = 0.01) populations. Parasitological analysis on wild-caught fish shows that parasite load is significantly higher on upland than on lowland fish, which suggests that large differences in selection intensity may indeed exist between populations. Based on the infection intensity, a substantial proportion of the upland fish would have suffered direct or indirect fitness consequences as a result of their high parasite loads. Selection by parasites plays a particularly important role in the evolution of guppies in the upland habitat, which has resulted in high levels of MHC diversity being maintained in this population despite considerable genetic drift.     

The evolutionary ecology of the major histocompatibility complex

The major histocompatibility complex (MHC) has become a paradigm for how selection can act to maintain adaptively important genetic diversity in natural populations. Here, we review the contribution of studies on the MHC in non-model species to our understanding of how selection affects MHC diversity, emphasising how ecological and ethological processes influence the tempo and mode of evolution at the MHC, and conversely, how variability at the MHC affects individual fitness, population dynamics and viability. We focus on three main areas: the types of information that have been used to detect the action of selection on MHC genes; the relative contributions of parasite-mediated and sexual selection on the maintenance of MHC diversity; and possible future lines of research that may help resolve some of the unanswered issues associated with MHC evolution.     

Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.     

中华白海豚3 个MHC 基因座位遗传变异的初步分析

主要组织相容性复合体(Major histocompatibility complex,MHC) 基因是由一组紧密连锁的基因组成,是哺乳动物免疫系统中最重要的组成部分。本文选择3 个MHC 基因座位的第二外元,即:MHC-I 类基因和II 类基因的DRA 和DQB 座位,初步调查濒危物种中华白海豚的遗传变异。共鉴定了2 个DRA、2 个DQB 和7 MHC-I等位基因。DRA 座位遗传变异非常低,而DQB 和MHC-I 座位具有相对较高水平的遗传变异。并且,在DQB 和MHC-I 基因座位的假定的抗原结合位点(Antigen binding sites,ABS),非同义替代明显大于同义替代,提示平衡选择(Balancing selection)维持这两个座位的多态性,而在DRA 座位上,并没有检测到平衡选择。系统发生分析表明中华白海豚的MHC 等位基因没有聚在一起,而是和其他的物种聚在一起,符合MHC 跨种进化(Transspecies evolution)的模式。     

EVALUATING PATTERNS OF CONVERGENT EVOLUTION AND TRANS‐SPECIES POLYMORPHISM AT MHC IMMUNOGENES IN TWO SYMPATRIC STICKLEBACK SPECIES

The immunologically important major histocompatibility complex (MHC) harbors some of the most polymorphic genes in vertebrates. These genes presumably evolve under parasite‐mediated selection and frequently show inconsistent allelic genealogies, where some alleles are more similar between species than within species. This phenomenon is thought to arise either from convergent evolution under parallel selection or from the preservation of ancient allelic lineages beyond speciation events (trans‐species polymorphism, TSP). Here, we examine natural populations of two sympatric stickleback species (Gasterosteus aculeatus and Pungitius pungitius) to investigate the contribution of these two mechanisms to the evolution of inconsistent allelic genealogies at the MHC. Overlapping parasite taxa between the two host species in three different habitats suggest contemporary parallel selection on the MHC genes. Accordingly, we detected a lack of species‐specific phylogenetic clustering in the immunologically relevant antigen‐binding residues of the MHC IIB genes which contrasted with the rest of the coding and noncoding sequence. However, clustering was not habitat‐specific and a codon‐usage analysis revealed patterns of similarity by descent. In this light, common descent via TSP, in combination with intraspecies gene conversion, rather than convergent evolution is the more strongly supported scenario for the inconsistent genealogy at the MHC.     

MHC class II variation in the endangered European mink <Emphasis Type="Italic">Mustela lutreola</Emphasis> (L. 1761)—consequences for species conservation

The polymorphic major histocompatibility complex (MHC) has gained a specific relevance in pathogen resistance and mate choice. Particularly the antigen-binding site (ABS), encoded by exon 2 of the DRB class II gene, exhibits numerous alleles and extensive sequence variations between alleles. A lack of MHC variability has attributed to instances such as bottleneck effects or relaxed selection pressure and has a certain impact on the long-term viability of the species concerned. As a result of seriously decreased population density during the last century, the current population of the endangered European mink (Mustela lutreola, L. 1761) has suffered from geographic isolation. In this study, we amplified a partial sequence of the MHC class II DRB exon 2 (229 bp), assessed the degree of genetic variation and compared the variability with those of other Mustelidae. As a result, nine alleles were detected in 20 investigated individuals, which differ from each other by four to 25 nucleotide substitutions (two to 11 amino acid substitutions). Whilst an equal ratio for synonymous and non-synonymous substitutions was found inside the ABS, synonymous substitutions were significantly higher than non-synonymous substitutions in the non-ABS region. Results might indicate that no positive selection exists within the ex situ population of M. lutreola, at least in the analysed fragment. In addition, phylogenetic analyses support the trans-species model of evolution. Becker and Nieberg have contributed equally to this work.     

Shared evolutionary origin of major histocompatibility complex polymorphism in sympatric lemurs

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

Sequence Polymorphism and Evolution of Three Cetacean MHC Genes

Sequence variability at three major histocompatibility complex (MHC) genes (DQB, DRA, and MHC-I) of cetaceans was investigated in order to get an overall understanding of cetacean MHC evolution. Little sequence variation was detected at the DRA locus, while extensive and considerable variability were found at the MHC-I and DQB loci. Phylogenetic reconstruction and sequence comparison revealed extensive sharing of identical MHC alleles among different species at the three MHC loci examined. Comparisons of phylogenetic trees for these MHC loci with the trees reconstructed only based on non-PBR sites revealed that allelic similarity/identity possibly reflected common ancestry and were not due to adaptive convergence. At the same time, trans-species evolution was also evidenced that the allelic diversity of the three MHC loci clearly pre-dated species divergence events according to the relaxed molecular clock. It may be the forces of balancing selection acting to maintain the high sequence variability and identical alleles in trans-specific manner at the MHC-I and DQB loci.     

Characterization of the major histocompatibility complex class II genes in miiuy croaker

Major histocompatibility complex (MHC) has a central role in the adaptive immune system by presenting foreign peptide to the T-cell receptor. In order to study the molecular function and genomic characteristic of class II genes in teleost, the full lengths of MHC class IIA and IIB cDNA and genomic sequence were cloned from miiuy croaker (Miichthys miiuy). As in other teleost, four exons and three introns were identified in miiuy croaker class IIA gene; but the difference is that six exons and five introns were identified in the miiuy croaker class IIB gene. The deduced amino acid sequence of class IIA and class IIB had 26.3-85.7% and 11.0-88.8% identity with those of mammal and teleost, respectively. Real-time quantitative RT-PCR demonstrated that the MHC class IIA and IIB were ubiquitously expressed in ten normal tissues; expression levels of MHC genes were found first upregulated and then downregulated, and finally by a recovery to normal level throughout the pathogenic bacteria infection process. In addition, we report on the underlying mechanism that maintains sequences diversity among many fish species. A series of site-model tests implemented in the CODEML program revealed that positive Darwinian selection is likely the cause of the molecular evolution in the fish MHC class II genes.     

Polymorphism and signature of selection in the MHC class I genes of the three-spined stickleback Gasterosteus aculeatus

The role and intensity of positive selection maintaining the polymorphism of major histocompatibility complex (MHC) class I genes in the three-spined stickleback Gasterosteus aculeatus was investigated. The highly polymorphic set of MHC class I genes found was organized in a single linkage group. Between 5 and 14 sequence variants per individual were identified by single-stranded conformation polymorphism (SSCP) analysis. Segregation analysis studied in 10 three-spined stickleback families followed the expected pattern of Mendelian inheritance. The gamete fusion in three-spined stickleback thus seems to be random with respect to the MHC class I genes. The DNA sequence analyses showed that the expressed MHC class I loci are under strong selection pressure, possibly mediated by parasites. Codons that were revealed to be under positive selection are potentially important in antigen binding. MHC class I sequences did not form significant supported clusters within a phylogenetic tree. Analogous to MHC class II genes, it was not possible to assign the class I sequences to a specific locus, suggesting that the class I genes may have been generated by recent gene duplication.