Evolution of a new nonclassical MHC class I locus in two Old World primate species

 HLA-G is a nonclassical major histocompatibility complex (MHC) class I molecule that is expressed only in the human placenta, suggesting that it plays an important role at the fetal-maternal interface. In rhesus monkeys, which have similar placentation to humans, the HLA-G orthologue is a pseudogene. However, rhesus monkeys express a novel placental MHC class I molecule, Mamu-AG, which has HLA-G-like characteristics. Phylogenetic analysis of AG alleles in two Old World primate species, the baboon and the rhesus macaque, revealed limited diversity characteristic of a nonclassical MHC class I locus. Gene trees constructed using classical and nonclassical primate MHC class I alleles demonstrated that the AG locus was most closely related to the classical A locus. Interestingly, gene tree analyses suggested that the AG alleles were most closely related to a subset of A alleles which are the products of an ancestral interlocus recombination event between the A and B loci. Calculation of the rates of synonymous and nonsynonymous substitution at the AG locus revealed that positive selection was not acting on the codons encoding the peptide binding region. In exon 4, however, the rate of nonsynonymous substitution was significantly lower than the rate of synonymous substitution, suggesting that negative selection was acting on these codons. Received: 22 April 1998 / Revised: 15 July 1998     

Diversity of MHC class I haplotypes in cynomolgus macaques

Cynomolgus macaques are widely used as a primate model for human diseases associated with an immunological process. Because there are individual differences in immune responsiveness, which are controlled by the polymorphic nature of the major histocompatibility (MHC) locus, it is important to reveal the diversity of MHC in the model animal. In this study, we analyzed 26 cynomolgus macaques from five families for MHC class I genes. We identified 32 Mafa-A, 46 Mafa-B, 6 Mafa-I, and 3 Mafa-AG alleles in which 14, 20, 3, and 3 alleles were novel. There were 23 MHC class I haplotypes and each haplotype was composed of one to three Mafa-A alleles and one to five Mafa-B alleles. Family studies revealed that there were two haplotypes which contained two Mafa-A1 alleles. These observations demonstrated further the complexity of MHC class I locus in the Old World monkey.     

Genome-wide association study identified the human leukocyte antigen region as a novel locus for plasma beta-2 microglobulin

Beta-2 microglobulin (B2M) is a component of the major histocompatibility complex (MHC) class I molecule and has been studied as a biomarker of kidney function, cardiovascular diseases and mortality. Little is known about the genes influencing its levels directly or through glomerular filtration rate (GFR). We conducted a genome-wide association study of plasma B2M levels in 6738 European Americans from the Atherosclerosis Risk in Communities study to identify novel loci for B2M and assessed its association with known estimated GFR (eGFR) loci. We identified 2 genome-wide significant loci. One was in the human leukocyte antigen (HLA) region on chromosome 6 (lowest p value = 1.8 × 10^?23 for rs9264638). At this locus, 6 index SNPs accounted for 3.2 % of log(B2M) variance, and their association with B2M could largely be explained by imputed classical alleles of the MHC class I genes: HLA-A, HLA-B, or HLA-C. The index SNPs at this locus were not associated with eGFR based on serum creatinine (eGFRcr). The other locus of B2M was on chromosome 12 (rs3184504 at SH2B3, beta = 0.02, p value = 3.1 × 10^?8), which was previously implicated as an eGFR locus. In conclusion, although B2M is known to be a component of MHC class I molecule, the association between HLA class I alleles and plasma B2M levels in a community-based population is novel. The identification of the two novel loci for B2M extends our understanding of its metabolism and informs its use as a kidney filtration biomarker.     

MHC class I genes in a New World primate, the cotton-top tamarin (Saguinus oedipus), have evolved by an active process of loci turnover

 Lymphocytes of a New World primate, the cotton-top tamarin (Saguinus oedipus), express classical G–related major histocompatibility complex (MHC) class I molecules with unusually limited polymorphism and variability. Three G-related loci, an F locus, an E locus, and two pseudogenes (So-N1 and So-N3) have been identified by cDNA library screening and extensive PCR analysis of both cDNA and genomic DNA from the cotton-top tamarin. Furthermore, each genus of the subfamily Callitrichinae (tamarins and marmosets) appears to express its own unique set of MHC class I genes, likely due to a rapid turnover of loci. The rapid emergence of unique MHC class I genes in the Callitrichinae genera, resulting from an active process of duplication and inactivation of loci, may account for the limited diversity of the MHC class I genes in the cotton-top tamarin. To determine the nature of the entire complement of MHC class I genes in the cotton-top tamarin, we synthesized a genomic DNA library and screened it with MHC class I-specific probes. We isolated nine new MHC class I pseudogenes from this library. These newly isolated tamarin G–related MHC class I pseudogenes are not closely related to any of their functional counterparts in the tamarin, suggesting that they do not share a recent common ancestral gene with the tamarin's currently expressed MHC class I loci. In addition, these tamarin sequences display a high rate of nonsynonymous substitutions in their putative peptide binding region. This indicates that the genes from which they have derived were likely subject to positive selection and, therefore, were once functional. Our data support the notion that an extremely high rate of loci turnover is largely responsible for the limited diversity of the MHC class I genes in the cotton-top tamarin. Received: 15 September 1997 / Revised: 2 July 1998     

Construction and Application of a Korean Reference Panel for Imputing Classical Alleles and Amino Acids of Human Leukocyte Antigen Genes

Genetic variations of human leukocyte antigen (HLA) genes within the major histocompatibility complex (MHC) locus are strongly associated with disease susceptibility and prognosis for many diseases, including many autoimmune diseases. In this study, we developed a Korean HLA reference panel for imputing classical alleles and amino acid residues of several HLA genes. An HLA reference panel has potential for use in identifying and fine-mapping disease associations with the MHC locus in East Asian populations, including Koreans. A total of 413 unrelated Korean subjects were analyzed for single nucleotide polymorphisms (SNPs) at the MHC locus and six HLA genes, including HLA-A, -B, -C, -DRB1, -DPB1, and -DQB1. The HLA reference panel was constructed by phasing the 5,858 MHC SNPs, 233 classical HLA alleles, and 1,387 amino acid residue markers from 1,025 amino acid positions as binary variables. The imputation accuracy of the HLA reference panel was assessed by measuring concordance rates between imputed and genotyped alleles of the HLA genes from a subset of the study subjects and East Asian HapMap individuals. Average concordance rates were 95.6% and 91.1% at 2-digit and 4-digit allele resolutions, respectively. The imputation accuracy was minimally affected by SNP density of a test dataset for imputation. In conclusion, the Korean HLA reference panel we developed was highly suitable for imputing HLA alleles and amino acids from MHC SNPs in East Asians, including Koreans.     

Nonequivalence of Classical MHC Class I Loci in Ability to Direct Effective Antiviral Immunity

Structural diversity in the peptide binding sites of the redundant classical MHC antigen presenting molecules is strongly selected in humans and mice. Although the encoded antigen presenting molecules overlap in antigen presenting function, differences in polymorphism at the MHC I A, B and C loci in humans and higher primates indicate these loci are not functionally equivalent. The structural basis of these differences is not known. We hypothesize that classical class I loci differ in their ability to direct effective immunity against intracellular pathogens. Using a picornavirus infection model and chimeric H-2 transgenes, we examined locus specific functional determinants distinguishing the ability of class I sister genes to direct effective anti viral immunity. Whereas, parental FVB and transgenic FVB mice expressing the H-2K^b gene are highly susceptible to persisting Theiler''s virus infection within the CNS and subsequent demyelination, mice expressing the D^b transgene clear the virus and are protected from demyelination. Remarkably, animals expressing a chimeric transgene, comprised primarily of K^b but encoding the peptide binding domain of D^b, develop a robust anti viral CTL response yet fail to clear virus and develop significant demyelination. Differences in expression of the chimeric K^bα1α2D^b gene (low) and D^b (high) in the CNS of infected mice mirror expression levels of their endogenous H-2^q counterparts in FVB mice. These findings demonstrate that locus specific elements other than those specifying peptide binding and T cell receptor interaction can determine ability to clear virus infection. This finding provides a basis for understanding locus-specific differences in MHC polymorphism, characterized best in human populations.     

The MHC E locus in macaques is polymorphic and is conserved between macaques and humans

Although the functions of the molecules encoded by the classical MHC class I loci are well defined, no function has been ascribed to the molecules encoded by the non-classical MHC class I loci. To investigate the evolution and conservation of the non-classical loci, we cloned and sequenced HLA-E homologues in macaques. We isolated four E locus alleles from five rhesus monkeys and two E locus alleles from one cynomolgus monkey, which indicated that the E locus in macaques is polymorphic. We also compared the rate of nucleotide substitution in the second intron of the macaque and human E locus alleles with that of exons two and three. The rate of nucleotide substitution was significantly higher in the introns, which suggested that the E locus has evolved under selective pressure. Additionally, comparison of the rates of synonymous and non-synonymous substitutions in the peptide binding region versus the remainder of the molecule suggested that the codons encoding the amino acids in the peptide binding region had been conserved in macaques and humans over the 36 million years since macaques and humans last shared a common ancestor.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers UO2976–UO2981     

Codon usage bias and base composition in MHC genes in humans and common chimpanzees

 Codon bias and base composition in major histocompatibility complex (MHC) sequences have been studied for both class I and II loci in Homo sapiens and Pan troglodytes. There is low to moderate codon bias for the MHC of humans and chimpanzees. In the class I loci, the same level of moderate codon bias is seen for HLA-B, HLA-C, Patr-A, Patr-B, and Patr-C, while at HLA-A the level of codon bias is lower. There is a correlation between codon usage bias and G+C content in the A and B loci in humans and chimps, but not at the C locus. To examine the effect of diversifying selection on codon bias, we subdivided class I alleles into antigen recognition site (ARS) and non-ARS codons. ARS codons had lower bias than non-ARS codons. This may indicate that the constraint of codon bias on nucleotide substitution may be selected against in ARS codons. At the class II loci, there are distinct differences between alpha and beta chain genes with respect to codon usage, with the beta chain genes being much more biased. Species-specific differences in base composition were seen in exon 2 at the DRB1 locus, with lower GC content in chimpanzees. Considering the complex evolutionary history of MHC genes, the study of codon usage patterns provides us with a better understanding of both the evolutionary history of these genes and the evolution of synonymous codon usage in genes under natural selection. Received: 2 April 1998 / Revised: 2 September 1998     

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     

The dominant MHC class I gene is adjacent to the polymorphic<Emphasis Type="Italic"> TAP2</Emphasis> gene in the duck,<Emphasis Type="Italic"> Anas platyrhynchos</Emphasis>

We are investigating the expression and linkage of major histocompatibility complex (MHC) class I genes in the duck (Anas platyrhynchos) with a view toward understanding the susceptibility of ducks to two medically important viruses: influenza A and hepatitis B. In mammals, there are multiple MHC class I loci, and alleles at a locus are polymorphic and co-dominantly expressed. In contrast, in lower vertebrates the expression of one locus predominates. Southern-blot analysis and amplification of genomic sequences suggested that ducks have at least four loci encoding MHC class I. To identify expressed MHC genes, we constructed an unamplified cDNA library from the spleen of a single duck and screened for MHC class I. We sequenced 44 positive clones and identified four MHC class I sequences, each sharing approximately 85% nucleotide identity. Allele-specific oligonucleotide hybridization to a Northern blot indicated that only two of these sequences were abundantly expressed. In chickens, the dominantly expressed MHC class I gene lies adjacent to the transporter of antigen processing (TAP2) gene. To investigate whether this organization is also found in ducks, we cloned the gene encoding TAP2 from the cDNA library. PCR amplification from genomic DNA allowed us to determine that the dominantly expressed MHC class I gene was adjacent to TAP2. Furthermore, we amplified two alleles of the TAP2 gene from this duck that have significant and clustered amino acid differences that may influence the peptides transported. This organization has implications for the ability of ducks to eliminate viral pathogens.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers AY294416–22     

Fine-Mapping the Genetic Association of the Major Histocompatibility Complex in Multiple Sclerosis: HLA and Non-HLA Effects

The major histocompatibility complex (MHC) region is strongly associated with multiple sclerosis (MS) susceptibility. HLA-DRB1*15:01 has the strongest effect, and several other alleles have been reported at different levels of validation. Using SNP data from genome-wide studies, we imputed and tested classical alleles and amino acid polymorphisms in 8 classical human leukocyte antigen (HLA) genes in 5,091 cases and 9,595 controls. We identified 11 statistically independent effects overall: 6 HLA-DRB1 and one DPB1 alleles in class II, one HLA-A and two B alleles in class I, and one signal in a region spanning from MICB to LST1. This genomic segment does not contain any HLA class I or II genes and provides robust evidence for the involvement of a non-HLA risk allele within the MHC. Interestingly, this region contains the TNF gene, the cognate ligand of the well-validated TNFRSF1A MS susceptibility gene. The classical HLA effects can be explained to some extent by polymorphic amino acid positions in the peptide-binding grooves. This study dissects the independent effects in the MHC, a critical region for MS susceptibility that harbors multiple risk alleles.     

MHC class I genes in the owl monkey: mosaic organisation, convergence and loci diversity

The MHC class I molecule plays an important role in immune response, pathogen recognition and response against vaccines and self- versus non-self-recognition. Studying MHC class I characteristics thus became a priority when dealing with Aotus to ensure its use as an animal model for biomedical research. Isolation, cloning and sequencing of exons 1–8 from 27 MHC class I alleles obtained from 13 individuals classified as belonging to three owl monkey species (A. nancymaae, A. nigriceps and A. vociferans) were carried out to establish similarities between Aotus MHC class I genes and those expressed by other New and Old World primates. Six Aotus MHC class I sequence groups (Ao-g1, Ao-g2, Ao-g3, Ao-g4, Ao-g5 and Ao-g6) weakly related to non-classical Catarrhini MHC were identified. An allelic lineage was also identified in one A. nancymaae and two A. vociferans monkeys, exhibiting a high degree of conservation, negative selection along the molecule and premature termination of the open reading frame at exon 5 (Ao-g5). These sequences high conservation suggests that they more likely correspond to a soluble form of Aotus MHC class I molecules than to a new group of processed pseudogenes. Another group, named Ao-g6, exhibited a strong relationship with Catarrhinis classical MHC-B-C loci. Sequence evolution and variability analysis indicated that Aotus MHC class I molecules experience inter-locus gene conversion phenomena, contributing towards their high variability.     

Single locus typing of MHC class I and class II B loci in a population of red jungle fowl

In species with duplicated major histocompatibility complex (MHC) genes, estimates of genetic variation often rely on multilocus measures of diversity. It is possible that such measures might not always detect more detailed patterns of selection at individual loci. Here, we describe a method that allows us to investigate classical MHC diversity in red jungle fowl (Gallus gallus), the wild ancestor of the domestic chicken, using a single locus approach. This is possible due to the well-characterised gene organisation of the ‘minimal essential’ MHC (BF/BL region) of the domestic chicken, which comprises two differentially expressed duplicated class I (BF) and two class II B (BLB) genes. Using a combination of reference strand-mediated conformation analysis, cloning and sequencing, we identify nine BF and ten BLB alleles in a captive population of jungle fowl. We show that six BF and five BLB alleles are from the more highly expressed locus of each gene, BF2 and BLB2, respectively. An excess of non-synonymous substitutions across the jungle fowl BF/BL region suggests that diversifying selection has acted on this population. Importantly, single locus screening reveals that the strength of selection is greatest on the highly expressed BF2 locus. This is the first time that a population of red jungle fowl has been typed at the MHC region, laying the basis for further research into the underlying processes acting to maintain MHC diversity in this and other species.     

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

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

A Mechanistic Basis for the Co-evolution of Chicken Tapasin and Major Histocompatibility Complex Class I (MHC I) Proteins

MHC class I molecules display peptides at the cell surface to cytotoxic T cells. The co-factor tapasin functions to ensure that MHC I becomes loaded with high affinity peptides. In most mammals, the tapasin gene appears to have little sequence diversity and few alleles and is located distal to several classical MHC I loci, so tapasin appears to function in a universal way to assist MHC I peptide loading. In contrast, the chicken tapasin gene is tightly linked to the single dominantly expressed MHC I locus and is highly polymorphic and moderately diverse in sequence. Therefore, tapasin-assisted loading of MHC I in chickens may occur in a haplotype-specific way, via the co-evolution of chicken tapasin and MHC I. Here we demonstrate a mechanistic basis for this co-evolution, revealing differences in the ability of two chicken MHC I alleles to bind and release peptides in the presence or absence of tapasin, where, as in mammals, efficient self-loading is negatively correlated with tapasin-assisted loading. We found that a polymorphic residue in the MHC I α3 domain thought to bind tapasin influenced both tapasin function and intrinsic peptide binding properties. Differences were also evident between the MHC alleles in their interactions with tapasin. Last, we show that a mismatched combination of tapasin and MHC alleles exhibit significantly impaired MHC I maturation in vivo and that polymorphic MHC residues thought to contact tapasin influence maturation efficiency. Collectively, this supports the possibility that tapasin and BF2 proteins have co-evolved, resulting in allele-specific peptide loading in vivo.     

Evidence for Directional Selection at a Novel Major Histocompatibility Class I Marker in Wild Common Frogs (Rana temporaria) Exposed to a Viral Pathogen (Ranavirus)

Whilst the Major Histocompatibility Complex (MHC) is well characterized in the anuran Xenopus, this region has not previously been studied in another popular model species, the common frog (Rana temporaria). Nor, to date, have there been any studies of MHC in wild amphibian host-pathogen systems. We characterise an MHC class I locus in the common frog, and present primers to amplify both the whole region, and specifically the antigen binding region. As no more than two expressed haplotypes were found in over 400 clones from 66 individuals, it is likely that there is a single class I locus in this species. This finding is consistent with the single class I locus in Xenopus, but contrasts with the multiple loci identified in axolotls, providing evidence that the diversification of MHC class I into multiple loci likely occurred after the Caudata/Anura divergence (approximately 350 million years ago) but before the Ranidae/Pipidae divergence (approximately 230 mya). We use this locus to compare wild populations of common frogs that have been infected with a viral pathogen (Ranavirus) with those that have no history of infection. We demonstrate that certain MHC supertypes are associated with infection status (even after accounting for shared ancestry), and that the diseased populations have more similar supertype frequencies (lower F_ST) than the uninfected. These patterns were not seen in a suite of putatively neutral microsatellite loci. We interpret this pattern at the MHC locus to indicate that the disease has imposed selection for particular haplotypes, and hence that common frogs may be adapting to the presence of Ranavirus, which currently kills tens of thousands of amphibians in the UK each year.     

Selection, trans-species polymorphism, and locus identification of major histocompatibility complex class IIβ alleles of New World ranid frogs

Genes encoded by the major histocompatibility complex (MHC) play key roles in the vertebrate immune system. However, our understanding of the evolutionary processes and underlying genetic mechanisms shaping these genes is limited in many taxa, including amphibians, a group currently impacted by emerging infectious diseases. To further elucidate the evolution of the MHC in frogs (anurans) and develop tools for population genetics, we surveyed allelic diversity of the MHC class II β1 domain in both genomic and complementary DNA of seven New World species in the genus Rana (Lithobates). To assign locus affiliation to our alleles, we used a “gene walking” technique to obtain intron 2 sequences that flanked MHC class IIβ exon 2. Two distinct intron sequences were recovered, suggesting the presence of at least two class IIβ loci in Rana. We designed a primer pair that successfully amplified an orthologous locus from all seven Rana species. In total, we recovered 13 alleles and documented trans-species polymorphism for four of the alleles. We also found quantitative evidence of selection acting on amino acid residues that are putatively involved in peptide binding and structural stability of the β1 domain of anurans. Our results indicated that primer mismatch can result in polymerase chain reaction (PCR) bias, which influences the number of alleles that are recovered. Using a single locus may minimize PCR bias caused by primer mismatch, and the gene walking technique was an effective approach for generating single-copy orthologous markers necessary for future studies of MHC allelic variation in natural amphibian populations.     

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.