The salmonid MHC class I: more ancient loci uncovered

An unprecedented level of sequence diversity has been maintained in the salmonid major histocompatibility complex (MHC) class I UBA gene, with between lineage AA sequence identities as low as 34%. The derivation of deep allelic lineages may have occurred through interlocus exon shuffling or convergence of ancient loci with the UBA locus, but until recently, no such ancient loci were uncovered. Herein, we document the existence of eight additional MHC class I loci in salmon (UCA, UDA, UEA, UFA, UGA, UHA, ULA, and ZE), six of which share exon 2 and 3 lineages with UBA, and three of which have not been described elsewhere. Half of the UBA exon 2 lineages and all UBA exon 3 lineages are shared with other loci. Two loci, UGA and UEA, share only a single exon lineage with UBA, likely generated through exon shuffling. Based on sequence homologies, we hypothesize that most exchanges and duplications occurred before or during tetraploidization (50 to 100 Ma). Novel loci that share no relationship with other salmonid loci are also identified (UHA and ZE). Each locus is evaluated for its potential to function as a class Ia gene based on gene expression, conserved residues and polymorphism. UBA is the only locus that can indisputably be classified as a class Ia gene, although three of the eight loci (ZE, UCA, and ULA) conform in three out of four measures. We hypothesize that these additional loci are in varying states of degradation to class Ib genes.     

Polymorphisms in major histocompatibility complex class IIα genes are associated with resistance to infectious hematopoietic necrosis in rainbow trout,Oncorhynchus mykiss (Walbaum, 1792)

Infectious hematopoietic necrosis is a serious viral disease of salmonids, including rainbow trout Oncorhynchus mykiss, and causes tremendous economic losses to the rainbow trout farming industry. Major histocompatibility complex (MHC) genes are crucial elements of adaptive immunity in vertebrate organisms and have been linked with the resistance to numerous pathogenic diseases. In this study, polymerase chain reaction‐single strand conformation polymorphism (PCR‐SSCP) followed by cloning and sequencing were used to examine polymorphisms in the DAA genes (specifically DAA exon 2 of MHC class IIα) of rainbow trout and investigate their association with the infectious hematopoietic necrosis virus (IHNV) resistance in rainbow trout. Seventeen alleles were resolved, including 13 novel alleles. Individuals possessed between two and five alleles, indicating that the genome harbours at least three closely‐related DAA exon 2 loci. The ratio of non‐synonymous to synonymous nucleotide substitutions suggested that DAA exon 2 is under positive selection. A greater variability of amino acids and non‐synonymous nucleotide substitution rate was evident in the peptide‐binding region (PBR) than in the non‐PBR (27.75%). Importantly, the analyses revealed that certain MHC class IIα alleles appear to confer resistance to IHNV in rainbow trout, while others confer susceptibility. The most common alleles in the resistant populations of rainbow trout, Onmy‐DAA*1301 and Onmy‐DAA*0304, confer resistance to IHNV and were not present in the susceptible population. Hence, these alleles may be ideal molecular markers that can assist the breeding of IHNV resistance in rainbow trout.     

Analysis of MHC class I genes across horse MHC haplotypes

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.     

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.     

Sequence analysis of MHC class I α2 from sockeye salmon (Oncorhynchus nerka)

Most studies assessing adaptive MHC diversity in salmon populations have focused on the classical class II DAB or DAA loci, as these have been most amenable to single PCR amplifications due to their relatively low level of sequence divergence. Herein, we report the characterization of the classical class I UBA α2 locus based on collections taken throughout the species range of sockeye salmon (Oncorhynchus nerka). Through use of multiple lineage-specific primer sets, denaturing gradient gel electrophoresis and sequencing, we identified thirty-four alleles from three highly divergent lineages. Sequence identity between lineages ranged from 30.0% to 56.8% but was relatively high within lineages. Allelic identity within the antigen recognition site (ARS) was greater than for the longer sequence. Global positive selection on UBA was seen at the sequence level (dN:dS = 1.012) with four codons under positive selection and 12 codons under negative selection.     

MHC polymorphism and disease resistance to Vibrio anguillarum in 12 selective Japanese flounder (Paralichthys olivaceus) families

Genetic variation in the major histocompatibility complex (MHC) class IIB was tested in Japanese flounder (Paralichthys olivaceus) for survival after challenge with bacterial infection. The material consisted of 6000 Japanese flounder from 60 families challenged with Vibrio anguillarum, which causes significantly different mortality in flounder families. Five individuals from each of six high-resistance (HR) and six low-resistance (LR) families were screened for their MHC class IIB genotypes using sequence analysis. High polymorphism of MHC IIB gene and at least three loci were discovered in Japanese flounder and the rate of d(N) occurred at a significantly higher frequency than that of d(S) in PBR. Among 60 individuals, 76 alleles were discovered and 15 alleles were used to study associations between alleles and resistance to disease. We found highly significant associations between resistance towards infectious disease caused by V. anguillarum and MHC class IIB polymorphism in Japanese flounder. Some alleles appeared in both HR and LR families, while some alleles were only discovered in HR or LR families. One allele, Paol-DAB*4301, was significantly more prevalent in HR families than in LR families (P=0.023). Paol-DAB*0601, Paol-DAB*0801, Paol-DAB*2001, Paol-DAB*3803 were discovered in two HR families with high frequency. One allele, Paol-DAB*1601, was discovered in three LR families. The steady heredity of MHC class IIB alleles was observed, and the family having Paol-DAB*4301 alleles was confirmed with higher resistance to V. anguillarum. This study confirmed the association between alleles of MHC class IIB gene and disease resistance, and also detected some alleles which might be correlated with high bacterial infection resistance. The disease resistance-related MHC markers could be used for molecular marker-assisted selective breeding in the flounder.     

Patterns of variation of the major histocompatibility complex class IIB loci in Chinese goose (<Emphasis Type="Italic">Anser Cygnoides</Emphasis>)

In order to understand the variations of genomic organization of the major histocompatibility complex (MHC) and provide data for the studies on disease resistance of avian species, the MHC class II polymorphism in Chinese Z-goose was investigated for the first time in the present study. Eight alleles, which probably came from different loci, were found in six different geese with only one obvious band in the restriction fragment length polymorphism data. The numbers of nonsynonymous substitutions (dN) in peptide binding region of exon 2 were higher than that of synonymous substitutions (dS), and no stop codons or frameshift mutations were found in this region, indicating that balance selection was in operation, and the sequences are not likely to be pseudogenes. In addition, we successfully obtained five different long MHC class II fragments (about 1,162 bp) in six geese and found that the length of intron 1 was longer than that in chicken and some other birds, but intron 2 seemed to be intermediate in length. The phylogenetic tree appeared to branch in an order consistent with accepted evolutionary pathway. X. Zhou and C. Li contributed equally to this work and shared the co-first authorship.     

HLA-DRA promoter polymorphism and diversity generation within the immune system

Multiple major histocompatibility complex (MHC) alleles exist at most class I and II loci. Polymorphism of MHC polypeptides may reflect either different levels of selective pressure operating on each molecule or different mutation rates at different loci. To gain further insight into this issue, we sequenced the non-coding promoter region of the HLA-DRA gene from several Epstein-Barr virus-transformed B cell lines and compared the extent of polymorphism found in this region with the known polymorphism of the HLA-DQB promoter. Our results indicate that the HLA-DRA promoter displays a low level of polymorphism while the promoter of HLA-DQB exhibits a nucleotide substitution rate fivefold greater than that of DRA. Moreover, through phylogenetic analysis, the HLA-DRA promoter was found to have diverged much less than the associated alleles of HLA-DRB1 and -DQA1. Taken together, these results suggest that the HLA-DRA promoter is highly conserved and may be under a stronger functional constraint than the promoter regions of other MHC class II genes. Received: 30 May 1996 / Revised: 24 January 1997     

Physical and genetic mapping of the rainbow trout major histocompatibility regions: evidence for duplication of the class I region

One of the most unexpected discoveries in MHC genetics came from studies dealing with the teleost MHC. Initially discovered in zebrafish, the MHC class I and II regions of all bony fish are not linked. Previous segregation analysis in trout suggested that the class I and II regions reside on completely different chromosomes. To learn more about MHC genomics in trout, we have isolated BAC clones harboring class Ia and Ib loci, a single BAC clone containing an MH class II gene ( DAB), as well as BAC clones containing the ABCB2 gene. Upon PCR and sequence confirmation, BAC clones were labeled and used as probes for in situ hybridization on rainbow trout metaphase chromosomes for determination of the physical locations of the trout MH regions. Finally, SNPs, RFLPs, and microsatellites found within the BAC clones allowed for these regions to be assigned to specific linkage groups on the OSU x Hotcreek (HC) and OSU x Arlee (ARL) genetic linkage maps. Our data demonstrate that the trout MH regions are located on at least four different chromosomes and the corresponding linkage groups, while also providing direct evidence for the partial duplication of the MH class I region in trout.     

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.     

Exchanges of short polymorphic DNA segments predating speciation in feline major histocompatibility complex class I genes

Sequence comparisons of 14 distinct MHC class I cDNA clones isolated from species representing the three major taxonomic lineages of Felidae (domestic cat lineage, ocelot lineage, and pantherine lineage) revealed that feline MHC class I alleles have highly mosaic structures with short polymorphic sequence motifs that are rearranged between alleles of individual MHC loci, between MHC class I genes within cat species, and between homologous MHC loci in different species. The pattern of sequence variation in felids supports the role of the following factors in production and maintenance of MHC variation: (1) gradual spontaneous mutation; (2) selective pressure to conserve certain residues but also to vary in hypervariable regions, notably residues that functionally participate in antigen recognition and presentation; and (3) recombination-mediated gene exchange between alleles and between related genes. The overall amount of genetic variation observed among MHC class I genes in the Felidae family is no greater than the amount of variation within any outbred cat species (i.e., domestic cat, ocelot). The occurrence of equivalent levels of polymorphism plus the simultaneous persistence of the same sequence motifs in divergent feline species suggest that most MHC class I nucleotide site polymorphism predated species divergences. Ancient polymorphisms have been transmitted through the speciation events and modern feline MHC class I alleles were derived by recombinational exchange of polymorphic sequence motifs. Moreover, some of these sequence motifs were found in other mammalian MHC class I genes, such as classical human HLA-B5, nonclassical human HLA-E class I genes, and bovine class I genes. These results raise the prospect of an ancient origin for some motifs, although the possibility of convergence in parallel mammalian radiations cannot be excluded. Correspondence to: N. Yuhki     

Genetic variation of the major histocompatibility complex (MHC class II β gene) in the threatened Gila trout, <Emphasis Type="Italic">Oncorhynchus gilae gilae</Emphasis>

Gila trout (Oncorhynchus gilae gilae) was federally protected in 1973 because of severe declines in abundance and geographic range size. At present, four relict genetic lineages of the species remain in mountain streams of New Mexico and Arizona, USA. Management actions aimed at species recovery, including hatchery production and restocking of formerly occupied streams, have been guided by information from non-functional genetic markers. In this study, we investigated genetic variation at exon 2 of the major histocompatibility complex (MHC) class II β gene that is involved in pathogen resistance and thus presumably under natural selection. Phylogenetic analysis revealed trans-species polymorphism and a significantly high ratio of non-synonymous to synonymous amino acid changes consistent with the action of historical balancing selection that maintained diversity at this locus in the past. However, Gila trout exhibited low allelic diversity (five alleles from 142 individuals assayed) compared to some other salmonid fishes, and populations that originated exclusively from hatcheries possessed three or fewer MHC alleles. Comparative analysis of genetic variation at MHC and six presumably neutrally evolving microsatellite loci revealed that genetic drift cannot be rejected as a primary force governing evolution of MHC in contemporary populations of Gila trout. Maintenance of diversity at MHC will require careful implementation of hatchery breeding protocols and continued protection of wild populations to prevent loss of allelic diversity due to drift.     

The major histocompatibility class I locus in Atlantic salmon (Salmo salar L.): polymorphism,linkage analysis and protein modelling

A cDNA library screening using the conserved exon 4 of Atlantic salmon Mhc class I as probe provided the basis for a study on Mhc class I polymorphism in a breeding population. Twelve different alleles were identified in the 82 dams and sires studied. No individual expressed more than two alleles, which corresponded to the diploid segregation patterns of the polymorphic marker residing within the 3'-untranslated tail. Close linkage between the Sasa-UBA and Sasa-TAP2B loci strengthens the claim that Sasa-UBA is the major Mhc class I locus in Atlantic salmon. We found no evidence for a second expressed classical or non-classical Mhc class I locus in Atlantic salmon. A phylogenetic analysis of salmonid Mhc class I sequences showed domains conserved between rainbow trout, brown trout and Atlantic salmon. Evidence for shuffling of the alpha(1) domain was identified and lineages of the remaining alpha(2) through the cytoplasmic tail gene segment can be defined. The coding sequence of one allele was found associated with two different markers, suggesting recombination within the 3'-tail dinucleotide repeat itself. Protein modelling of several Sasa-UBA alleles shows distinct differences in their peptide binding domains and enables a further understanding of the functionality of the high polymorphism.     

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.     

Class I MHC expression in the yellow baboon

MHC class I molecules play a crucial role in the immune response to pathogens and vaccines and in self/non-self recognition. Therefore, characterization of MHC class I gene expression of Papio subspecies is a prerequisite for studies of immunology and transplantation in the baboon (papio hamadryas). To elucidate MHC class I expression and variation within Papio subspecies and to further investigate the evolution of A and B loci in Old World primates, we have characterized the expressed class I repertoire of the yellow baboon (Papio hamadryas cynocephalus) by cDNA library screening. A total of nine distinct MHC class I cDNAs were isolated from a spleen cDNA library. The four A alleles and four B alleles obtained represent four distinct loci indicating that a duplication of the A and B loci has taken place in the lineage leading to these Old World primates. No HLA--C homologue/orthologue was found. In addition a single, nonclassical homologue of HLA--E was characterized. Examination of nucleotide and extrapolated protein sequences indicates that alleles at the two B loci are much more diversified than the alleles at the A loci. One of the A loci in particular appears to display very limited polymorphism in both Papio hamadryas cynocephalus and Papio hamadryas anubis subspecies. The failure to detect a homologue of HLA--C in the baboon provides additional evidence for the more recent origin of this locus in the pongidae and hominidae: Further comparative analysis with MHC sequences among the primate species reveals specific patterns of divergence and conservation within class I molecules of the yellow baboon.