HLAs, TCRs, and KIRs, a Triumvirate of Human Cell-Mediated Immunity.

In all human cells, human leukocyte antigen (HLA) class I glycoproteins assemble with a peptide and take it to the cell surface for surveillance by lymphocytes. These include natural killer (NK) cells and γδ T cells of innate immunity and αß T cells of adaptive immunity. In healthy cells, the presented peptides derive from human proteins, to which lymphocytes are tolerant. In pathogen-infected cells, HLA class I expression is perturbed. Reduced HLA class I expression is detected by KIR and CD94:NKG2A receptors of NK cells. Almost any change in peptide presentation can be detected by αß CD8+ T cells. In responding to extracellular pathogens, HLA class II glycoproteins, expressed by specialized antigen-presenting cells, present peptides to αß CD4+ T cells. In comparison to the families of major histocompatibility complex (MHC) class I, MHC class II and αß T cell receptors, the antigenic specificity of the γδ T cell receptors is incompletely understood.

ELM-MHC: An Improved MHC Identification Method with Extreme Learning Machine Algorithm.

The major histocompatibility complex (MHC) is a term for all gene groups of a major histocompatibility antigen. It binds to peptide chains derived from pathogens and displays pathogens on the cell surface to facilitate T-cell recognition and perform a series of immune functions. MHC molecules are critical in transplantation, autoimmunity, infection, and tumor immunotherapy. Combining machine learning algorithms and making full use of bioinformatics analysis technology, more accurate recognition of MHC is an important task. The paper proposed a new MHC recognition method compared with traditional biological methods and used the built classifier to classify and identify MHC I and MHC II. The classifier used the SVMProt 188D, bag-of-ngrams (BonG), and information theory (IT) mixed feature representation methods and used the extreme learning machine (ELM), which selects lin-kernel as the activation function and used 10-fold cross-validation and the independent test set validation to verify the accuracy of the constructed classifier and simultaneously identify the MHC and identify the MHC I and MHC II, respectively. Through the 10-fold cross-validation, the proposed algorithm obtained 91.66% accuracy when identifying MHC and 94.442% accuracy when identifying MHC categories. Furthermore, an online identification Web site named ELM-MHC was constructed with the following URL: http://server.malab.cn/ELM-MHC/ .

Allelic nomenclature for the duplicated MHC class II DQ genes in sheep.

The principal MHC class II molecules involved in the presentation of peptides to the antigen specific receptors on CD4+ T cells genes in sheep are derived from DR and DQ genes. Allelic nomenclature systems for the DRB1 and its partner DRA loci are available for Ovid's; however, no official nomenclature is available for the DQ genes which creates ambiguity within the research community. Ovine MHC haplotypes include at least two pairs of DQA and DQB genes, termed DQA1, DQB1 and DQA2, DQB2 and both sets are polymorphic and both seem to be functional. In a number of haplotypes, the DQA1 locus appears to be absent (DQA1-null) and is replaced by a second locus termed DQA2-like. Here, we identify families of alleles based on sequence similarity and phylogenetic clustering which correspond to each of the DQA and DQB genes identified in previous genomic and transcript analyses of homozygous animals. Using such criteria to cluster sequences, we have named 82 full-length and partial cDNA transcripts derived from domestic sheep (Ovis aries) which correspond to alleles at the Ovar-DQA1, DQA2, DQA2-like, DQB1, DQB2 and DQB2-like genes and provide associated sequence resources available to the research community through the IPD-MHC Database. This sets the basis for naming and annotation of DQ genes within the ovine MHC and may be used as a template for DQ genes in other ruminant species which will ultimately support research in livestock infectious disease.

Functional variation at an expressed MHC class IIß locus associates with Ranavirus infection intensity in larval anuran populations.

Infectious diseases are causing catastrophic losses to global biodiversity. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian pathogen susceptibility. MHC genes encode two classes of polymorphic cell-surface molecules that can recognize and bind to diverse pathogen peptides. While MHC class I genes are the classic mediators of viral-acquired immunity, larval amphibians do not express them. Consequently, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Rana sylvatica), which are highly susceptible to Ranavirus, across 17 ponds and 2 years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIß locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if supertypes or alleles influenced host responses to Ranavirus. Among 381 sampled individuals, 26% were infected with Ranavirus. We recovered 20 unique MHC class IIß alleles that fell into two deeply diverged clades and seven supertypes. MHC genotypes were associated with Ranavirus infection intensity, but not prevalence. Specifically, MHC heterozygotes and supertype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and other supertypes. We conclude that MHC class IIß functional genetic variation is an important component of Ranavirus susceptibility. Identifying immunogenetic signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

Identification of the major histocompatibility complex class-II DM and DO alleles in a cohort of northern pig-tailed macaques (Macaca leonina).

The northern pig-tailed macaque (Macaca leonina) has been considered as an independent species from the pig-tailed macaque group. We have previously reported that this species macaque has the potential to be a useful animal model in HIV/AIDS pathogenesis and vaccine studies due to its susceptibility to HIV-1. To develop this animal into a potential HIV/AIDS model, we have studied the classical MHC genes of this animal. In this study, the non-classical MHC genes Malo-DM and Malo-DO alleles were first characterized by sequencing and cloning in 12 unrelated northern pig-tailed macaques. A total of 20 full-length sequences identified include 4 Malo-DMA, 5 Malo-DMB, 7 Malo-DOA, and 4 Malo-DOB alleles. Most of these allele sequences were shared between northern pig-tailed macaque and other macaque species in exon 2. The full-length MHC-DM and MHC-DO sequences provide more comprehensive analysis of immunogenetics of northern pig-tailed macaques and increase the value of the macaques in further biomedical studies.

Transcriptome analyses of immune tissues from three Japanese frogs (genus Rana ) reveals their utility in characterizing major histocompatibility complex class II.

BACKGROUND: In Japan and East Asia, endemic frogs appear to be tolerant or not susceptible to chytridiomycosis, a deadly amphibian disease caused by the chytrid fungus Batrachochytridium dendrobatidis (Bd). Japanese frogs may have evolved mechanisms of immune resistance to pathogens such as Bd. This study characterizes immune genes expressed in various tissues of healthy Japanese Rana frogs. RESULTS: We generated transcriptome data sets of skin, spleen and blood from three adult Japanese Ranidae frogs (Japanese brown frog Rana japonica, the montane brown frog Rana ornativentris, and Tago's brown frog Rana tagoi tagoi) as well as whole body of R. japonica and R. ornativentris tadpoles. From this, we identified tissue- and stage-specific differentially expressed genes; in particular, the spleen was most enriched for immune-related genes. A specific immune gene, major histocompatibility complex class IIB (MHC-IIB), was further characterized due to its role in pathogen recognition. We identified a total of 33 MHC-IIB variants from the three focal species (n = 7 individuals each), which displayed evolutionary signatures related to increased MHC variation, including balancing selection. Our supertyping analyses of MHC-IIB variants from Japanese frogs and previously studied frog species identified potential physiochemical properties of MHC-II that may be important for recognizing and binding chytrid-related antigens. CONCLUSIONS: This is one of the first studies to generate transcriptomic resources for Japanese frogs, and contributes to further understanding the immunogenetic factors associated with resistance to infectious diseases in amphibians such as chytridiomycosis. Notably, MHC-IIB supertyping analyses identified unique functional properties of specific MHC-IIB alleles that may partially contribute to Bd resistance, and such properties provide a springboard for future experimental validation.

Genetic diversity comparison of the DQA gene in European rabbit (Oryctolagus cuniculus) populations.

The European rabbit (Oryctolagus cuniculus) natural populations within the species native region, the Iberian Peninsula, are considered a reservoir of genetic diversity. Indeed, the Iberia was a Pleistocene refuge to the species and currently two subspecies are found in the peninsula (Oryctolagus cuniculus cuniculus and Oryctolagus cuniculus algirus). The genes of the major histocompatibility complex (MHC) have been substantially studied in wild populations due to their exceptional variability, believed to be pathogen driven. They play an important function as part of the adaptive immune system affecting the individual fitness and population viability. In this study, the MHC variability was assessed by analysing the exon 2 of the DQA gene in several European rabbit populations from Portugal, Spain and France and in domestic breeds. Twenty-eight DQA alleles were detected, among which 18 are described for the first time. The Iberian rabbit populations are well differentiated from the French population and domestic breeds. The Iberian populations retained the higher allelic diversity with the domestic breeds harbouring the lowest; in contrast, the DQA nucleotide diversity was higher in the French population. Signatures of positive selection were detected in four codons which are putative peptide-binding sites and have been previously detected in other mammals. The evolutionary relationships showed instances of trans-species polymorphism. Overall, our results suggest that the DQA in European rabbits is evolving under selection and genetic drift.

Copy number variation and genetic diversity of MHC Class IIb alleles in an alien population of Xenopus laevis.

Xenopus laevis (the African clawed frog), which originated through hybridisation and whole genome duplication, has been used as a model for genetics and development for many years, but surprisingly little is known about immune gene variation in natural populations. The purpose of this study was to use an isolated population of X. laevis that was introduced to Wales, UK in the past 50 years to investigate how variation at the MHC compares to that at other loci, following a severe population bottleneck. Among 18 individuals, we found nine alleles based on exon 2 sequences of the Class IIb region (which includes the peptide binding region). Individuals carried from one to three of the loci identified from previous laboratory studies. Genetic variation was an order of magnitude higher at the MHC compared with three single-copy nuclear genes, but all loci showed high levels of heterozygosity and nucleotide diversity and there was not an excess of homozygosity or decrease in diversity over time that would suggest extensive inbreeding in the introduced population. Tajima's D was positive for all loci, which is consistent with a bottleneck. Moreover, comparison with published sequences identified the source of the introduced population as the Western Cape region of South Africa, where most commercial suppliers have obtained their stocks. These factors suggest that despite founding by potentially already inbred individuals, the alien population in Wales has maintained substantial genetic variation at both adaptively important and neutral genes.

The diversity of major histocompatibility complex class II DRB1 gene in sheep breeds from Xinjiang, China.

Exon 2 of the ovine leukocyte antigen OLA-DRB1 locus was examined in sheep from the Xinjiang Karakul Ram and Bashibai populations, and three generations of hybrids were derived from a cross between Bashibai and Altai Argali wild sheep. This identified 12 novel alleles and 30 previously reported alleles. A neighbor-joining tree of the amino acid sequences of these 42 alleles revealed allelic clusters shared across the study populations. There were significant differences in allelic frequency between Karakul Ram and Bashibai sheep. DRB1*K18cC was the most frequent allele in Kararul Ram with a frequency of 21.2%, while DRB1*2F10c8 (13.2%) and DRB1*0803 (13.2%) were the most frequent alleles found in Bashibai sheep; the alleles DRB1*2F16c2, DRB1*1601, and DRB1*0803 occurred most frequently in F1, F2, and F3 populations, with frequencies of 17.6%, 14.3%, and 20%, respectively. Although many alleles were shared by Bashibai and hybrid sheep, some alleles differed between them, especially in the F1 generation of the Bashibai × Altai Argali cross. The hybrid-specific alleles indicated the introgression of Altai Argali alleles into hybrid flocks. A population tree based on the OLA-DRB1 allelic frequency in each population indicated that the Bashibai sheep and three hybrid populations were similar, with Karakul Ram being genetically distinct.

Assessment of biodiversity in Chilean cattle using the distribution of major histocompatibility complex class II BoLA-DRB3 allele.

Bovine leukocyte antigens (BoLAs) are used extensively as markers for bovine disease and immunological traits. In this study, we estimated BoLA-DRB3 allele frequencies using 888 cattle from 10 groups, including seven cattle breeds and three crossbreeds: 99 Red Angus, 100 Black Angus, 81 Chilean Wagyu, 49 Hereford, 95 Hereford × Angus, 71 Hereford × Jersey, 20 Hereford × Overo Colorado, 113 Holstein, 136 Overo Colorado, and 124 Overo Negro cattle. Forty-six BoLA-DRB3 alleles were identified, and each group had between 12 and 29 different BoLA-DRB3 alleles. Overo Negro had the highest number of alleles (29); this breed is considered in Chile to be an 'Old type' European Holstein Friesian descendant. By contrast, we detected 21 alleles in Holstein cattle, which are considered to be a 'Present type' Holstein Friesian cattle. Chilean cattle groups and four Japanese breeds were compared by neighbor-joining trees and a principal component analysis (PCA). The phylogenetic tree showed that Red Angus and Black Angus cattle were in the same clade, crossbreeds were closely related to their parent breeds, and Holstein cattle from Chile were closely related to Holstein cattle in Japan. Overall, the tree provided a thorough description of breed history. It also showed that the Overo Negro breed was closely related to the Holstein breed, consistent with historical data indicating that Overo Negro is an 'Old type' Holstein Friesian cattle. This allelic information will be important for investigating the relationship between major histocompatibility complex (MHC) and disease.

Sequence variations of the locus-specific 5' untranslated regions of SLA class I genes and the development of a comprehensive genomic DNA-based high-resolution typing method for SLA-2.

The genetic diversity of the major histocompatibility complex (MHC) class I molecules of pigs has not been well characterized. Therefore, the influence of MHC genetic diversity on the immune-related traits of pigs, including disease resistance and other MHC-dependent traits, is not well understood. Here, we attempted to develop an efficient method for systemic analysis of the polymorphisms in the epitope-binding region of swine leukocyte antigens (SLA) class I genes. We performed a comparative analysis of the last 92 bp of the 5' untranslated region (UTR) to the beginning of exon 4 of six SLA classical class I-related genes, SLA-1, -2, -3, -4, -5, and -9, from 36 different sequences. Based on this information, we developed a genomic polymerase chain reaction (PCR) and direct sequencing-based comprehensive typing method for SLA-2. We successfully typed SLA-2 from 400 pigs and 8 cell lines, consisting of 9 different pig breeds, and identified 49 SLA-2 alleles, including 31 previously reported alleles and 18 new alleles. We observed differences in the composition of SLA-2 alleles among different breeds. Our method can be used to study other SLA class I loci and to deepen our knowledge of MHC class I genes in pigs.

Eco-immunology of fish invasions: the role of MHC variation.

The relationship between invaders and the pathogens encountered in their new environment can have a large effect on invasion success. Invaders can become free from their natural pathogens and reallocate costly immune resources to growth and reproduction, thereby increasing invasion success. Release from enemies and relaxation of selective pressures could render newly founded populations more variable at immune-related genes, such as the major histocompatibility complex (MHC), particularly when they have different origins. Using rainbow and brown trout, two of the world's most successful fish invaders, we tested the general hypothesis that invaders should display high intrapopulation immunogenetic diversity and interpopulation divergence, due to the interplay between genetic drift and successive waves of genetically divergent introductions. We analysed genetic diversity and signatures of selection at the MHC class II ß immune-related locus. In both species, MHC diversity (allelic richness and heterozygosity) for southern hemisphere populations was similar to values reported for populations at their native range. However, MHC functional diversity was limited, and population immunogenetic structuring weaker than that observed using neutral markers. Depleted MHC functional diversity could reflect a decrease in immune response, immune-related assortative mating or selection for resistance to newly encountered parasites. Given that the role of MHC diversity in the survival of these populations remains unclear, depleted functional diversity of invasive salmonids could compromise their long-term persistence. A better understanding of the eco-immunology of invaders may help in managing and preventing the impact of biological invasions, a major cause of loss of biodiversity worldwide.

Nucleotide sequence and polymorphism of MHC class II DQB exon 2 alleles in Chinese yakow (Bos grunniens × Bos taurus).

The objective of this study is to describe the gene characteristics of the bovine lymphocyte antigen (BoLA)-DQB exon 2 locus in Chinese yakow (Bos grunniens × Bos taurus) and to compare it with previously reported patterns in other bovidae species to investigate genetic factors for disease resistance. The exon 2 of the MHC class II DQB gene was cloned and sequenced. It was revealed by sequence analyses that there are 36 DQB exon 2 alleles among 44 Chinese yakow. These alleles exhibited a high degree of nucleotide and amino acid polymorphism with most amino acid variations occurring at positions forming the peptide-binding sites (PBS). The DQB loci were analysed for patterns of synonymous (dS ) and nonsynonymous (dN ) substitution. The Chinese yakow was observed to be under strong positive selection in the DQB exon 2 peptide-binding sites (dN = 0.147, P < 0.01). It appears that this variability among Chinese yakow confers the ability to mount immune responses to a wide variety of peptides or pathogens.

Consensus classification of human leukocyte antigen class II proteins.

Class II human leukocyte antigens (HLA II) are proteins involved in the human immunological adaptive response by binding and exposing some pre-processed, non-self peptides in the extracellular domain in order to make them recognizable by the CD4+ T lymphocytes. However, the understanding of HLA-peptide binding interaction is a crucial step for designing a peptide-based vaccine because the high rate of polymorphisms in HLA class II molecules creates a big challenge, even though the HLA II proteins can be grouped into supertypes, where members of different class bind a similar pool of peptides. Hence, first we performed the supertype classification of 27 HLA II proteins using their binding affinities and structural-based linear motifs to create a stable group of supertypes. For this purpose, a well-known clustering method was used, and then, a consensus was built to find the stable groups and to show the functional and structural correlation of HLA II proteins. Thus, the overlap of the binding events was measured, confirming a large promiscuity within the HLA II-peptide interactions. Moreover, a very low rate of locus-specific binding events was observed for the HLA-DP genetic locus, suggesting a different binding selectivity of these proteins with respect to HLA-DR and HLA-DQ proteins. Secondly, a predictor based on a support vector machine (SVM) classifier was designed to recognize HLA II-binding peptides. The efficiency of prediction was estimated using precision, recall (sensitivity), specificity, accuracy, F-measure, and area under the ROC curve values of random subsampled dataset in comparison with other supervised classifiers. Also the leave-one-out cross-validation was performed to establish the efficiency of the predictor. The availability of HLA II-peptide interaction dataset, HLA II-binding motifs, high-quality amino acid indices, peptide dataset for SVM training, and MATLAB code of the predictor is available at http://sysbio.icm.edu.pl/HLA .

NetMHCIIpan-3.0, a common pan-specific MHC class II prediction method including all three human MHC class II isotypes, HLA-DR, HLA-DP and HLA-DQ.

Major histocompatibility complex class II (MHCII) molecules play an important role in cell-mediated immunity. They present specific peptides derived from endosomal proteins for recognition by T helper cells. The identification of peptides that bind to MHCII molecules is therefore of great importance for understanding the nature of immune responses and identifying T cell epitopes for the design of new vaccines and immunotherapies. Given the large number of MHC variants, and the costly experimental procedures needed to evaluate individual peptide-MHC interactions, computational predictions have become particularly attractive as first-line methods in epitope discovery. However, only a few so-called pan-specific prediction methods capable of predicting binding to any MHC molecule with known protein sequence are currently available, and all of them are limited to HLA-DR. Here, we present the first pan-specific method capable of predicting peptide binding to any HLA class II molecule with a defined protein sequence. The method employs a strategy common for HLA-DR, HLA-DP and HLA-DQ molecules to define the peptide-binding MHC environment in terms of a pseudo sequence. This strategy allows the inclusion of new molecules even from other species. The method was evaluated in several benchmarks and demonstrates a significant improvement over molecule-specific methods as well as the ability to predict peptide binding of previously uncharacterised MHCII molecules. To the best of our knowledge, the NetMHCIIpan-3.0 method is the first pan-specific predictor covering all HLA class II molecules with known sequences including HLA-DR, HLA-DP, and HLA-DQ. The NetMHCpan-3.0 method is available at http://www.cbs.dtu.dk/services/NetMHCIIpan-3.0 .

Identification and characterization of the major histocompatibility complex class II DQB (MhcMath-DQB1) alleles in Tibetan macaques (Macaca thibetana).

Tibetan macaque (Macaca thibetana), an endangered primate species endemic to China, have been used as experimental animal model for various human diseases. Major histocompatibility complex (MHC) genes play a crucial role in the susceptibility and/or resistance to many human diseases, but little is known about Tibetan macaques. To gain an insight into the MHC background and to facilitate the experimental use of Tibetan macaques, the second exon of Mhc-DQB1 gene was sequenced in a cohort of wild Tibetan macaques living in the Sichuan province of China. A total of 23 MhcMath-DQB1 alleles were identified for the first time, illustrating a marked allelic polymorphism at the DQB1 locus for these macaques. Most of the sequences (74%) observed in this study belong to DQB1*06 (9 alleles) and DQB1*18 (8 alleles) lineages, and the rest (26%) belong to DQB1*15 (3 alleles) and DQB1*17 (3 alleles) lineages. The most frequent alleles detected among these macaques were MhcMath-DQB1*15:02:02 (17.9%), followed by Math-DQB1*06:06, 17:03 and 18:01, which were detected in 9 (16.1%) of the monkeys, respectively. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the peptide-binding region, suggesting balancing selection for maintaining polymorphisms at the MHC class II DQB1 locus. Phylogenetic analyses confirms the trans-species model of evolution of the Mhc-DQB1 genes in non-human primates, and in particular, the extensive allele sharing is observed between Tibetan and other macaque species.

Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (Rhinopithecus roxellana).

BACKGROUND: Small, isolated populations often experience loss of genetic variation due to random genetic drift. Unlike neutral or nearly neutral markers (such as mitochondrial genes or microsatellites), major histocompatibility complex (MHC) genes in these populations may retain high levels of polymorphism due to balancing selection. The relative roles of balancing selection and genetic drift in either small isolated or bottlenecked populations remain controversial. In this study, we examined the mechanisms maintaining polymorphisms of MHC genes in small isolated populations of the endangered golden snub-nosed monkey (Rhinopithecus roxellana) by comparing genetic variation found in MHC and microsatellite loci. There are few studies of this kind conducted on highly endangered primate species. RESULTS: Two MHC genes were sequenced and sixteen microsatellite loci were genotyped from samples representing three isolated populations. We isolated nine DQA1 alleles and sixteen DQB1 alleles and validated expression of the alleles. Lowest genetic variation for both MHC and microsatellites was found in the Shennongjia (SNJ) population. Historical balancing selection was revealed at both the DQA1 and DQB1 loci, as revealed by excess non-synonymous substitutions at antigen binding sites (ABS) and maximum-likelihood-based random-site models. Patterns of microsatellite variation revealed population structure. FST outlier analysis showed that population differentiation at the two MHC loci was similar to the microsatellite loci. CONCLUSIONS: MHC genes and microsatellite loci showed the same allelic richness pattern with the lowest genetic variation occurring in SNJ, suggesting that genetic drift played a prominent role in these isolated populations. As MHC genes are subject to selective pressures, the maintenance of genetic variation is of particular interest in small, long-isolated populations. The results of this study may contribute to captive breeding and translocation programs for endangered species.

Organization of Mhc class II A and B genes in the tilapiine fish Oreochromis.

Perch-like fishes of the family Cichlidae are models for the study of speciation. An important tool in these studies is the major histocompatibility complex (Mhc) and its organization. The present study takes the first step toward the elucidation of the Mhc class II gene organization in the tilapiine fish Oreochromis niloticus (Orni). Using class II A- and class II B-specific probes, Mhc-bearing clones were identified and isolated from a bacterial artificial chromosome (BAC) library. The analysis of these clones by a combination of molecular, genetic-mapping, and phylogenetic methods led to the identification of nine class II A and 15 class II B loci. Genes at these loci constitute two families, which we designate as class IIa and class IIb families. Each of the families contains A and B loci. Some genes in both families are expressed and functional. The two families differ in their chromosomal location (they are unlinked) and their mode of evolution. The class IIa family genes are conserved across different teleost taxonomical orders, whereas the class IIb family genes are apparently products of multiple, more recent, rounds of gene duplications. The rounds established at least five monophyletic groups of genes. The founding unit of each monophyletic group might have been a pair of class II A and B loci.

MHC class II genes in the European badger (Meles meles): characterization, patterns of variation, and transcription analysis.

The major histocompatibility complex (MHC) comprises many genes, some of which are polymorphic with numerous alleles. Sequence variation among alleles is most pronounced in exon 2 of the class II genes, which encodes the α1 and ß1 domains that form the antigen-binding site (ABS) for the presentation of peptides. The MHC thus plays an important role in pathogen defense. European badgers (Meles meles) are a good species in which to study the MHC, as they harbor a variety of pathogens. We present the first characterization of MHC class II genes, isolated from genomic DNA (gDNA) and complementary DNA (cDNA), in the European badger. Examination of seven individuals revealed four DRB, two DQB, two DQA, and two DRA putatively functional gDNA sequences. All of these sequences, except DRA, exhibited high variability in exon 2; DRB had the highest variability. The ABS codons demonstrated high variability, due potentially to balancing selection, while non-ABS codons had lower variability. Positively selected sites were detected in DRB and DQA. Phylogenetic analysis demonstrated trans-species polymorphism of class II genes. Comparison with cDNA from whole blood revealed that only DRB had a transcription pattern reflecting the alleles that were present in the gDNA, while the other three genes had disparities between gDNA and cDNA. Only one sequence was transcribed, even though two gDNA sequences were present, from each of both DQB and DRA. Our characterization of badger MHC sequences forms a basis for further studies of MHC variability, mate choice, and pathogen resistance in this, and other, species.

454 sequencing reveals extreme complexity of the class II Major Histocompatibility Complex in the collared flycatcher.

BACKGROUND: Because of their functional significance, the Major Histocompatibility Complex (MHC) class I and II genes have been the subject of continuous interest in the fields of ecology, evolution and conservation. In some vertebrate groups MHC consists of multiple loci with similar alleles; therefore, the multiple loci must be genotyped simultaneously. In such complex systems, understanding of the evolutionary patterns and their causes has been limited due to challenges posed by genotyping. RESULTS: Here we used 454 amplicon sequencing to characterize MHC class IIB exon 2 variation in the collared flycatcher, an important organism in evolutionary and immuno-ecological studies. On the basis of over 152,000 sequencing reads we identified 194 putative alleles in 237 individuals. We found an extreme complexity of the MHC class IIB in the collared flycatchers, with our estimates pointing to the presence of at least nine expressed loci and a large, though difficult to estimate precisely, number of pseudogene loci. Many similar alleles occurred in the pseudogenes indicating either a series of recent duplications or extensive concerted evolution. The expressed alleles showed unambiguous signals of historical selection and the occurrence of apparent interlocus exchange of alleles. Placing the collared flycatcher's MHC sequences in the context of passerine diversity revealed transspecific MHC class II evolution within the Muscicapidae family. CONCLUSIONS: 454 amplicon sequencing is an effective tool for advancing our understanding of the MHC class II structure and evolutionary patterns in Passeriformes. We found a highly dynamic pattern of evolution of MHC class IIB genes with strong signals of selection and pronounced sequence divergence in expressed genes, in contrast to the apparent sequence homogenization in pseudogenes. We show that next generation sequencing offers a universal, affordable method for the characterization and, in perspective, genotyping of MHC systems of virtually any complexity.