A hypothetical model of the foreign antigen binding site of class II histocompatibility molecules

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle)

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

MHC-B Diversity and Signs of Respiratory Illness in Wild, East African Chimpanzees ( Pan troglodytes schweinfurthii )

Many traits, intrinsic and extrinsic to an organism, contribute to interindividual variation in immunity in wild habitats. The vertebrate Major Histocompatibility Complex (MHC) includes genes encoding antigen-presenting molecules that are highly variable, and that variation often predicts susceptibility/resistance to and recovery from pathogen infection. I compare MHC-B variation at two long-term chimpanzee research sites, Kibale National Park in Uganda and Gombe National Park in Tanzania. Using decades of respiratory health data available for these chimpanzees, I test hypotheses associated with maintenance of diversity at MHC loci, including heterozygote, divergent allele, and rare allele advantage hypotheses, and predictions for unique function of MHC-B in great apes. I found, despite confirmation of recent shared ancestry between Kibale and Gombe chimpanzees, including an overlapping MHC-B allele repertoire and similar MHC-B phenotype compositions, chimpanzees from the two research sites experienced differences in the occurrence of respiratory signs and had different associations of MHC-B diversity with signs of respiratory illness. Kibale chimpanzees with heterozygous genotypes and different peptide-binding supertypes were observed less often with respiratory signs than those homozygous or possessing the same supertypes, but this same association was not observed among Gombe chimpanzees. Gombe chimpanzees with specific MHC-B phenotypes that enable engagement of Natural Killer (NK) cells were observed more often with respiratory signs than chimpanzees with other phenotypes, but this was not observed at Kanyawara. This study emphasizes local adaptation in shaping genetic and phenotypic traits in different infectious disease contexts, even among close genetic relatives of the same subspecies, and highlights utility for continued and simultaneous tracking of host immune genes and specific pathogens in wild species.

Juvenile idiopathic arthritis-associated genetic loci exhibit spatially constrained gene regulatory effects across multiple tissues and immune cell types

Juvenile idiopathic arthritis (JIA) is an autoimmune, inflammatory joint disease with complex genetic etiology. Previous GWAS have found many genetic loci associated with JIA. However, the biological mechanism behind JIA remains unknown mainly because most risk loci are located in non-coding genetic regions. Interestingly, increasing evidence has found that regulatory elements in the non-coding regions can regulate the expression of distant target genes through spatial (physical) interactions. Here, we used information on the 3D genome organization (Hi-C data) to identify target genes that physically interact with SNPs within JIA risk loci. Subsequent analysis of these SNP-gene pairs using data from tissue and immune cell type-specific expression quantitative trait loci (eQTL) databases allowed the identification of risk loci that regulate the expression of their target genes. In total, we identified 59 JIA-risk loci that regulate the expression of 210 target genes across diverse tissues and immune cell types. Functional annotation of spatial eQTLs within JIA risk loci identified significant overlap with gene regulatory elements (i.e., enhancers and transcription factor binding sites). We found target genes involved in immune-related pathways such as antigen processing and presentation (e.g., ERAP2, HLA class I and II), the release of pro-inflammatory cytokines (e.g., LTBR, TYK2), proliferation and differentiation of specific immune cell types (e.g., AURKA in Th17 cells), and genes involved in physiological mechanisms related to pathological joint inflammation (e.g., LRG1 in arteries). Notably, many of the tissues where JIA-risk loci act as spatial eQTLs are not classically considered central to JIA pathology. Overall, our findings highlight the potential tissue and immune cell type-specific regulatory changes contributing to JIA pathogenesis. Future integration of our data with clinical studies can contribute to the development of improved JIA therapy.

Incorporating the Molecular Mimicry of Environmental Antigens into the Causality of Autoimmune Hepatitis

Molecular mimicry between foreign and self-antigens has been implicated as a cause of autoimmune hepatitis in experimental models and cross-reacting antibodies in patients. This review describes the experimental and clinical evidence for molecular mimicry as a cause of autoimmune hepatitis, indicates the limitations and uncertainties of this premise, and encourages investigations that assess diverse environmental antigens as sources of disease-relevant molecular mimics. Pertinent articles were identified in PubMed using multiple search phrases. Several pathogens have linear or conformational epitopes that mimic the self-antigens of autoimmune hepatitis. The occurrence of an acute immune-mediated hepatitis after vaccination for severe acute respiratory syndrome (SARS)-associated coronavirus 2 (SARS-CoV-2) has suggested that vaccine-induced peptides may mimic disease-relevant tissue antigens. The intestinal microbiome is an under-evaluated source of gut-derived antigens that could also engage in molecular mimicry. Chaperone molecules may enhance the pathogenicity of molecular mimics, and they warrant investigation. Molecular mimics of immune dominant epitopes within cytochrome P450 IID6, the autoantigen most closely associated with autoimmune hepatitis, should be sought in diverse environmental antigens and assessed for pathogenicity. Avoidance strategies, dietary adjustments, vaccine improvement, and targeted manipulation of the intestinal microbiota may emerge as therapeutic possibilities. In conclusion, molecular mimicry may be a missing causality of autoimmune hepatitis. Molecular mimics of key immune dominant epitopes of disease-specific antigens must be sought in diverse environmental antigens. The ubiquity of molecular mimicry compels rigorous assessments of peptide mimics for immunogenicity and pathogenicity in experimental models. Molecular mimicry may complement epigenetic modifications as causative mechanisms of autoimmune hepatitis.

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

Relationship between genome-wide and MHC class I and II genetic diversity and complementarity in a nonhuman primate

Although mate choice is expected to favor partners with advantageous genetic properties, the relative importance of genome-wide characteristics, such as overall heterozygosity or kinship, versus specific loci, is unknown. To disentangle genome-wide and locus-specific targets of mate choice, we must first understand congruence in global and local variation within the same individual. This study compares genetic diversity, both absolute and relative to other individuals (i.e., complementarity), assessed across the genome to that found at the major histocompatibility complex (MHC), a hyper-variable gene family integral to immune system function and implicated in mate choice across species. Using DNA from 22 captive olive baboons (Papio anubis), we conducted double digest restriction site-associated DNA sequencing to estimate genome-wide heterozygosity and kinship, and sequenced two class I and two class II MHC loci. We found that genome-wide diversity was not associated with MHC diversity, and that diversity at class I MHC loci was not correlated with diversity at class II loci. Additionally, kinship was a significant predictor of the number of MHC alleles shared between dyads at class II loci. Our results provide further evidence of the strong selective pressures maintaining genetic diversity at the MHC in comparison to other randomly selected sites throughout the genome. Furthermore, our results indicate that class II MHC disassortative mate choice may mediate inbreeding avoidance in this population. Our study suggests that mate choice favoring genome-wide genetic diversity is not always synonymous with mate choice favoring MHC diversity, and highlights the importance of controlling for kinship when investigating MHC-associated mate choice.

Using de novo genome assembly and high-throughput sequencing to characterize the MHC region in a non-model bird, the Eurasian coot

Genes of the Major Histocompatibility Complex (MHC) form a key component of vertebrate adaptive immunity, as they code for molecules which bind antigens of intra- and extracellular pathogens (MHC class I and II, respectively) and present them to T cell receptors. In general, MHC genes are hyper-polymorphic and high MHC diversity is often maintained within natural populations (via balancing selection) and within individuals (via gene duplications). Because of its complex architecture with tandems of duplicated genes, characterization of MHC region in non-model vertebrate species still poses a major challenge. Here, we combined de novo genome assembly and high-throughput sequencing to characterize MHC polymorphism in a rallid bird species, the Eurasian coot Fulica atra. An analysis of genome assembly indicated high duplication rate at MHC-I, which was also supported by targeted sequencing of peptide-binding exons (at least five MHC-I loci genotyped). We found high allelic richness at both MHC-I and MHC-II, although signature of diversifying selection and recombination (gene conversion) was much stronger at MHC-II. Our results indicate that Eurasian coot retains extraordinary polymorphism at both MHC classes (when compared to other non-passerine bird species), although they may be subject to different evolutionary mechanism.

Neutral and Selective Processes Shape MHC Diversity in Roe Deer in Slovenia

Simple Summary

Disease prevention and appropriate wildlife management are among the major challenges in wildlife conservation. In the present study, we made a first assessment of the variability of major histocompatibility complex (MHC) genes in roe deer in Slovenia and evaluated local population adaptation by comparing MHC variability with neutral microsatellites. We discovered three new MHC DRB exon 2 alleles in addition to seven previously described in the literature. Moreover, we found evidence of historical positive selection, as selection analysis indicated that approx. 10% of the encoded amino acids were subjected to episodic positive selection. This study provides the basis for further research on immunogenetic variation in roe deer and highlights opportunities to incorporate genetic data into science-based population management.

Abstract

Disease control and containment in free-ranging populations is one of the greatest challenges in wildlife management. Despite the importance of major histocompatibility complex (MHC) genes for immune response, an assessment of the diversity and occurrence of these genes is still rare in European roe deer, the most abundant and widespread large mammal in Europe. Therefore, we examined immunogenetic variation in roe deer in Slovenia to identify species adaptation by comparing the genetic diversity of the MHC genes with the data on neutral microsatellites. We found ten MHC DRB alleles, three of which are novel. Evidence for historical positive selection on the MHC was found using the maximum likelihood codon method. Patterns of MHC allelic distribution were not congruent with neutral population genetic findings. The lack of population genetic differentiation in MHC genes compared to existing structure in neutral markers suggests that MHC polymorphism was influenced primarily by balancing selection and, to a lesser extent, by neutral processes such as genetic drift, with no clear evidence of local adaptation. Selection analyses indicated that approx. 10% of amino acids encoded under episodic positive selection. This study represents one of the first steps towards establishing an immunogenetic map of roe deer populations across Europe, aiming to better support science-based management of this important game species.

Discovery of an ancient MHC category with both class I and class II features

Two classes of major histocompatibility complex (MHC) molecules, MHC class I and MHC class II, constitute the basis of our elaborate, adaptive immune system as antigen-presenting molecules. They perform distinct, critical functions: especially, MHC class I in case of antivirus and antitumor defenses, and MHC class II, in case of effective antibody responses. This important class diversification has long been enigmatic, as vestiges of the evolutionary molecular changes have not been found. The revealed ancient MHC category represents a plausible intermediate group between the two classes, and the data suggest that class II preceded class I in molecular evolution. Fundamental understanding of the molecular evolution of MHC molecules should contribute to understanding the basis of our complex biological defense system.

Two classes of major histocompatibility complex (MHC) molecules, MHC class I and class II, play important roles in our immune system, presenting antigens to functionally distinct T lymphocyte populations. However, the origin of this essential MHC class divergence is poorly understood. Here, we discovered a category of MHC molecules (W-category) in the most primitive jawed vertebrates, cartilaginous fish, and also in bony fish and tetrapods. W-category, surprisingly, possesses class II–type α- and β-chain organization together with class I–specific sequence motifs for interdomain binding, and the W-category α2 domain shows unprecedented, phylogenetic similarity with β₂-microglobulin of class I. Based on the results, we propose a model in which the ancestral MHC class I molecule evolved from class II–type W-category. The discovery of the ancient MHC group, W-category, sheds a light on the long-standing critical question of the MHC class divergence and suggests that class II type came first.

Structural Comparison Between MHC Classes I and II; in Evolution, a Class-II-Like Molecule Probably Came First

Structures of peptide-loaded major histocompatibility complex class I (pMHC-I) and class II (pMHC-II) complexes are similar. However, whereas pMHC-II complexes include similar-sized IIα and IIβ chains, pMHC-I complexes include a heavy chain (HC) and a single domain molecule β₂-microglobulin (β₂-m). Recently, we elucidated several pMHC-I and pMHC-II structures of primitive vertebrate species. In the present study, a comprehensive comparison of pMHC-I and pMHC-II structures helps to understand pMHC structural evolution and supports the earlier proposed—though debated—direction of MHC evolution from class II-type to class I. Extant pMHC-II structures share major functional characteristics with a deduced MHC-II-type homodimer ancestor. Evolutionary establishment of pMHC-I presumably involved important new functions such as (i) increased peptide selectivity by binding the peptides in a closed groove (ii), structural amplification of peptide ligand sequence differences by binding in a non-relaxed fashion, and (iii) increased peptide selectivity by syngeneic heterotrimer complex formation between peptide, HC, and β₂-m. These new functions were associated with structures that since their establishment in early pMHC-I have been very well conserved, including a shifted and reorganized P1 pocket (aka A pocket), and insertion of a β₂-m hydrophobic knob into the peptide binding domain β-sheet floor. A comparison between divergent species indicates better sequence conservation of peptide binding domains among MHC-I than among MHC-II, agreeing with more demanding interactions within pMHC-I complexes. In lungfishes, genes encoding fusions of all MHC-IIα and MHC-IIβ extracellular domains were identified, and although these lungfish genes presumably derived from classical MHC-II, they provide an alternative mechanistic hypothesis for how evolution from class II-type to class I may have occurred.

The strength of selection is consistent across both domains of the MHC class I peptide-binding groove in birds

Background

The Major Histocompatibility Complex (MHC) codes for the key vertebrate immune receptors responsible for pathogen recognition. Foreign antigens are recognized via their compatibility to hyper-variable region of the peptide-binding groove (PBR), which consists of two separate protein domains. Specifically, the PBR of the MHC class I receptors, which recognize intra-cellular pathogens, has two α domains encoded by exon 2 (α₁) and exon 3 (α₂) of the same gene. Most research on avian MHC class I polymorphism has traditionally focused exclusively on exon 3 and comparisons of selection between the two domains have been hampered by the scarcity of molecular data for exon 2. Thus, it is not clear whether the two domains vary in their specificity towards different antigens and whether they are subject to different selective pressure.

Results

Here, we took advantage of rapidly accumulating genomic resources to test for the differences in selection patterns between both MHC class I domains of the peptide-binding groove in birds. For this purpose, we compiled a dataset of MHC class I exon 2 and 3 sequences for 120 avian species from 46 families. Our phylogenetically-robust approach provided strong evidence for highly consistent levels of selection on the α₁ and α₂ domains. There were strong correlations in all selection measures (number of positively/negatively selected residues and dN/dS ratios) between both PBR exons. Similar positive associations were found for the level of amino acid polymorphism across the two domains.

Conclusions

We conclude that the strength of selection and the level of polymorphism are highly consistent between both peptide-binding domains (α₁ and α₂) of the avian MHC class I.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01812-x.

Genetic variation and selection in the major histocompatibility complex Class II gene in the Guizhou pony

The Guizhou pony (GZP) is an indigenous species of equid found in the mountains of the Guizhou province in southwest China. We selected four regions of the equine leukocyte antigen (ELA), including DQA, DRA, DQB, and DRB, and used them to assess the diversity of the major histocompatibility complex (MHC) class II gene using direct sequencing technology. DRA had the lowest d_N/d_S ratio (0.560) compared with the other three loci, indicating that DRA was conserved and could be conserved after undergoing selective processes. Nine DQA, five DQB, nine DRA, and seven DRB codons were under significant positive selection at the antigen binding sites (ABS), suggesting that the selected residues in ABS may play a significant role in the innate immune system of the GZP. Two GZP alleles were shared with Przewalski’s horse, and six older GZP haplotypes had a better relationship with other horse species by one or two mutational steps, indicating that the GZP may be a natural ancient variety of equid. The specific diversity of ABS and the numbers of unique haplotypes in the evolutionary process affords this species a better genetic fitness and ability to adapt to the native environment.

Allelic diversity and patterns of selection at the major histocompatibility complex class I and II loci in a threatened shorebird, the Snowy Plover (Charadrius nivosus)

Background

Understanding the structure and variability of adaptive loci such as the major histocompatibility complex (MHC) genes is a primary research goal for evolutionary and conservation genetics. Typically, classical MHC genes show high polymorphism and are under strong balancing selection, as their products trigger the adaptive immune response in vertebrates. Here, we assess the allelic diversity and patterns of selection for MHC class I and class II loci in a threatened shorebird with highly flexible mating and parental care behaviour, the Snowy Plover (Charadrius nivosus) across its broad geographic range.

Results

We determined the allelic and nucleotide diversity for MHC class I and class II genes using samples of 250 individuals from eight breeding population of Snowy Plovers. We found 40 alleles at MHC class I and six alleles at MHC class II, with individuals carrying two to seven different alleles (mean 3.70) at MHC class I and up to two alleles (mean 1.45) at MHC class II. Diversity was higher in the peptide-binding region, which suggests balancing selection. The MHC class I locus showed stronger signatures of both positive and negative selection than the MHC class II locus. Most alleles were present in more than one population. If present, private alleles generally occurred at very low frequencies in each population, except for the private alleles of MHC class I in one island population (Puerto Rico, lineage tenuirostris).

Conclusion

Snowy Plovers exhibited an intermediate level of diversity at the MHC, similar to that reported in other Charadriiformes. The differences found in the patterns of selection between the class I and II loci are consistent with the hypothesis that different mechanisms shape the sequence evolution of MHC class I and class II genes. The rarity of private alleles across populations is consistent with high natal and breeding dispersal and the low genetic structure previously observed at neutral genetic markers in this species.

Comparative genetics of the major histocompatibility complex in humans and nonhuman primates

The major histocompatibility complex (MHC) is one of the most gene-dense regions of the mammalian genome. Multiple genes within the human MHC (HLA) show extensive polymorphism, and currently, more than 26,000 alleles divided over 39 different genes are known. Nonhuman primate (NHP) species are grouped into great and lesser apes and Old and New World monkeys, and their MHC is studied mostly because of their important role as animal models in preclinical research or in connection with conservation biology purposes. The evolutionary equivalents of many of the HLA genes are present in NHP species, and these genes may also show abundant levels of polymorphism. This review is intended to provide a comprehensive comparison relating to the organization and polymorphism of human and NHP MHC regions.

HLA-DM Stabilizes the Empty MHCII Binding Groove: A Model Using Customized Natural Move Monte Carlo

MHC class II molecules bind peptides derived from extracellular proteins that have been ingested by antigen-presenting cells and display them to the immune system. Peptide loading occurs within the antigen-presenting cell and is facilitated by HLA-DM. HLA-DM stabilizes the open conformation of the MHCII binding groove when no peptide is bound. While a structure of the MHCII/HLA-DM complex exists, the mechanism of stabilization is still largely unknown. Here, we applied customized Natural Move Monte Carlo to investigate this interaction. We found a possible long-range mechanism that implicates the configuration of the membrane-proximal globular domains in stabilizing the open state of the empty MHCII binding groove.

Molecular genetic characterization and haplotype diversity of swine leukocyte antigen in Chinese Rongshui miniature pigs

The Rongshui miniature pig is an important model animal for studying livestock disease prevention and control in China. The highly polymorphic swine leukocyte antigen (SLA) has been the focus of considerable interest because of the strong, reproducible associations between particular SLA haplotypes and infectious diseases. In this study, we identified 42 alleles at eight polymorphic SLA loci (SLA-1, SLA-3, SLA-2, SLA-6, DRA, DRB1, DQA, and DQB1) representing seven class I and six class II haplotypes using reverse transcription-polymerase chain reaction (RT-PCR) sequence-based typing and PCR-sequence specific primers in Rongshui miniature pigs. The official names were designated by the SLA Nomenclature Committee of the International Society for Animal Genetics. Seven class I haplotypes, Hp-5b.0, 86.0, 87.0, 88.0, 89.0, 90.0 and 91.0, and four class II haplotypes, Hp-0.18b, 0.19c, 0.41 and 0.47, had not previously been reported in other pig breeds. We also comprehensively analyzed the molecular genetic characterization and phylogenies of the identified alleles and the SLA haplotype diversity in Rongshui miniature pigs. SLA-1 and SLA-6 genes were under positive selection, while SLA-2 was under neutral selection, and the other five genes were under purifying selection. The highly polymorphic new alleles may be derived by nucleotide mutations, insertions and deletions, and fragment recombination, and alleles segregated based on sequence differences and peptide-binding motifs, rather than on pig breed. SLA haplotype diversity was generated by allele/gene conversion and recombination. These results will be helpful for elucidating the molecular genetic mechanisms influencing differential disease resistance among pigs with different SLA haplotypes.

HLA alleles, especially amino-acid signatures of HLA-DPB1, might contribute to the molecular pathogenesis of early-onset autoimmune thyroid disease

The major histocompatibility complex region has been suggested to play an important role in the development of autoimmune thyroid disease (AITD). In this study, we investigated the associations of human leukocyte antigen (HLA) alleles and amino acid variants of HLA with early-onset AITD. HLA class I and class II genes were analyzed in 116 Korean children with AITDs (Graves’ disease [GD]: 71, Hashimoto’s disease [HD]: 45) and 142 healthy controls. HLA-B*46:01 (OR = 3.96, Pc = 0.008), -C*01:02 (OR = 2.51 Pc = 0.04), -DPB1*02:02 (OR = 3.99, Pc = 0.04), and -DPB1*05:01 (OR = 4.6, Pc = 0.003) were significantly associated with GD, and HLA-A*02:07 (OR = 4.68, Pc = 0.045) and -DPB1*02:02 (OR = 6.57, Pc = 0.0001) with HD. The frequency of HLA-DPB1*05:01 was significantly higher in GD patients than in HD patients (Pc = 0.0005). Furthermore, differences were found between patients with Thyroid associated ophthalmopathy (TAO) and those without TAO in the distribution of HLA-B*54:01 (8.6% vs. 30.6%, P = 0.04) and -C*03:03 (37.1% vs. 11.1%, P = 0.02). In the analysis of amino acid variants of HLA molecules, both Leu35 (OR = 23.38, P = 0.0002) and Glu55 (OR = 23.38, P = 0.0002) of HLA-DPB1 were strongly associated with GD and showed different distributions between GD and HD (P = 0.001). Our results suggest that HLA alleles, especially amino-acid signatures of the HLA-DP β chain, might contribute to the molecular pathogenesis of early-onset AITD.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore

Background

Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci.

Results

Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection.

Conclusion

Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1420-5) contains supplementary material, which is available to authorized users.

A global analysis of selection at the avian MHC

Recent advancements in sequencing technology have resulted in rapid progress in the study of the major histocompatibility complex (MHC) in non-model avian species. Here, we analyze a global dataset of avian MHC class I and class II sequences (ca. 11,000 sequences from over 250 species) to gain insight into the processes that govern macroevolution of MHC genes in birds. Analysis of substitution rates revealed striking differences in the patterns of diversifying selection between passerine and non-passerine birds. Non-passerines showed stronger selection at MHC class II, which is primarily involved in recognition of extracellular pathogens, while passerines showed stronger selection at MHC class I, which is involved in recognition of intracellular pathogens. Positions of positively selected amino-acid residues showed marked discrepancies with peptide-binding residues (PBRs) of human MHC molecules, suggesting that using a human classification of PBRs to assess selection patterns at the avian MHC may be unjustified. Finally, our analysis provided evidence that indel mutations can make a substantial contribution to adaptive variation at the avian MHC.

Low genetic variation in the MHC class II DRB gene and MHC-linked microsatellites in endangered island populations of the leopard cat (Prionailurus bengalensis) in Japan

Isolated populations of the leopard cat (Prionailurus bengalensis) on Tsushima and Iriomote islands in Japan are classified as subspecies P. b. euptilurus and P. b. iriomotensis, respectively. Because both populations have decreased to roughly 100, an understanding of their genetic diversity is essential for conservation. We genotyped MHC class II DRB exon 2 and MHC-linked microsatellite loci to evaluate the diversity of MHC genes in the Tsushima and Iriomote cat populations. We detected ten and four DRB alleles in these populations, respectively. A phylogenetic analysis showed DRB alleles from both populations to be closely related to those in other felid DRB lineages, indicating trans-species polymorphism. The MHC-linked microsatellites were more polymorphic in the Tsushima than in the Iriomote population. The MHC diversity of both leopard cat populations is much lower than in the domestic cat populations on these islands, probably due to inbreeding associated with founder effects, geographical isolation, or genetic drift. Our results predict low resistance of the two endangered populations to new pathogens introduced to the islands.

Bonobos Maintain Immune System Diversity with Three Functional Types of MHC-B

Fast-evolving MHC class I polymorphism serves to diversify NK cell and CD8 T cell responses in individuals, families, and populations. Because only chimpanzee and bonobo have strict orthologs of all HLA class I, their study gives unique perspectives on the human condition. We defined polymorphism of Papa-B, the bonobo ortholog of HLA-B, for six wild bonobo populations. Sequences for Papa-B exon 2 and 3 were determined from the genomic DNA in 255 fecal samples, minimally representing 110 individuals. Twenty-two Papa-B alleles were defined, each encoding a different Papa-B protein. No Papa-B is identical to any chimpanzee Patr-B, human HLA-B, or gorilla Gogo-B. Phylogenetic analysis identified a clade of MHC-B, defined by residues 45-74 of the α₁ domain, which is broadly conserved among bonobo, chimpanzee, and gorilla. Bonobo populations have 3-14 Papa-B allotypes. Three Papa-B are in all populations, and they are each of a different functional type: allotypes having the Bw4 epitope recognized by killer cell Ig-like receptors of NK cells, allotypes having the C1 epitope also recognized by killer cell Ig-like receptors, and allotypes having neither epitope. For population Malebo, these three Papa-B are the only Papa-B allotypes. Although small in number, their sequence divergence is such that the nucleotide diversity (mean proportional distance) of Papa-B in Malebo is greater than in the other populations and is also greater than expected for random combinations of three Papa-B Overall, Papa-B has substantially less diversity than Patr-B in chimpanzee subspecies and HLA-B in indigenous human populations, consistent with bonobo having experienced narrower population bottlenecks.

Evidence for selection maintaining MHC diversity in a rodent species despite strong density fluctuations

Strong spatiotemporal variation in population size often leads to reduced genetic diversity limiting the adaptive potential of individual populations. Key genes of adaptive variation are encoded by the immune genes of the major histocompatibility complex (MHC) playing an essential role in parasite resistance. How MHC variation persists in rodent populations that regularly experience population bottlenecks remains an important topic in evolutionary genetics. We analysed the consequences of strong population fluctuations on MHC class II DRB exon 2 diversity in two distant common vole (Microtus arvalis) populations in three consecutive years using a high-throughput sequencing approach. In 143 individuals, we detected 25 nucleotide alleles translating into 14 unique amino acid MHC alleles belonging to at least three loci. Thus, the overall allelic diversity and amino acid distance among the remaining MHC alleles, used as a surrogate for the range of pathogenic antigens that can be presented to T-cells, are still remarkably high. Both study populations did not show significant population differentiation between years, but significant differences were found between sites. We concluded that selection processes seem to be strong enough to maintain moderate levels of MHC diversity in our study populations outcompeting genetic drift, as the same MHC alleles were conserved between years. Differences in allele frequencies between populations might be the outcome of different local parasite pressures and/or genetic drift. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in cyclic populations.

Organization and characteristics of the major histocompatibility complex class II region in the Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis)

Little is known about the major histocompatibility complex (MHC) in the genome of Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis) (YFP) or other cetaceans. In this study, a high-quality YFP bacterial artificial chromosome (BAC) library was constructed. We then determined the organization and characterization of YFP MHC class II region by screening the BAC library, followed by sequencing and assembly of positive BAC clones. The YFP MHC class II region consists of two segregated contigs (218,725 bp and 328,435 bp respectively) that include only eight expressed MHC class II genes, three pseudo MHC genes and twelve non-MHC genes. The YFP has fewer MHC class II genes than ruminants, showing locus reduction in DRB, DQA, DQB, and loss of DY. In addition, phylogenic and evolutionary analyses indicated that the DRB, DQA and DQB genes might have undergone birth-and-death evolution, whereas the DQB gene might have evolved under positive selection in cetaceans. These findings provide an essential foundation for future work, such as estimating MHC genetic variation in the YFP or other cetaceans. This work is the first report on the MHC class II region in cetaceans and offers valuable information for understanding the evolution of MHC genome in cetaceans.

Unconventional T-cell recognition of an arthritogenic epitope of proteoglycan aggrecan released from degrading cartilage

It has been proposed that peptide epitopes bind to MHC class II molecules to form distinct structural conformers of the same MHC II-peptide complex termed type A and type B, and that the two conformers of the same peptide-MHC II complex are recognized by distinct CD4 T cells, termed type A and type B T cells. Both types recognize short synthetic peptides but only type A recognize endosomally processed intact antigen. Type B T cells that recognize self peptides from exogenously degraded proteins have been shown to escape negative selection during thymic development and so have the potential to contribute to the pathogenesis of autoimmunity. We generated and characterized mouse CD4 T cells specific for an arthritogenic epitope of the candidate joint autoantigen proteoglycan aggrecan. Cloned T-cell hybridomas specific for a synthetic peptide containing the aggrecan epitope showed two distinct response patterns based on whether they could recognize processed intact aggrecan. Fine mapping demonstrated that both types of T-cell recognized the same core epitope. The results are consistent with the generation of aggrecan-specific type A and type B T cells. Type B T cells were activated by supernatants released from degrading cartilage, indicating the presence of antigenic extracellular peptides or fragments of aggrecan. Type B T cells could play a role in the pathogenesis of proteoglycan-induced arthritis in mice, a model for rheumatoid arthritis, by recognizing extracellular peptides or protein fragments of joint autoantigens released by inflamed cartilage.

Signature Patterns of MHC Diversity in Three Gombe Communities of Wild Chimpanzees Reflect Fitness in Reproduction and Immune Defense against SIVcpz

Major histocompatibility complex (MHC) class I molecules determine immune responses to viral infections. These polymorphic cell-surface glycoproteins bind peptide antigens, forming ligands for cytotoxic T and natural killer cell receptors. Under pressure from rapidly evolving viruses, hominoid MHC class I molecules also evolve rapidly, becoming diverse and species-specific. Little is known of the impact of infectious disease epidemics on MHC class I variant distributions in human populations, a context in which the chimpanzee is the superior animal model. Population dynamics of the chimpanzees inhabiting Gombe National Park, Tanzania have been studied for over 50 years. This population is infected with SIVcpz, the precursor of human HIV-1. Because HLA-B is the most polymorphic human MHC class I molecule and correlates strongly with HIV-1 progression, we determined sequences for its ortholog, Patr-B, in 125 Gombe chimpanzees. Eleven Patr-B variants were defined, as were their frequencies in Gombe's three communities, changes in frequency with time, and effect of SIVcpz infection. The growing populations of the northern and central communities, where SIVcpz is less prevalent, have stable distributions comprising a majority of low-frequency Patr-B variants and a few high-frequency variants. Driving the latter to high frequency has been the fecundity of immigrants to the northern community, whereas in the central community, it has been the fecundity of socially dominant individuals. In the declining population of the southern community, where greater SIVcpz prevalence is associated with mortality and emigration, Patr-B variant distributions have been changing. Enriched in this community are Patr-B variants that engage with natural killer cell receptors. Elevated among SIVcpz-infected chimpanzees, the Patr-B*06:03 variant has striking structural and functional similarities to HLA-B*57, the human allotype most strongly associated with delayed HIV-1 progression. Like HLA-B*57, Patr-B*06:03 correlates with reduced viral load, as assessed by detection of SIVcpz RNA in feces.

T Cells: Soldiers and Spies--The Surveillance and Control of Effector T Cells by Regulatory T Cells

Traditionally, T cells were CD4+ helper or CD8+ cytotoxic T cells, and with antibodies, they were the soldiers of immunity. Now, many functionally distinct subsets of activated CD4+ and CD8+ T cells have been described, each with distinct cytokine and transcription factor expression. For CD4+ T cells, these include Th1 cells expressing the transcription factor T-bet and cytokines IL-2, IFN-γ, and TNF-β; Th2 cells expressing GATA-3 and the cytokines IL-4, IL-5, and IL-13; and Th17 cells expressing RORγt and cytokines IL-17A, IL-17F, IL-21, and IL-22. The cytokines produced determine the immune inflammation that they mediate. T cells of the effector lineage can be naïve T cells, recently activated T cells, or memory T cells that can be distinguished by cell surface markers. T regulatory cells or spies were characterized as CD8+ T cells expressing I-J in the 1970s. In the 1980s, suppressor cells fell into disrepute when the gene for I-J was not present in the mouse MHC I region. At that time, a CD4+ T cell expressing CD25, the IL-2 receptor-α, was identified to transfer transplant tolerance. This was the same phenotype of activated CD4+ CD25+ T cells that mediated rejection. Thus, the cells that could induce tolerance and undermine rejection had similar badges and uniforms as the cells effecting rejection. Later, FOXP3, a transcription factor that confers suppressor function, was described and distinguishes T regulatory cells from effector T cells. Many subtypes of T regulatory cells can be characterized by different expressions of cytokines and receptors for cytokines or chemokines. In intense immune inflammation, T regulatory cells express cytokines characteristic of effector cells; for example, Th1-like T regulatory cells express T-bet, and IFN-γ-like Th1 cells and effector T cells can change sides by converting to T regulatory cells. Effector T cells and T regulatory cells use similar molecules to be activated and mediate their function, and thus, it can be very difficult to distinguish soldiers from spies.

Maintaining microendemic primate species along an environmental gradient - parasites as drivers for species differentiation

Understanding the drivers of species adaptations to changing environments on the one hand and the limits for hybridization on the other hand is among the hottest questions in evolutionary biology. Parasites represent one of the major selective forces driving host evolution and at least those with free-living stages are at the same time dependent on the ecological conditions of their host's habitat. Local immunological adaptations of host species to varying parasite pressure are therefore expected and might represent the genetic basis for ecological speciation and the maintenance of recently diverged species. Madagascar provides one of the rare examples where two partially sympatric primate species (Microcebus griseorufus, M. murinus) and their hybrids, as well as an allopatric species (M. cf rufus) live in close proximity along a very steep environmental gradient ranging from southern dry spiny bush to gallery forest to evergreen eastern humid rain forest, thus mimicking the situation encountered during extensions and retreats of vegetation formations under changing climatic conditions. This system was used to study parasite infection and immune gene (MHC) adaptations to varying parasite pressure that might provide selective advantages to pure species over hybrids. Parasite burdens increased with increasing humidity. M. griseorufus, M. murinus, and their hybrids but not M. rufus shared the same MHC alleles, indicating either retention of ancestral polymorphism or recent gene flow. The hybrids had much higher prevalence of intestinal parasites than either of the parent species living under identical environmental conditions. The different representation of parasites can indicate a handicap for hybrids that maintains species identities.

Characterization of swine leukocyte antigen (SLA) polymorphism by sequence-based and PCR-SSP methods in Chinese Bama miniature pigs

The highly polymorphic swine leukocyte antigen (SLA) genes have been repeatedly shown to influence swine immune traits, disease resistance, vaccine responsiveness and tumour penetrance. Analysis of the SLA diversity in as many pig breeds as possible is important to clarify the relationships between SLA genes and diseases or traits, and develop these pigs as valuable animal models for biomedical research. The Chinese Bama miniature pig breed is an economically significant breed that is available at several research institutions in China. In this study, we identified a total of 32 alleles at five polymorphic SLA loci (SLA-1, SLA-3, SLA-2, DRB1 and DQB1) representing nine class I and seven class II haplotypes using the reverse transcription polymerase chain reaction (RT-PCR) sequence-based typing (SBT) method. The possible functional sites of the SLA genes were predicted and analyzed by comparison with those of the human and mouse. Based on the sequence information, we subsequently developed a rapid PCR-based typing assay using sequence-specific primers (PCR-SSP) to efficiently follow the SLA types of the progeny. In the studied cohort (2n = 562), the most prevalent Haplotype Hp-35.6 (SLA-1(∗)1201, SLA-1(∗)1301-SLA-3(∗)0502-SLA-2(∗)1001-DRB1(∗)0501-DQB1(∗)0801) was identified in 182 Bama pigs with a frequency of 32.38%. The presence of the duplicated SLA-1 locus was confirmed in five of the class I haplotypes. Moreover, we identified two crossovers within the class I region and one between the class I and class II regions, which corresponded to recombination frequencies of 0.36% and 0.18%, respectively. The information of this study is essential for an understanding of the SLA allelic architecture and diversity, and it will be helpful for studying the adaptive immune response and further developing the more effective vaccines in the context of SLA specificities.

Development of MHC-Linked Microsatellite Markers in the Domestic Cat and Their Use to Evaluate MHC Diversity in Domestic Cats, Cheetahs, and Gir Lions

Diversity within the major histocompatibility complex (MHC) reflects the immunological fitness of a population. MHC-linked microsatellite markers provide a simple and an inexpensive method for studying MHC diversity in large-scale studies. We have developed 6 MHC-linked microsatellite markers in the domestic cat and used these, in conjunction with 5 neutral microsatellites, to assess MHC diversity in domestic mixed breed (n = 129) and purebred Burmese (n = 61) cat populations in Australia. The MHC of outbred Australian cats is polymorphic (average allelic richness = 8.52), whereas the Burmese population has significantly lower MHC diversity (average allelic richness = 6.81; P < 0.01). The MHC-linked microsatellites along with MHC cloning and sequencing demonstrated moderate MHC diversity in cheetahs (n = 13) and extremely low diversity in Gir lions (n = 13). Our MHC-linked microsatellite markers have potential future use in diversity and disease studies in other populations and breeds of cats as well as in wild felid species.

Autoantibody and human leukocyte antigen profiles in children with autoimmune liver disease and their first-degree relatives

OBJECTIVE

Familial clustering of juvenile autoimmune liver disease (AILD), including autoimmune hepatitis and autoimmune sclerosing cholangitis (ASC), is rare, despite a high prevalence of autoimmune disorders in AILD families.

METHODS

To investigate this discrepancy, we measured autoantibodies diagnostic for AILD, anti-nuclear, anti-smooth muscle, anti-liver kidney microsomal type 1, anti-liver cytosol type 1, and anti-soluble liver antigen antibodies, and human leukocyte antigen profiles in 31 patients and 65 of their first-degree relatives (FDR). The autoantibody profile was compared with that of 42 healthy subjects (HS).

RESULTS

Autoantibodies were detected in 71% (22/31) patients. Anti-nuclear antibody or anti-smooth muscle antibody were present in 4/65 FDR (6.2%). HS were negative for all autoantibodies. The frequencies of homozygous HLA DRB1*0301 (DR3) genes and haplotype A1-B8-DR3 were higher in the patients (25% and 43%) than in FDR (9% and 27%) and HS (0% and 16%). The frequencies of disease-protective genes DR4 and/or DR15 were lower in the patients (25%) than in FDR (42%) and HS (42%). Only 1 family contained 2 patients with AILD, 1 with ASC and 1 with primary sclerosing cholangitis. Both patients possessed A1-B8-DR3 genes, the ASC being homozygous and the primary sclerosing cholangitis heterozygous. Six FDR had nonhepatic autoimmune disorders, none being autoantibody positive.

CONCLUSIONS

Homozygosity for DR3 plays a major role in the predisposition to juvenile AILD. Diagnostic autoantibodies for AILD are rare among patients' FDR and not linked to clinical manifestation of AILD.

Extraordinary MHC class II B diversity in a non-passerine, wild bird: the Eurasian Coot Fulica atra (Aves: Rallidae)

The major histocompatibility complex (MHC) hosts the most polymorphic genes ever described in vertebrates. The MHC triggers the adaptive branch of the immune response, and its extraordinary variability is considered an evolutionary consequence of pathogen pressure. The last few years have witnessed the characterization of the MHC multigene family in a large diversity of bird species, unraveling important differences in its polymorphism, complexity, and evolution. Here, we characterize the first MHC class II B sequences isolated from a Rallidae species, the Eurasian Coot Fulica atra. A next-generation sequencing approach revealed up to 265 alleles that translated into 251 different amino acid sequences (β chain, exon 2) in 902 individuals. Bayesian inference identified up to 19 codons within the presumptive peptide-binding region showing pervasive evidence of positive, diversifying selection. Our analyses also detected a significant excess of high-frequency segregating sites (average Tajima's D = 2.36, P < 0.05), indicative of balancing selection. We found one to six different alleles per individual, consistent with the occurrence of at least three MHC class II B gene duplicates. However, the genotypes comprised of three alleles were by far the most abundant in the population investigated (49.4%), followed by those with two (29.6%) and four (17.5%) alleles. We suggest that these proportions are in agreement with the segregation of MHC haplotypes differing in gene copy number. The most widespread segregating haplotypes, according to our findings, would contain one single gene or two genes. The MHC class II of the Eurasian Coot is a valuable system to investigate the evolutionary implications of gene copy variation and extensive variability, the greatest ever found, to the best of our knowledge, in a wild population of a non-passerine bird.

MHC class I recognition by monocyte-/macrophage-specific receptors

The most important transplantation antigens in the discrimination between "self" and "nonself" are encoded by genes in the major histocompatibility complex (MHC) locus (H-2 in mice). It has been assumed that T lymphocytes are the effector cells for allograft rejection, as athymic nude rodents fail to reject allografts. In 1988, we i.p. transplanted Meth A (H-2D(d)K(d)) tumor cells into C57BL/6 (H-2D(b)K(b)) mice and found macrophages to be cytotoxic against the allografts. In 1996, several groups using CD4 or CD8 knockout mice reported that non-T cells were the effector cells for the rejection of skin or organ allografts. In 1998, we ascertained that macrophages were the effector cells of skin allograft rejection. Recently, we isolated cDNA clones encoding monocyte/macrophage MHC receptors (MMRs) for H-2D(d) and H-2K(d); established H-2D(d)- and/or H-2K(d)-transgenic mice and lymphoma cells; and found, using MMR-deficient mice, that MMR and T-cell receptor were essential for the rejection of transgenic skin and lymphoma, respectively.

Comprehensive analysis of medaka major histocompatibility complex (MHC) class II genes: implications for evolution in teleosts

The major histocompatibility complex (MHC) class II molecules play central roles in adaptive immunity by regulating immune response via the activation of CD4 T cells. The full complement of the MHC class II genes has been elucidated only in mammalian species to date. To understand the evolution of these genes, we performed their first comprehensive analysis in nonmammalian species using a teleost, medaka (Oryzias latipes). Based on a database search, cDNA cloning, and genomic PCR, medaka was shown to possess five pairs of expressed class II genes, comprising one IIA and one IIB gene. Each pair was located on a different chromosome and was not linked to the class I genes. Only one pair showed a high degree of polymorphism and was considered to be classical class II genes, whereas the other four pairs were nonclassical. Phylogenetic analysis of all medaka class II genes and most reported teleost class II genes revealed that the IIA and IIB genes formed separate clades, each containing three well-corresponding lineages. One lineage contained three medaka genes and all known classical class II genes of Ostariophysi and Euteleostei and was presumed to be an original lineage of the teleost MHC class II genes. The other two lineages contained one nonclassical medaka gene each and some Euteleostei genes. These results indicate that multiple lineages of the teleost MHC class II genes have been conserved for hundreds of millions of years and that the tightly linked IIA and IIB genes have undergone concerted evolution.

MHC variability in an isolated wolf population in Italy

Small, isolated populations may experience increased extinction risk due to reduced genetic variability at important functional genes, thus decreasing the population's adaptive potential. The major histocompatibility complex (MHC), a key immunological gene cluster, usually shows high variability maintained by positive or balancing selection in response to challenges by pathogens. Here we investigated for the first time, the variability of 3 MHC class II genes (DRB1, DQA1, and DQB1) in 94 samples collected from Italian wolves. The Italian wolf population has been long isolated south of the Alps and is presently recovering from a recent bottleneck that decreased the population to less than 100 individuals. Despite the bottleneck, Italian wolves show remarkable MHC variability with 6-9 alleles per locus, including 2 recently described alleles at DRB1. MHC sequences show signatures of historical selective pressures (high d N/d S ratio, ω > 1.74) but no evidence of ongoing selection. Variation at the MHC genes and 12 background microsatellite loci were not apparently affected by the recent bottleneck. Although MHC alleles of domestic dog origin were detected in 8 genetically admixed individuals, these alleles were rare or absent in nonadmixed wolves. Thus, despite known hybridization events between domestic dogs and Italian wolves, the Italian wolf population does not appear affected by deep introgression of domestic dog MHC alleles.

Comprehensive identification of high-frequency and co-occurring Mafa-DPA1, Mafa-DQA1, Mafa-DRA, and Mafa-DOA alleles in Vietnamese cynomolgus macaques

High-frequency alleles and/or co-occurring human leukocyte antigen alleles across loci appear to be more important than individual alleles as markers of disease risk and have clinical value as biomarkers for targeted screening or the development of new disease therapies. To better elucidate the major histocompatibility complex (MHC) background and to facilitate the experimental use of cynomolgus macaques, Mafa-DPA1, Mafa-DQA1, Mafa-DRA, and Mafa-DOA alleles were characterized, and their combinations were investigated in 30 Vietnamese macaques by gene cloning and sequencing. A total of 26 Mafa-DPA1, 18 Mafa-DQA1, 9 Mafa-DRA, and 15 Mafa-DOA alleles, including 7 high-frequency alleles, were identified in this study, respectively. In addition, 15 Mafa-DQA1, 17 Mafa-DPA1, 15 Mafa-DOA, and 2 Mafa-DRA alleles represented novel sequences that had not been documented in earlier studies. Our results also showed that the Vietnamese macaques might be valuable because no less than 30% of the test animals possessed Mafa-DRA*01:02:01 (90%), -DQA1*26:01:03 (37%), -DOA*01:02:07 (34%), and -DQA1*01:03:03 (30%). We previously reported that the combinations of MHC class II alleles, including the combination of DOA*01:02:07-DPA1*02:09 and DOA*01:02:07-DQA1*01:03:03, were detected in 17 and 14% of the animals, respectively. Interestingly, more than two Mafa-DQA1 and Mafa-DPA1 alleles were detected in one animal in this study, which suggested that they might be caused by a chromosomal duplication. If our findings can be validated by other studies, it will further enrich the number of known Mafa-DPA1 and Mafa-DQA1 polymorphisms. Our results identified the co-occurring MHC alleles across loci in a cohort of Vietnamese cynomolgus macaques, which emphasized the value of this species as a model for biomedical research.

For many but not for all: how the conformational flexibility of the peptide/MHCII complex shapes epitope selection

The adaptive immune response starts when CD4+ T cells recognize peptide antigens presented by class II molecules of the Major Histocompatibility Complex (MHCII). Two outstanding features of MHCII molecules are their polymorphism and the ability of each allele to bind a large panoply of peptides. The ability of each MHCII molecule to interact with a limited, though broad, range of amino acid sequences, or "permissive specificity" of binding, is the result of structural flexibility. This flexibility has been identified through biochemical and biophysical studies, and molecular dynamic simulations have modeled the conformational rearrangements that the peptide and the MHCII undergo during interaction. Moreover, there is evidence that the structural flexibility of the peptide/MHCII complex correlates with the activity of the "peptide-editing" molecule DM. In light of the impact that these recent findings have on our ability to predict MHCII epitopes, a review of the structural and thermodynamic determinants of peptide binding to MHCII is proposed.

Extensive polymorphism of the major histocompatibility complex DRA gene in Balkan donkeys: perspectives on selection and genealogy

The major histocompatibility complex (MHC) contains genes important for immune response in mammals, and these genes exhibit high polymorphism and diversity. The DRA gene, a member of the MHC class II family, is highly conserved across a large number of mammalian species, but it displays exceptionally rich sequence variations in Equidae members. We analyzed allelic polymorphism of the DRA locus in 248 donkeys sampled across the Balkan Peninsula (Albania, Bulgaria, Croatia, Macedonia, Greece and Montenegro). Five known alleles and two new alleles were identified. The new allele Eqas-DRA*0601 was found to carry a synonymous mutation, and new allele Eqas-DRA*0701, a non-synonymous mutation. We further analyzed the historical selection and allele genealogy at the DRA locus in equids. Signals of positive selection obtained by various tests were ambiguous. A conservative conclusion is that DRA polymorphism occurred relatively recently and that positive selection has been acting on the DRA locus for a relatively brief period.

Evolution of major histocompatibility complex class I and class II genes in the brown bear

BACKGROUND

Major histocompatibility complex (MHC) proteins constitute an essential component of the vertebrate immune response, and are coded by the most polymorphic of the vertebrate genes. Here, we investigated sequence variation and evolution of MHC class I and class II DRB, DQA and DQB genes in the brown bear Ursus arctos to characterise the level of polymorphism, estimate the strength of positive selection acting on them, and assess the extent of gene orthology and trans-species polymorphism in Ursidae.

RESULTS

We found 37 MHC class I, 16 MHC class II DRB, four DQB and two DQA alleles. We confirmed the expression of several loci: three MHC class I, two DRB, two DQB and one DQA. MHC class I also contained two clusters of non-expressed sequences. MHC class I and DRB allele frequencies differed between northern and southern populations of the Scandinavian brown bear. The rate of nonsynonymous substitutions (dN) exceeded the rate of synonymous substitutions (dS) at putative antigen binding sites of DRB and DQB loci and, marginally significantly, at MHC class I loci. Models of codon evolution supported positive selection at DRB and MHC class I loci. Both MHC class I and MHC class II sequences showed orthology to gene clusters found in the giant panda Ailuropoda melanoleuca.

CONCLUSIONS

Historical positive selection has acted on MHC class I, class II DRB and DQB, but not on the DQA locus. The signal of historical positive selection on the DRB locus was particularly strong, which may be a general feature of caniforms. The presence of MHC class I pseudogenes may indicate faster gene turnover in this class through the birth-and-death process. South-north population structure at MHC loci probably reflects origin of the populations from separate glacial refugia.

DLA class II alleles and haplotypes are associated with risk for and protection from chronic hepatitis in the English Springer spaniel

Chronic hepatitis (CH) is common in dogs in the United Kingdom. An increased prevalence of the disease is seen in the English Springer spaniel (ESS), and this breed suffer from a severe form with young to middle aged female dogs being predisposed. The disease shares histological features with those of human viral hepatitis, although the specific aetiological agent has not yet been identified. The aim of the current study was to investigate whether dog leucocyte antigen (DLA) class II alleles and haplotypes are associated with susceptibility/resistance to CH in the ESS. Sequence-based genotyping of the polymorphic exon 2 from DLA-DRB1, -DQA1 and -DQB1 class II loci were performed in 66 ESSs with CH and 84 healthy controls. There was a significant difference in the distribution of the protective alleles DRB1*00501 (3.0% vs. 12.0%, odds ratio [OR] = 0.23, 95% confidence interval [CI] = 0.06–0.74) and DQB1*00501 (3.8% vs. 12.0%, OR = 0.29, 95% CI = 0.09–0.85) between cases and controls. The haplotype DLA-DRB1*00501/DQA1*00301/DQB1*00501 was present in 11.9% of controls and 3.0% of cases and was significantly associated with protection against disease development (OR = 0.26, 95% CI = 0.08–0.80). There was a significant difference in the distribution of the risk alleles DRB1*00601 (14.4% vs. 6.5%, OR = 2.40, 95% CI = 1.10–5.63) and DQB1*00701 (14.4% vs. 6.5%, OR = 2.40, 95% CI = 1.10–5.63) between cases and controls. A risk haplotype (DLA-DRB1*00601/DQA1*005011/DQB1*00701) was present in 14.4% of cases and 6.5% of controls and conferred an elevated risk of developing CH with an OR of 3.13 (95% CI = 1.20–8.26). These results demonstrate that DLA class II is significantly associated with risk and protection from developing CH in ESSs.

Molecular characterization of MHC-DRB cDNA in water buffalo (Bubalus bubalis)

In the present study, water buffalo MHC (Bubu)-DRB cDNA was cloned and characterized. The 1022 base long-amplified cDNA product encompassed a single open reading frame of 801 bases that coded for 266 amino acids. The Bubu-DRB sequence showed maximum homology with the BoLA-DRB3*0101 allele of cattle. A total of seven amino acid residues were found to be unique for the Bubu-DRB sequence. The majority of amino acid substitutions was observed in the β₁ domain. Residues associated with important functions were mostly conserved. Water buffalo DRB was phylogenetically closer to goat DRB*A.

Impact of historical founder effects and a recent bottleneck on MHC variability in Commander Arctic foxes (Vulpes lagopus)

Populations of Arctic foxes (Vulpes lagopus) have been isolated on two of the Commander Islands (Bering and Mednyi) from the circumpolar distributed mainland population since the Pleistocene. In 1970-1980, an epizootic outbreak of mange caused a severe population decline on Mednyi Island. Genes of the major histocompatibility complex (MHC) play a primary role in infectious disease resistance. The main objectives of our study were to compare contemporary variation of MHC class II in mainland and island Arctic foxes, and to document the effects of the isolation and the recent bottleneck on MHC polymorphism by analyzing samples from historical and contemporary Arctic foxes. In 184 individuals, we found 25 unique MHC class II DRB and DQB alleles, and identified evidence of balancing selection maintaining allelic lineages over time at both loci. Twenty different MHC alleles were observed in mainland foxes and eight in Bering Island foxes. The historical Mednyi population contained five alleles and all contemporary individuals were monomorphic at both DRB and DQB. Our data indicate that despite positive and diversifying selection leading to elevated rates of amino acid replacement in functionally important antigen-binding sites, below a certain population size, balancing selection may not be strong enough to maintain genetic diversity in functionally important genes. This may have important fitness consequences and might explain the high pathogen susceptibility in some island populations. This is the first study that compares MHC diversity before and after a bottleneck in a wild canid population using DNA from museum samples.