Major histocompatibility complex variation and evolution at a single, expressed DQA locus in two genera of elephants

Genes of the vertebrate major histocompatibility complex (MHC) are crucial to defense against infectious disease, provide an important measure of functional genetic diversity, and have been implicated in mate choice and kin recognition. As a result, MHC loci have been characterized for a number of vertebrate species, especially mammals; however, elephants are a notable exception. Our study is the first to characterize patterns of genetic diversity and natural selection in the elephant MHC. We did so using DNA sequences from a single, expressed DQA locus in elephants. We characterized six alleles in 30 African elephants (Loxodonta africana) and four alleles in three Asian elephants (Elephas maximus). In addition, for two of the African alleles and three of the Asian alleles, we characterized complete coding sequences (exons 1–5) and nearly complete non-coding sequences (introns 2–4) for the class II DQA loci. Compared to DQA in other wild mammals, we found moderate polymorphism and allelic diversity and similar patterns of selection; patterns of non-synonymous and synonymous substitutions were consistent with balancing selection acting on the peptides involved in antigen binding in the second exon. In addition, balancing selection has led to strong trans-species allelism that has maintained multiple allelic lineages across both genera of extant elephants for at least 6 million years. We discuss our results in the context of MHC diversity in other mammals and patterns of evolution in elephants.     

Characterization of the major histocompatibility complex class II <Emphasis Type="Italic">DOB,DPB1</Emphasis>, and <Emphasis Type="Italic">DQB1</Emphasis> alleles in cynomolgus macaques of Vietnamese origin

Major histocompatibility complex (MHC) molecules play an important role in the susceptibility and/or resistance to many diseases. To gain an insight into the MHC background and to facilitate the experimental use of cynomolgus macaques, the second exon of the MhcMafa-DOB, -DPB1, and -DQB1 genes from 143 cynomolgus macaques were characterized by cloning to sequencing. A total of 16 Mafa-DOB, 16 Mafa-DPB1, and 34 Mafa-DQB1 alleles were identified, which revealed limited, moderate, and marked allelic polymorphism at DOB, DPB1, and DQB1, respectively, in a cohort of cynomolgus macaques of Vietnamese origin. In addition, 16 Mafa-DOB, 5 Mafa-DPB1, and 8 Mafa-DQB1 alleles represented novel sequences that had not been reported in earlier studies. Almost of the sequences detected at the DOB and DQB1 locus in the present study belonged to DOB*01 (100%) and DQB1*06 (62%) lineages, respectively. Interestingly, four, three, and one high-frequency alleles were detected at Mafa-DOB, -DPB1, and -DQB1, respectively, in this monkeys. The alleles with the highest frequency among these monkeys were Mafa-DOB*010102, Mafa-DPB1*13, and Mafa-DQB1*0616, and these were found in 33 (25.6%) of 129 monkeys, 32 (31.37%) of 102 monkeys, and 30 (31%) of 143 monkeys, respectively. The high-frequency alleles may represent high priority targets for additional characterization of immune function. We also carried out evolutionary and population analyses using these sequences to reveal population-specific alleles. This information will not only promote the understanding of MHC diversity and polymorphism in the cynomolgus macaque but will also increase the value of this species as a model for biomedical research.     

Recent divergence of the<Emphasis Type="Italic"> HLA-DRB1*04</Emphasis> allelic lineage from the<Emphasis Type="Italic"> DRB1*0701</Emphasis> lineage after the separation of the human and chimpanzee species

Conventional phylogenetic trees for the human leukocyte antigen (HLA)-DRB1 alleles constructed by the neighbor-joining (Saitou and Nei 1987) and UPGMA (Sneath and Sokal 1973) methods using nucleotide sequences of the DRB1 alleles suggest that DRB1*0701 may have diverged from other DRB1 alleles before the separation of the human and chimpanzee species, because of a large number of nucleotide changes in DRB1*0701 compared with any of the other DRB1 alleles. Here we show new evidence that the haplotypes centering on DRB1*0701 and DRB1*04 alleles are the most homologous. This suggests that these haplotypes have derived from the common ancestral haplotype, and that they have likely retained complete linkage disequilibrium even after the divergence of the DRB1*0701 and DRB1*04 allelic lineages. Together with the corresponding haplotype carrying chimpanzee DRB1*0701, which has a high sequence homology to HLA-DRB1*0701, these haplotypes reveal that: (1) the DRB1*04 allelic lineage may have been generated from the DRB1*0701 lineage after the separation of the human and chimpanzee species; (2) the DRB1*04 allelic lineage possibly has a higher substitution rate of DRB1 compared with pseudogene and neutral region; (3) there could be a significant difference in the substitution rate of DRB1 between the DRB1*0701 and DRB1*04 allelic lineages. Based on the difference between the present and previous results, we would like to propose that phylogenetic studies using not only nucleotide sequences of the DRB1 alleles but also haplotypes centering on the alleles should be conducted for understanding detailed phylogenetic relationships of the DRB1 alleles.     

Allelic differences in hemolytic activity and protein concentration of BF molecules are found in association with particular HLA haplotypes

After separating the *F and *S alleles by electrophoresis the allele-specific hemolytic activity was detected by agarose overlay method using the programmable densitometer for scanning. The hemolytic activity of BF allotypes was analyzed from 81 individuals. In thirteen FS heterozygous serum samples BF F had lower hemolysis than BF S. Four FF homozygous samples also exhibited lower hemolysis than a homozygous control sample. The low hemolytic activity of F in FS heterozygotes was not due to decreased protein concentrations relative to S. On the contrary, BF F was associated with higher protein concentration than BF S. The relative quantitation of the allele specific BF protein was done by crossed immunoelectrophoresis. BF F with low hemolytic activity but with high protein concentration associated strongly with HLA B35 phenotype and the family material confirmed the association with the haplotypes A3, Cw4, B35, DR1, BFFB, C4A3BQO (or A2BQO, A3,2BQO). The results suggest that particular MHC haplotypes contain a factor B allele with encoding for poor hemolytic activity or that MHC haplotype specific regulatory elements affect pre- or post-translational activity levels.     

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     

The role of major histocompatibility complex diversity in vigour of fish males (Abramis brama L.) and parasite selection

The major histocompatibility complex (MHC) presents a group of genes with highly polymorphic loci involved in specific immune responses. The factors maintaining extensive MHC polymorphism have been questioned, considering three possible hypotheses of parasite‐mediated selection driving an extensive MHC diversity (i.e. heterozygote advantage, rare‐allele advantage, and favouring optimal MHC diversity). The patterns of MHC diversity of class IIB genes were investigated following two noncontradicting hypotheses, parasite‐driven selection and MHC‐based mating preferences, using males of common bream collected in the spawning period. Two allelic groups DAB1 and DAB3 were recognized from the phylogenetic analyses. Individuals expressed one or two alleles of the same or different allelic groups. Several individuals shared identical alleles; however, the presence of parasite species was not associated with the occurrence of a particular allele. The presence of different allelic groups (only DAB1, only DAB3, or both DAB1 and DAB3) in individuals was not associated with parasite presence or diversity. The expression of two DAB1 alleles was associated with higher endoparasite abundance. Moreover, nucleotide diversity in individuals expressing a single type of alleles (DAB1 or DAB3) increased with the abundance of ectoparasitic Dactylogyrus spp. (Monogenea) and Ergasilus sp. (Crustacea). This suggests that the expression of two alleles of a single allelic type is related to high metazoan parasite infection whereas no significant influence of parasitism on the combined allelic form (the presence of both DAB1 and DAB3 alleles) was found. Moreover, the expression of two alleles of a single allelic type was related to decreased immunocompetence measured by spleen size. The condition factor was higher in fish expressing the combined allelic type. Thus, the presence of alleles of different lineages in individuals appears to be advantageous for individual male fitness. The expression of a single allelic type was related to higher sexual ornamentation, which could support the role of MHC in the hypothesis of the sexual selection of ‘good genes’. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society, 2007, 90 , 525–538.     

2004 Nomenclature for the chicken major histocompatibility (<Emphasis Type="Italic">B</Emphasis> and <Emphasis Type="Italic">Y</Emphasis> ) complex

The first standard nomenclature for the chicken (Gallus gallus) major histocompatibility (B) complex published in 1982 describing chicken major histocompatibility complex (MHC) variability is being revised to include subsequent findings. Considerable progress has been made in identifying the genes that define this polymorphic region. Allelic sequences for MHC genes are accumulating at an increasing rate without a standard system of nomenclature in place. The recommendations presented here were derived in workshops held during International Society of Animal Genetics and Avian Immunology Research Group meetings. A nomenclature for B and Y (Rfp-Y) loci and alleles has been developed that can be applied to existing and newly defined haplotypes including recombinants. A list of the current standard B haplotypes is provided with reference stock, allele designations, and GenBank numbers for corresponding MHC class I and class II sequences. An updated list of proposed names for B recombinant haplotypes is included, as well as a list of over 17 Y haplotypes designated to date.     

Large deletions within the first intron in <Emphasis Type="Italic">VRN-1</Emphasis> are associated with spring growth habit in barley and wheat

The broad adaptability of wheat and barley is in part attributable to their flexible growth habit, in that spring forms have recurrently evolved from the ancestral winter growth habit. In diploid wheat and barley growth habit is determined by allelic variation at the VRN-1 and/or VRN-2 loci, whereas in the polyploid wheat species it is determined primarily by allelic variation at VRN-1. Dominant Vrn-A1 alleles for spring growth habit are frequently associated with mutations in the promoter region in diploid wheat and in the A genome of common wheat. However, several dominant Vrn-A1, Vrn-B1, Vrn-D1 (common wheat) and Vrn-H1 (barley) alleles show no polymorphisms in the promoter region relative to their respective recessive alleles. In this study, we sequenced the complete VRN-1 gene from these accessions and found that all of them have large deletions within the first intron, which overlap in a 4-kb region. Furthermore, a 2.8-kb segment within the 4-kb region showed high sequence conservation among the different recessive alleles. PCR markers for these deletions showed that similar deletions were present in all the accessions with known Vrn-B1 and Vrn-D1 alleles, and in 51 hexaploid spring wheat accessions previously shown to have no polymorphisms in the VRN-A1 promoter region. Twenty-four tetraploid wheat accessions had a similar deletion in VRN-A1 intron 1. We hypothesize that the 2.8-kb conserved region includes regulatory elements important for the vernalization requirement. Epistatic interactions between VRN-H2 and the VRN-H1 allele with the intron 1 deletion suggest that the deleted region may include a recognition site for the flowering repression mediated by the product of the VRN-H2 gene of barley.     

Locus specificity of polymorphic alleles and evolution by a birth-and-death process in mammalian MHC genes

We have conducted an extensive phylogenetic analysis of polymorphic alleles from human and mouse major histocompatibility complex (MHC) class I and class II genes. The phylogenetic tree obtained for 212 complete human class I allele sequences (HLA-A, -B, and -C) has shown that all alleles from the same locus form a single cluster, which is highly supported by bootstrap values, except for one HLA-B allele (HLA-B*7301). Mouse MHC class I loci did not show locus-specific clusters of polymorphic alleles. This was considered to be because of either interlocus genetic exchange or the confusing designation of loci in different haplotypes at the present time. The locus specificity of polymorphic alleles was also observed in human and mouse MHC class II loci. It was therefore concluded that interlocus recombination or gene conversion is not very important for generating MHC diversity, with a possible exception of mouse class I loci. According to the phylogenetic trees of complete coding sequences, we classified human MHC class I (HLA-A, -B, and -C) and class II (DRB1) alleles into three to five major allelic lineages (groups), which were monophyletic with high bootstrap values. Most of these allelic groups remained unchanged even in phylogenetic trees based on individual exons, though this does not exclude the possibility of intralocus recombination involving short DNA segments. These results, together with the previous observation that MHC loci are subject to frequent duplication and deletion, as well as to balancing selection, indicate that MHC evolution in mammals is in agreement with the birth-and-death model of evolution, rather than with the model of concerted evolution.     

Allelic variation at the <Emphasis Type="Italic">VRN-1</Emphasis> promoter region in polyploid wheat

Vernalization, the requirement of a long exposure to low temperatures to induce flowering, is an essential adaptation of plants to cold winters. We have shown recently that the vernalization gene VRN-1 from diploid wheat Triticum monococcum is the meristem identity gene APETALA1, and that deletions in its promoter were associated with spring growth habit. In this study, we characterized the allelic variation at the VRN-1 promoter region in polyploid wheat. The Vrn-A1a allele has a duplication including the promoter region. Each copy has similar foldback elements inserted at the same location and is flanked by identical host direct duplications (HDD). This allele was found in more than half of the hexaploid varieties but not among the tetraploid lines analyzed here. The Vrn-A1b allele has two mutations in the HDD region and a 20-bp deletion in the 5 UTR compared with the winter allele. The Vrn-A1b allele was found in both tetraploid and hexaploid accessions but at a relatively low frequency. Among the tetraploid wheat accessions, we found two additional alleles with 32 bp and 54 bp deletions that included the HDD region. We found no size polymorphisms in the promoter region among the winter wheat varieties. The dominant Vrn-A1 allele from two spring varieties from Afghanistan and Egypt (Vrn-A1c allele) and all the dominant Vrn-B1 and Vrn-D1 alleles included in this study showed no differences from their respective recessive alleles in promoter sequences. Based on these results, we concluded that the VRN-1 genes should have additional regulatory sites outside the promoter region studied here.     

Extensive sharing of MHC class II alleles between rhesus and cynomolgus macaques

In contrast to rhesus monkeys, substantial knowledge on cynomolgus monkey major histocompatibility complex (MHC) class II haplotypes is lacking. Therefore, 17 animals, including one pedigreed family, were thoroughly characterized for polymorphic Mhc class II region genes as well as their mitochondrial DNA (mtDNA) sequences. Different cynomolgus macaque populations appear to exhibit unique mtDNA profiles reflecting their geographic origin. Within the present panel, 10 Mafa-DPB1, 14 Mafa-DQA1, 12 Mafa-DQB1, and 35 Mafa-DRB exon 2 sequences were identified. All of these alleles cluster into lineages that were previously described for rhesus macaques. Moreover, about half of the Mafa-DPB1, Mafa-DQA1, and Mafa-DQB1 alleles and one third of the Mafa-DRB exon 2 sequences are identical to rhesus macaque orthologues. Such a high level of Mhc class II allele sharing has not been reported for primate species. Pedigree analysis allowed the characterization of nine distinct Mafa class II haplotypes, and seven additional ones could be deduced. Two of these haplotypes harbor a duplication of the Mafa-DQB1 locus. Despite extensive allele sharing, rhesus and cynomolgus monkeys do not appear to possess identical Mhc class II haplotypes, thus illustrating that new haplotypes were generated after speciation by recombination-like processes.     

Genetic markers of adaptive processes in the Far Eastern pink salmon <Emphasis Type="Italic">Oncorhynchus gorbuscha</Emphasis>: Allelic diversity at the locus of major histocompatibility complex MHC I-<Emphasis Type="Italic">A1</Emphasis>

To clarify allelic diversity at the locus of major histocompatibility complex MHC class I-A1 in the Far Eastern pink salmon Oncorhynchus gorbuscha, sequencing of the electrophoretic alleles isolated from the gel (DGGE alleles) was performed. In 47 individuals, the genotypes of which consisted of ten DGGE alleles, 18 MHC I-A1 nucleotide sequences were revealed, and thus, eight cryptic alleles not detected by electrophoresis were identified. Eleven of these alleles were identified earlier in pink salmon from Hokkaido, Alaska, and British Columbia, and seven, possibly, were unique to the populations from some Far Eastern regions. Six of the previously determined DGGE alleles corresponded to more than one nucleotide sequence. However, the sequences attributed to the same DGGE allele differed on average by less than 1 nucleotide. These findings point to sufficient sensitivity of the DGGE method, although the genetic diversity and differentiation estimates obtained with it will obviously be somewhat underestimated. Considerable predominance of nonsynonymous substitutions over the synonymous ones in the codons of the MHC I-A1 antigen-binding site confirms the presence of positive selection aimed at providing the population resistance to local spectrum of pathogens. Refinement of the allelic composition of the adaptively important MHC genetic marker will contribute to more complete understanding of the adaptive genetic structure of pink salmon as an important element of the overall population structure of the species.     

Relic of ancient recombinations in gibbon ABO blood group genes deciphered through phylogenetic network analysis

The primate ABO blood group gene encodes a glycosyl transferase (either A or B type), and is known to have large coalescence times among the allelic lineages in human. We determined nucleotide sequences of ca. 2.2 kb of this gene for 23 individuals of three gibbon species (agile gibbon, white-handed gibbon, and siamang), and observed a total of 24 haplotypes. We found relics of five ancient intragenic recombinations, occurred during ca. 2–7 million years ago, through a phylogenetic network analysis. The coalescence time between A and B alleles estimate precede the divergence (ca. 8 MYA) of siamang and common gibbon lineages. This establishes the coexistence of divergent allelic lineages of the ABO blood group gene for a long period in the ancestral gibbon species, and strengthens the non-neutral evolution for this gene.     

MHC-linked class III genes. Analysis of C4 gene frequencies,complotypes and associations with distinct HLA haplotypes in German Caucasians

The class III complement components, C4, C2 and factor B (BF), are encoded in the human major histocompatibility complex (MHC). The two genes determining C4 (C4A and C4B) display considerable polymorphism and, thus, are important markers for HLA. In combination with alleles of C2 and BF they can be grouped into unique complotypes. We have analyzed the C4 alleles in a panel of 204 unrelated German Caucasians and studied their segregation with HLA haplotypes in 24 normal families. Inclusion of the class III markers with the class I and 11 alleles provides a more refined picture of the genetic structure of the MHC in these families. When charted according to the HLA-B locus specificities the MHCs can be clustered into groups showing distinctly homogenous or heterogenous complotypes. The identification of such groups is valuable for the selection of genetic material to analyze the molecular genetics of the human MHC.Abbreviations BF factor B - C2 second component of complement - C4 fourth component of complement - EDTA ethylenediamine tetraacetate - GLO glyoxalase-I - MHC major histocompatibility complex     

Two ancient allelic lineages at the single classical class I locus in the Xenopus MHC

Unlike all other vertebrates examined to date, there is only one detectable class I locus in the Xenopus MHC. On the bases of a nearly ubiquitous and high tissue expression, extensive polymorphism, and MHC linkage, this gene is of the classical or class Ia type. Sequencing analysis of class Ia cDNAs encoded by eight defined MHC haplotypes reveals two very old allelic lineages that perhaps emerged when humans and mice diverged from a common ancestor up to 100 million years ago. The unprecedented age of these lineages suggests that different class Ia genes from ancestors of the laboratory model Xenopus laevis are now expressed as alleles in this species. The lineages are best defined by their cytoplasmic and alpha2 peptide-binding domains, and there are highly diverse alleles (defined by the alpha1 peptide-binding domain) in each lineage. Surprisingly, the alpha3 domains are homogenized in both lineages, suggesting that interallelic gene conversion/recombination maintains the high sequence similarity.     

Allelic Combinations of Promoter and Exon 2 in <Emphasis Type="Italic">DQB1</Emphasis> in Dogs and Wolves

Polymorphism of PBRs of the major histocompatibility complex (MHC) genes is well recognized, but the polymorphism also extends to proximal promoter regions. Examining DQB1 variability in dogs and wolves, we identified 7 promoter variants and 13 exon 2 alleles among 89 dogs, including a previously unknown DQB1 exon 2 allele, and 8 promoter variants and 9 exon 2 alleles among 85 wolves. As expected from previous studies and from a close chromosomal location, strong linkage disequilibrium was demonstrated in both wolves and dogs by having significantly fewer promoter/exon 2 combinations than expected from simulations of randomized data sets. Interestingly, we noticed weaker haplotypic associations in dogs than in wolves. Dogs had twice as many promoter/exon 2 combinations as wolves and an almost 2-fold difference in the number of exon 2 alleles per promoter variant. This difference was not caused by an admixture of breeds in our group of dogs because the high ratio of observed to expected number of haplotypes persisted within a single dog breed, the German Shepherd. Ewens-Watterson tests indicated that both the promoter and exon 2 are under the balancing selection, and both regions appear to be more recently derived in the dog than in the wolf. Hence, although reasons for the differences are unknown, they may relate to altered selection pressure on patterns of expression. Deviations from normal MHC expression patterns have been associated with autoimmune diseases, which occur frequently in several dog breeds. Further knowledge about these deviations may help us understand the source of such diseases.     

Identification of MhcMafa-DRB Alleles in a Cohort of Cynomolgus Macaques of Vietnamese Origin

Cynomolgus macaques have been used widely to build a research model of infectious and chronic diseases, as well as in transplantation studies, where disease susceptibility and/or resistance are associated with the major histocompatibility complex (MHC). To better elucidate polymorphisms and genetic differences in the Mafa‐DRB locus, and facilitate the experimental use of cynomolgus macaques, we used pool screening combined with cloning and direct sequencing of polymerase chain reaction products to characterize MhcMafa‐DRB gene alleles in 153 Vietnamese cynomolgus macaques. We identified 30 Mafa‐DRB alleles belonging to 17 allelic lineages, including four novel sequences that had not been documented in earlier reports. The highest frequency allele was Mafa‐DRB*W27:04, which was present in 7 of 35 (20%) monkeys. The next most frequent alleles were Mafa‐DRB*3:07 and Mafa‐DRB*W7:01, which were detected in 5 of 35 (14.3%) and 4 of 35 (11.4%) of the monkeys, respectively. The high‐frequency alleles in this Vietnamese population may be high priority targets for additional characterization of immune functions. Only the DRB1*03 and DRB1*10 lineages were also present in humans, whereas the remaining alleles were monkey‐specific lineages. We found 25 variable sites by aligning the deduced amino acid sequences of 29 identified alleles. Evolutionary and population analyses based on these sequences showed that human, rhesus, and cynomolgus macaques share several Mhc‐DRB lineages and the shared polymorphisms in the DRB region may be attributable to the existence of interbreeding between rhesus and cynomolgus macaques. This information will promote the understanding of MHC diversity and polymorphism in cynomolgus macaques and increase the value of this species as a model for biomedical research. Am. J. Primatol. 74:958‐966, 2012. © 2012 Wiley Periodicals, Inc.     

High-resolution typing for chicken BF2 (MHC class I) alleles by automated sequencing

Sequence-based typing (SBT) was developed for major histocompatibility complex (MHC) class I and class II alleles in humans. We report here the development and application of a SBT method for alleles of the chicken BF2 locus (the more polymorphic of the two MHC class I loci in chickens). Exon 2 of the BF2 gene was selectively amplified from genomic DNA using a BF2 locus-specific PCR primer. Exon 2 sequences were sufficient to identify the 21 distinct BF2 alleles described in standard B haplotypes of Leghorns and in commercial broiler-breeder lines. Sixty-six samples from MHC typed, pedigreed chickens were tested, including 50 different heterozygous combinations. BF2 sequences from all B homozygotes were successfully amplified, and all combinations of BF2 alleles in heterozygotes were co-amplified equally. The two different BF2 alleles in heterozygotes could be identified unambiguously by distinct sequence motif patterns. In tests of samples of unknown B genotype in commercial broiler-breeder flocks, we identified expected BF2 alleles as well as an allele not previously encountered in one of the lines.