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.     

In silico peptide-binding predictions of passerine MHC class I reveal similarities across distantly related species, suggesting convergence on the level of protein function

The major histocompatibility complex (MHC) genes are the most polymorphic genes found in the vertebrate genome, and they encode proteins that play an essential role in the adaptive immune response. Many songbirds (passerines) have been shown to have a large number of transcribed MHC class I genes compared to most mammals. To elucidate the reason for this large number of genes, we compared 14 MHC class I alleles (α1–α3 domains), from great reed warbler, house sparrow and tree sparrow, via phylogenetic analysis, homology modelling and in silico peptide-binding predictions to investigate their functional and genetic relationships. We found more pronounced clustering of the MHC class I allomorphs (allele specific proteins) in regards to their function (peptide-binding specificities) compared to their genetic relationships (amino acid sequences), indicating that the high number of alleles is of functional significance. The MHC class I allomorphs from house sparrow and tree sparrow, species that diverged 10 million years ago (MYA), had overlapping peptide-binding specificities, and these similarities across species were also confirmed in phylogenetic analyses based on amino acid sequences. Notably, there were also overlapping peptide-binding specificities in the allomorphs from house sparrow and great reed warbler, although these species diverged 30 MYA. This overlap was not found in a tree based on amino acid sequences. Our interpretation is that convergent evolution on the level of the protein function, possibly driven by selection from shared pathogens, has resulted in allomorphs with similar peptide-binding repertoires, although trans-species evolution in combination with gene conversion cannot be ruled out.     

Evolution by recombination and transspecies polymorphism in the MHC class I gene of Xenopus laevis

The patterns of major histocompatibility complex (MHC) evolution involve duplications, deletions, and independent divergence of loci during episodes punctuated by natural selection. Major differences in MHC evolution among taxa have previously been attributed to variation in linkage patterns of class I and class II MHC genes. Here we characterize patterns of evolution in the MHC class Ia gene of Xenopus laevis in terms of polymorphism, recombination, and extent of transspecies polymorphism. We also compare these patterns to see if a correlation exists with linkage or separation of the MHC class I and class II regions as seen in amphibians and teleost fishes. In X. laevis, we find high levels of polymorphism. Also, genetic exchange is relatively frequent and occurs in intron II, reshuffling allelic forms of exons 2 and 3. Evolutionary relationships among class I alleles show an intermingling of alleles from divergent Xenopus species rather than a species-specific clustering. Results indicate that the patterns of evolution are similar to those found in salmonid fishes and are different from the mode of evolution seen in primates. Similar patterns of class Ia evolution in salmonid fishes and X. laevis suggest that nonlinkage of class I and class II regions alone is insufficient to explain some patterns of MHC evolution in salmonids.     

Polymorphism, natural selection, and structural modeling of class Ia MHC in the African clawed frog (Xenopus laevis)

In the African clawed frog (Xenopus laevis), two deeply divergent allelic lineages of multiple genes of the class I MHC region have been discovered. For the MHC class I UAA locus, functional differences and the molecular basis for lineages maintenance are unknown. Alleles of linked class I region genes also exhibit strong disequilibrium with specific MHC alleles, but the underlying cause is not clear. We use MHC class Ia sequence data to estimate substitution rates and investigate structural differences between allelic lineages from protein models. Results indicate the operation of natural selection, and differences in the steric properties in the F pocket of the peptide-binding region among lineages. Variability in this pocket likely enables allelic lineages to bind very different sets of peptides and to interact differently with MHC chaperones in the endoplasmic reticulum. These results constitute evidence of the molecular evolutionary basis for 1) the maintenance of allelic lineages, 2) functional differences among lineages, and 3) strong linkage disequilibrium of allelic variants of class I region genes in X. laevis.Electronic Supplementary Material Supplementary material is available for this article at     

Polymorphisms and tissue expression of the feline leukocyte antigen class I loci FLAI-E, FLAI-H, and FLAI-K

Cytotoxic CD8+ T-cell immunosurveillance for intracellular pathogens, such as viruses, is controlled by classical major histocompatibility complex (MHC) class Ia molecules, and ideally, these antiviral T-cell populations are defined by the specific peptide and restricting MHC allele. Surprisingly, despite the utility of the cat in modeling human viral immunity, little is known about the feline leukocyte antigen class I complex (FLAI). Only a few coding sequences with uncertain locus origin and expression patterns have been reported. Of 19 class I genes, three loci—FLAI-E, FLAI-H, and FLAI-K—are predicted to encode classical molecules, and our objective was to evaluate their status by analyzing polymorphisms and tissue expression. Using locus-specific, PCR-based genotyping, we amplified 33 FLAI-E, FLAI-H, and FLAI-K alleles from 12 cats of various breeds, identifying, for the first time, alleles across three distinct loci in a feline species. Alleles shared the expected polymorphic and invariant sites in the α1/α2 domains, and full-length cDNA clones possessed all characteristic class Ia exons. Alleles could be assigned to a specific locus with reasonable confidence, although there was evidence of potentially confounding interlocus recombination between FLAI-E and FLAI-K. Only FLAI-E, FLAI-H, and FLAI-K origin alleles were amplified from cDNAs of multiple tissue types. We also defined hypervariable regions across these genes, which permitted the assignment of names to both novel and established alleles. As predicted, FLAI-E, FLAI-H, and FLAI-K fulfill the major criteria of class Ia genes. These data represent a necessary prerequisite for studying epitope-specific antiviral CD8+ T-cell responses in cats.     

Ancestral organization of the MHC revealed in the amphibian Xenopus

With the advent of the Xenopus tropicalis genome project, we analyzed scaffolds containing MHC genes. On eight scaffolds encompassing 3.65 Mbp, 122 MHC genes were found of which 110 genes were annotated. Expressed sequence tag database screening showed that most of these genes are expressed. In the extended class II and class III regions the genomic organization, excluding several block inversions, is remarkably similar to that of the human MHC. Genes in the human extended class I region are also well conserved in Xenopus, excluding the class I genes themselves. As expected from previous work on the Xenopus MHC, the single classical class I gene is tightly linked to immunoproteasome and transporter genes, defining the true class I region, present in all nonmammalian jawed vertebrates studied to date. Surprisingly, the immunoproteasome gene PSMB10 is found in the class III region rather than in the class I region, likely reflecting the ancestral condition. Xenopus DMalpha, DMbeta, and C2 genes were identified, which are not present or not clearly identifiable in the genomes of any teleosts. Of great interest are novel V-type Ig superfamily (Igsf) genes in the class III region, some of which have inhibitory motifs (ITIM) in their cytoplasmic domains. Our analysis indicates that the vertebrate MHC experienced a vigorous rearrangement in the bony fish and bird lineages, and a translocation and expansion of the class I genes in the mammalian lineage. Thus, the amphibian MHC is the most evolutionary conserved MHC so far analyzed.     

Is a mutator analogous to the Ig hypermutator of the sheep ileal Peyer's patch active on MHC class I genes in the germ line?

 Polymorphic sequence variation in the peptide-binding domains of MHC class I molecules appears to have been driven largely by the constructive action of natural selection on the specificity of the peptide-binding groove. Similar features are displayed by the variable domains of immunoglobulins generated in the sheep ileal Peyer's patch, but in this case there is evidence that the action of a targeted hypermutator acting on a selected substrate rather than antigen-driven selection is responsible for the pattern of variation in the system. Such a hypermutator acting in the germ line would not only mimic the action of natural selection but also, by convergent mutation, generate similar patterns of variation in unrelated alleles that could be interpreted as evidence for short-tract gene conversion. We analyzed human class I MHC alleles in the light of these data, but failed to find evidence of the action of a similar hypermutator. A search for other mutationally driven patterns of variation also failed, even in hypervariable residues from parsimonious phylogenies. Single-nucleotide variation at these residues is also frequent in recent allelic variants, but the data are as consistent with short-tract gene conversion as with base mutation. We conclude that the patterns of allelic variation in MHC molecules are not driven by mutational pressure, but rather by conventional mutational processes, accompanied by short-tract gene conversion and intense natural selection. Received: 6 October 1999 / Revised: 30 December 1999     

Classical MHC class I genes composed of highly divergent sequence lineages share a single locus in rainbow trout (Oncorhynchus mykiss).

The classical MHC class I genes have been known to be highly polymorphic in various vertebrates. To date, putative allelic sequences of the classical MHC class I genes in teleost fish have been reported in several studies. However, the establishment of their allelic status has been hampered in most cases by the lack of appropriate genomic information. In the present study, using heterozygous and homozygous fish, we obtained classical-type MHC class I sequences of rainbow trout (Oncorhynchus mykiss) and investigated their allelic relationship by gene amplification and Southern and Northern hybridization analyses. The results indicated that all MHC class I sequences we obtained were derived from a single locus. Based on this, a unique polymorphic nature of the MHC class I locus of rainbow trout has been revealed. The mosaic combination of highly divergent ancient sequences in the peptide-binding domains is notable, and the variable nature around the boundary between the alpha3 and transmembrane domains is unprecedented.     

Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline Mhc-DRB peptide-binding region

The long-term evolution of major histocompatibility complex (MHC) involves the birth-and-death process and independent divergence of loci during episodes punctuated by natural selection. Here, we investigated the molecular signatures of natural selection at exon-2 of MHC class II DRB gene which includes a part of the peptide-binding region (PBR) in seven of eight putative extant Felidae lineages. The DRB alleles in felids can be mainly divided into five lineages. Signatures of trans-species polymorphism among major allelic lineages indicate that balancing selection has maintained the MHC polymorphism for a long evolutionary time. Analysis based on maximum likelihood models of codon substitution revealed overall purifying selection acting on the feline DRB. Sites that have undergone positive selection and those that are under divergent selective pressure among lineages were detected and found to fall within the putative PBR. This study increased our understanding of the nature of selective forces acting on DRB during feline radiation.     

Crystal structures of two rat MHC class Ia (RT1-A) molecules that are associated differentially with peptide transporter alleles TAP-A and TAP-B

Antigenic peptides are loaded onto class I MHC molecules in the endoplasmic reticulum (ER) by a complex consisting of the MHC class I heavy chain, beta(2)-microglobulin, calreticulin, tapasin, Erp57 (ER60) and the transporter associated with antigen processing (TAP). While most mammalian species transport these peptides into the ER via a single allele of TAP, rats have evolved different TAPs, TAP-A and TAP-B, that are present in different inbred strains. Each TAP delivers a different spectrum of peptides and is associated genetically with distinct subsets of MHC class Ia alleles, but the molecular basis for the conservation (or co-evolution) of the two transporter alleles is unknown. We have determined the crystal structures of a representative of each MHC subset, viz RT1-A(a) and RT1-A1(c), in association with high-affinity nonamer peptides. The structures reveal how the chemical properties of the two different rat MHC F-pockets match those of the corresponding C termini of the peptides, corroborating biochemical data on the rates of peptide-MHC complex assembly. An unusual sequence in RT1-A1(c) leads to a major deviation from the highly conserved beta(3)/alpha(1) loop (residues 40-59) conformation in mouse and human MHC class I structures. This loop change contributes to profound changes in the shape of the A-pocket in the peptide-binding groove and may explain the function of RT1-A1(c) as an inhibitory natural killer cell ligand.     

A tangled history: patterns of major histocompatibility complex evolution in the African mole-rats (Family: Bathyergidae)

Phylogenetic trees based upon major histocompatibility complex (MHC) gene sequences, particularly those encompassing sites encoding the antigen recognition site, are often discordant with the species tree. It has been argued that the principal cause of such discordance is the presence of ancestrally derived polymorphisms persisting through speciation events as a consequence of selection. In the present study, we examine the evolution of the MHC class II DQα1 gene in an unusual family of hystricomorph rodents, the African mole-rats (Family: Bathyergidae). We show that there is a high level of trans-species polymorphism and that this is a result of positive selection. Furthermore, the major lineages of the gene tree are characterized by allelic motifs occurring in regions that coincide with the pocket domains of the putative antigen recognition site, a region that has been shown to be under positive selection in a number of MHC genes from a range of species. Finally, these alleles may have been retained for at least 48 million years. This is significantly older than the estimate for the equivalent primate locus and appears to be one of the oldest documented sets of MHC alleles. We suggest that these allelic motifs possess polymorphisms that have been immunologically important to African mole-rats over long periods of evolutionary history.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 493–503.     

ERp57 interacts with conserved cysteine residues in the MHC class I peptide-binding groove

The oxidoreductase ERp57 is a component of the major histocompatibility complex (MHC) class I peptide-loading complex. ERp57 can interact directly with MHC class I molecules, however, little is known about which of the cysteine residues within the MHC class I molecule are relevant to this interaction. MHC class I molecules possess conserved disulfide bonds between cysteines 101-164, and 203-259 in the peptide-binding and alpha3 domain, respectively. By studying a series of mutants of these conserved residues, we demonstrate that ERp57 predominantly associates with cysteine residues in the peptide-binding domain, thus indicating ERp57 has direct access to the peptide-binding groove of MHC class I molecules during assembly.     

Unusual evolutionary conservation and further species-specific adaptations of a large family of nonclassical MHC class Ib genes across different degrees of genome ploidy in the amphibian subfamily Xenopodinae

Nonclassical MHC class Ib (class Ib) genes are a family of highly diverse and rapidly evolving genes wherein gene numbers, organization, and expression markedly differ even among closely related species rendering class Ib phylogeny difficult to establish. Whereas among mammals there are few unambiguous class Ib gene orthologs, different amphibian species belonging to the anuran subfamily Xenopodinae exhibit an unusually high degree of conservation among multiple class Ib gene lineages. Comparative genomic analysis of class Ib gene loci of two divergent (~65 million years) Xenopodinae subfamily members Xenopus laevis (allotetraploid) and Xenopus tropicalis (diploid) shows that both species possess a large cluster of class Ib genes denoted as Xenopus/Silurana nonclassical (XNC/SNC). Our study reveals two distinct phylogenetic patterns among these genes: some gene lineages display a high degree of flexibility, as demonstrated by species-specific expansion and contractions, whereas other class Ib gene lineages have been maintained as monogenic subfamilies with very few changes in their nucleotide sequence across divergent species. In this second category, we further investigated the XNC/SNC10 gene lineage that in X. laevis is required for the development of a distinct semi-invariant T cell population. We report compelling evidence of the remarkable high degree of conservation of this gene lineage that is present in all 12 species of the Xenopodinae examined, including species with different degrees of ploidy ranging from 2, 4, 8 to 12 N. This suggests that the critical role of XNC10 during early T cell development is conserved in amphibians.     

Genetic variability of MHC class II DQB exon 2 alleles in yak (Bos grunniens)

The major histocompatibility class (MHC) DQ molecules are dimeric glycoproteins revealing antigen presentation to CD⁴⁺ T cells. In the present study, the exon 2 of the MHC class II DQB gene from 32 yaks (Bos grunniens) was cloned, sequenced and compared with previously reported patterns for other bovidae. It was revealed by sequence analyses that there are 25 DQB exon 2 alleles among 32 yaks, all alleles are found to belong to DQB1 loci. These alleles exhibited a high degree of nucleotide and amino acid polymorphisms with most amino acid variations occurring at positions forming the peptide-binding sites. The DQB loci were analyzed for patterns of synonymous (d _S) and non-synonymous (d _N) substitution. The yak was observed to be under strong positive selection in the DQB exon 2 peptide-binding sites (d _N = 0.15, P < 0.001). It appears that this variability among yaks confers the ability to mount immune responses to a wide variety of peptides or pathogens.     

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.     

Structure, regulatory polymorphisms, and allelic hypervariability regions in murine I-A alpha

Class II major histocompatibility complex (MHC) molecules, the Ia antigens, are intimately involved in regulating the intensity and specificity of the cellular and humoral responses to T cell-dependent antigens. One approach to understanding the mechanism of this regulation is to analyze the structure and allelic polymorphism of Ia molecules. In addition there are regulatory polymorphisms in the expression of the I-E alpha and I-E beta class II MHC polypeptide chains. Analysis of the cDNA sequence indicates that I-A and I-E alpha chains are similar with short stretches of homology and other regions of nonhomology. Analysis of Northern blots of mRNA indicates that at least three separate types of regulatory polymorphisms result in failure of expression of I-E alpha. Comparison of allelic sequences of six alleles of the I-A alpha chain shows that almost all of the allelic polymorphism is in the first domain and that within the first domain it is clustered in three allelic hypervariable regions within the first domain of I-A alpha. The structural and functional implications of these findings are discussed.     

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.