Sequence variability analysis of human class I and class II MHC molecules: functional and structural correlates of amino acid polymorphisms

OptMAVEn--a new framework for the de novo design of antibody variable region models targeting specific antigen epitopes

Antibody-based therapeutics provides novel and efficacious treatments for a number of diseases. Traditional experimental approaches for designing therapeutic antibodies rely on raising antibodies against a target antigen in an immunized animal or directed evolution of antibodies with low affinity for the desired antigen. However, these methods remain time consuming, cannot target a specific epitope and do not lead to broad design principles informing other studies. Computational design methods can overcome some of these limitations by using biophysics models to rationally select antibody parts that maximize affinity for a target antigen epitope. This has been addressed to some extend by OptCDR for the design of complementary determining regions. Here, we extend this earlier contribution by addressing the de novo design of a model of the entire antibody variable region against a given antigen epitope while safeguarding for immunogenicity (Optimal Method for Antibody Variable region Engineering, OptMAVEn). OptMAVEn simulates in silico the in vivo steps of antibody generation and evolution, and is capable of capturing the critical structural features responsible for affinity maturation of antibodies. In addition, a humanization procedure was developed and incorporated into OptMAVEn to minimize the potential immunogenicity of the designed antibody models. As case studies, OptMAVEn was applied to design models of neutralizing antibodies targeting influenza hemagglutinin and HIV gp120. For both HA and gp120, novel computational antibody models with numerous interactions with their target epitopes were generated. The observed rates of mutations and types of amino acid changes during in silico affinity maturation are consistent with what has been observed during in vivo affinity maturation. The results demonstrate that OptMAVEn can efficiently generate diverse computational antibody models with both optimized binding affinity to antigens and reduced immunogenicity.

Nonsynonymous substitution rate heterogeneity in the peptide-binding region among different HLA-DRB1 lineages in humans

An extraordinary diversity of amino acid sequences in the peptide-binding region (PBR) of human leukocyte antigen [HLA; human major histocompatibility complex (MHC)] molecules has been maintained by balancing selection. The process of accumulation of amino acid diversity in the PBR for six HLA genes (HLA-A, B, C, DRB1, DQB1, and DPB1) shows that the number of amino acid substitutions in the PBR among alleles does not linearly correlate with the divergence time of alleles at the six HLA loci. At these loci, some pairs of alleles show significantly less nonsynonymous substitutions at the PBR than expected from the divergence time. The same phenomenon was observed not only in the HLA but also in the rat MHC. To identify the cause for this, DRB1 sequences, a representative case of a typical nonlinear pattern of substitutions, were examined. When the amino acid substitutions in the PBR were placed with maximum parsimony on a maximum likelihood tree based on the non-PBR substitutions, heterogeneous rates of nonsynonymous substitutions in the PBR were observed on several branches. A computer simulation supported the hypothesis that allelic pairs with low PBR substitution rates were responsible for the stagnation of accumulation of PBR nonsynonymous substitutions. From these observations, we conclude that the nonsynonymous substitution rate at the PBR sites is not constant among the allelic lineages. The deceleration of the rate may be caused by the coexistence of certain pathogens for a substantially long time during HLA evolution.

Prediction and characterization of helper T-cell epitopes from pneumococcal surface adhesin A

Pneumococcal surface adhesin A (PsaA) is a multifunctional lipoprotein known to bind nasopharyngeal epithelial cells, and is significantly involved in bacterial adherence and virulence. Identification of PsaA peptides that optimally bind human leucocyte antigen (HLA) and elicit a potent immune response would be of great importance to vaccine development. However, this is hindered by the multitude of HLA polymorphisms in humans. To identify the conserved immunodominant epitopes, we used an experimental dataset of 28 PsaA synthetic peptides and in silico methods to predict specific peptide-binding to HLA and murine MHC class II molecules. We also characterized spleen and cervical lymph node (CLN) -derived T helper (Th) lymphocyte cytokine responses to these peptides after Streptococcus pneumoniae strain EF3030 challenge in mice. Individual, yet overlapping, peptides 15 amino acids in length revealed residues of PsaA that consistently caused the highest interferon-γ, interleukin-2 (IL-2), IL-5 and IL-17 responses and proliferation as well as moderate IL-10 and IL-4 responses by ex vivo re-stimulated splenic and CLN CD4⁺ T cells isolated from S. pneumoniae strain EF3030-challenged F1 (B6 × BALB/c) mice. In silico analysis revealed that peptides from PsaA may interact with a broad range of HLA-DP, -DQ and -DR alleles, due in part to regions lacking β-turns and asparagine endopeptidase sites. These data suggest that Th cell peptides (7, 19, 20, 22, 23 and 24) screened for secondary structures and MHC class II peptide-binding affinities can elicit T helper cytokine and proliferative responses to PsaA peptides.

Study of MHC class II region polymorphism in the Filipino cynomolgus macaque population

The cynomolgus macaque (Macaca fascicularis) is currently used as an animal model in various fields of immunology especially in the development of innovative vaccines for the prevention and treatment of infectious diseases. The polymorphism of the major histocompatibility complex (MHC) influences the development of adaptive immune responses and it is crucial to characterize the polymorphism of cynomolgus MHC genes. We present here a systematic study of the MHC class II haplotypes in the Filipino macaque population. By the study of a large sample of Filipino animals (N = 353), we have characterized 18 MHC class II haplotypes by means of genotyping seven microsatellites. The animals were DRB genotyped by means of PCR-SSO or DGGE-sequencing on genomic amplified fragments. We cloned and sequenced the complementary DNA (cDNA) of DQA, DQB, DPA, and DPB genes of 117 animals. Combining the microsatellite genotyping and cDNA characterized in the 117 animals, we defined genetic association between the cDNA and the microsatellites and characterized 18 MHC class II haplotypes. For 104 animals out of the 353 studied, the presence of a recombinant haplotype was highly probable. Thirty-four percent of recombination was located in 256 kb segment between D6S2876 and D6S2747 microsatellites, a region encompassing several hot spots of recombination in the human MHC.

Mutant OmpF porins of Yersinia pseudotuberculosis with deletions of external loops: structure-functional and immunochemical properties

Recombinant mutant OmpF porins from Yersinia pseudotuberculosis outer membrane were obtained using site-directed mutagenesis. Here we used four OmpF mutants where single extracellular loops L1, L4, L6, and L8 were deleted one at a time. The proteins were expressed in Escherichia coli at levels comparable to full-sized recombinant OmpF porin and isolated from the inclusion bodies. Purified trimers of the mutant porins were obtained after dialysis and consequent ion-exchange chromatography. Changes in molecular and spatial structure of the mutants obtained were studied using SDS-PAGE and optical spectroscopy (circular dichroism and intrinsic protein fluorescence). Secondary and tertiary structure of the mutant proteins was found to have some features in comparison with that of the full-sized recombinant OmpF. As shown by bilayer lipid membrane technique, the pore-forming activity of purified mutant porins was identical to OmpF porin isolated from the bacterial outer membrane. Lacking of the external loops mentioned above influenced significantly upon the antigenic structure of the porin as demonstrated using ELISA.

Quantifying Significance of MHC II Residues

The major histocompatibility complex (MHC), a cell-surface protein mediating immune recognition, plays important roles in the immune response system of all higher vertebrates. MHC molecules are highly polymorphic and they are grouped into serotypes according to the specificity of the response. It is a common belief that a protein sequence determines its three dimensional structure and function. Hence, the protein sequence determines the serotype. Residues play different levels of importance. In this paper, we quantify the residue significance with the available serotype information. Knowing the significance of the residues will deepen our understanding of the MHC molecules and yield us a concise representation of the molecules. In this paper we propose a linear programming-based approach to find significant residue positions as well as quantifying their significance in MHC II DR molecules. Among all the residues in MHC II DR molecules, 18 positions are of particular significance, which is consistent with the literature on MHC binding sites, and succinct pseudo-sequences appear to be adequate to capture the whole sequence features. When the result is used for classification of MHC molecules with serotype assigned by WHO, a 98.4 percent prediction performance is achieved. The methods have been implemented in java (http://code.google.com/p/quassi/).

A new HLA allele, HLA-B*08:108, described in two unrelated Spanish individuals

HLA-B*08:108 shows one nucleotide difference regarding B*08:01:01 at codon 109 (CTC>TTC, L109>F109).

Selection of conserved epitopes from hepatitis C virus for pan-populational stimulation of T-cell responses

The hepatitis C virus (HCV) is able to persist as a chronic infection, which can lead to cirrhosis and liver cancer. There is evidence that clearance of HCV is linked to strong responses by CD8 cytotoxic T lymphocytes (CTLs), suggesting that eliciting CTL responses against HCV through an epitope-based vaccine could prove an effective means of immunization. However, HCV genomic plasticity as well as the polymorphisms of HLA I molecules restricting CD8 T-cell responses challenges the selection of epitopes for a widely protective vaccine. Here, we devised an approach to overcome these limitations. From available databases, we first collected a set of 245 HCV-specific CD8 T-cell epitopes, all known to be targeted in the course of a natural infection in humans. After a sequence variability analysis, we next identified 17 highly invariant epitopes. Subsequently, we predicted the epitope HLA I binding profiles that determine their potential presentation and recognition. Finally, using the relevant HLA I-genetic frequencies, we identified various epitope subsets encompassing 6 conserved HCV-specific CTL epitopes each predicted to elicit an effective T-cell response in any individual regardless of their HLA I background. We implemented this epitope selection approach for free public use at the EPISOPT web server.

NMR spectroscopy reveals unexpected structural variation at the protein-protein interface in MHC class I molecules

β2-Microglobulin (β2m) is a small, monomorphic protein non-covalently bound to the heavy chain (HC) in polymorphic major histocompatibility complex (MHC) class I molecules. Given the high evolutionary conservation of structural features of β2m in various MHC molecules as shown by X-ray crystallography, β2m is often considered as a mere scaffolding protein. Using nuclear magnetic resonance (NMR) spectroscopy, we investigate here whether β2m residues at the interface to the HC exhibit changes depending on HC polymorphisms and the peptides bound to the complex in solution. First we show that human β2m can effectively be produced in deuterated form using high-cell-density-fermentation and we employ the NMR resonance assignments obtained for triple-labeled β2m bound to the HLA-B*27:09 HC to examine the β2m-HC interface. We then proceed to compare the resonances of β2m in two minimally distinct subtypes, HLA-B*27:09 and HLA-B*27:05, that are differentially associated with the spondyloarthropathy Ankylosing Spondylitis. Each of these subtypes is complexed with four distinct peptides for which structural information is already available. We find that only the resonances at the β2m-HC interface show a variation of their chemical shifts between the different complexes. This indicates the existence of an unexpected plasticity that enables β2m to accommodate changes that depend on HC polymorphism as well as on the bound peptide through subtle structural variations of the protein-protein interface.

Mate choice for genetic compatibility in the house mouse

In house mice, genetic compatibility is influenced by the t haplotype, a driving selfish genetic element with a recessive lethal allele, imposing fundamental costs on mate choice decisions. Here, we evaluate the cost of genetic incompatibility and its implication for mate choice in a wild house mice population. In laboratory reared mice, we detected no fertility (number of embryos) or fecundity (ability to conceive) costs of the t, and yet we found a high cost of genetic incompatibility: heterozygote crosses produced 40% smaller birth litter sizes because of prenatal mortality. Surprisingly, transmission of t in crosses using +/t males was influenced by female genotype, consistent with postcopulatory female choice for + sperm in +/t females. Analysis of paternity patterns in a wild population of house mice showed that +/t females were more likely than +/+ females to have offspring sired by +/+ males, and unlike +/+ females, paternity of their offspring was not influenced by +/t male frequency, further supporting mate choice for genetic compatibility. As the major histocompatibility complex (MHC) is physically linked to the t, we investigated whether females could potentially use variation at the MHC to identify male genotype at the sperm or individual level. A unique MHC haplotype is linked to the t haplotype. This MHC haplotype could allow the recognition of t and enable pre- and postcopulatory mate choice for genetic compatibility. Alternatively, the MHC itself could be the target of mate choice for genetic compatibility. We predict that mate choice for genetic compatibility will be difficult to find in many systems, as only weak fertilization biases were found despite an exceptionally high cost of genetic incompatibility.

PREDIVAC: CD4+ T-cell epitope prediction for vaccine design that covers 95% of HLA class II DR protein diversity

Background

CD4+ T-cell epitopes play a crucial role in eliciting vigorous protective immune responses during peptide (epitope)-based vaccination. The prediction of these epitopes focuses on the peptide binding process by MHC class II proteins. The ability to account for MHC class II polymorphism is critical for epitope-based vaccine design tools, as different allelic variants can have different peptide repertoires. In addition, the specificity of CD4+ T-cells is often directed to a very limited set of immunodominant peptides in pathogen proteins. The ability to predict what epitopes are most likely to dominate an immune response remains a challenge.

Results

We developed the computational tool Predivac to predict CD4+ T-cell epitopes. Predivac can make predictions for 95% of all MHC class II protein variants (allotypes), a substantial advance over other available methods. Predivac bases its prediction on the concept of specificity-determining residues. The performance of the method was assessed both for high-affinity HLA class II peptide binding and CD4+ T-cell epitope prediction. In terms of epitope prediction, Predivac outperformed three available pan-specific approaches (delivering the highest specificity). A central finding was the high accuracy delivered by the method in the identification of immunodominant and promiscuous CD4+ T-cell epitopes, which play an essential role in epitope-based vaccine design.

Conclusions

The comprehensive HLA class II allele coverage along with the high specificity in identifying immunodominant CD4+ T-cell epitopes makes Predivac a valuable tool to aid epitope-based vaccine design in the context of a genetically heterogeneous human population.The tool is available at: http://predivac.biosci.uq.edu.au/.

A rationally engineered anti-HIV peptide fusion inhibitor with greatly reduced immunogenicity

Peptides derived from the C-terminal heptad repeat 2 (HR2) region of the HIV-1 gp41 envelope glycoprotein, so-called C peptides, are very efficient HIV-1 fusion inhibitors. We previously developed innovative gene therapeutic approaches aiming at the direct in vivo production of C peptides from genetically modified host cells and found that T cells expressing membrane-anchored or secreted C peptides are protected from HIV-1 infection. However, an unwanted immune response against such antiviral peptides may significantly impair clinical efficacy and pose safety risks to patients. To overcome this problem, we engineered a novel C peptide, V2o, with greatly reduced immunogenicity and excellent antiviral activity. V2o is based on the chimeric C peptide C46-EHO, which is derived from the HR2 regions of HIV-2(EHO) and HIV-1(HxB2) and has broad anti-HIV and anti-simian immunodeficiency virus activity. Antibody and major histocompatibility complex class I epitopes within the C46-EHO peptide sequence were identified by in silico and in vitro analyses. Using rational design, we removed these epitopes by amino acid substitutions and thus minimized antigenicity and immunogenicity considerably. At the same time, the antiviral activity of the deimmunized peptide V2o was preserved or even enhanced compared to that of the parental C46-EHO peptide. Thus, V2o is an ideal candidate, especially for those novel therapeutic approaches for HIV infection that involve direct in vivo production of antiviral C peptides.

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.

Sequencing of two novel HLA-A*29 alleles, A*29:39 and A*29:40

A*29:39 differs from A*29:02:01 by three clustered amino acid replacements at α1 domain, T73>I73, A76>V76 and N77>D77. A*29:40 shows one nucleotide difference regarding A*29:02:01 allele, resulting in one amino acid substitution at position 154, E154>G154.

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.

Global major histocompatibility class II β (mh-IIβ)-polymorphism in Arctic charr Salvelinus alpinus

This study explored the use of the gene encoding the β subunit of the major histocompatibility (MH) receptor as a population marker in Arctic charr Salvelinus alpinus. The use of this polymorphic marker allowed differentiation of the S. alpinus lineages previously defined using mitochondrial DNA (mtDNA) but also allowed differentiation between the populations studied within those lineages. The majority of the variation observed here occurred prior to the last glaciation event. Nevertheless, all S. alpinus populations were differentiated using both MH Class II β (mh-IIβ) sequences and allelic frequencies. The fact that all the populations studied presented high rates of non-synonymous: synonymous substitutions and high levels of interpopulation variation, suggested mh-IIβ as an ideal marker to assess differentiation among S. alpinus populations in ways that may represent divergence both by genetic drift and natural adaptation to the local environment.

Insight into gene evolution within Cervidae and Bovidae through genetic variation in MHC-DQA in the black muntjac (Muntiacus crinifrons)

The critical role that the major histocompatibility complex plays in the immune recognition of parasites and pathogens makes its evolutionary dynamics exceptionally relevant to ecology, population biology, and conservation studies. The black muntjac is a rare deer endemic to a small mountainous region in eastern China. We found that this species has two DQA loci through cDNA expression and sequence variation analysis. The level of variation at both DQA loci was found to be extremely low (three alleles for DQA1 and four alleles for DQA2), possibly because of past bottlenecks and the species' relatively solitary behavior pattern. The ratio of d(N)/d(S) in the putative peptide binding region of the DQA2 locus (13.36, P = 0.012) was significantly larger than one but not that of DQA1 (0.94, P = 0.95), suggesting strong positive selection at the DQA2 but not at the DQA1 locus. This difference might reflect different sets of evolutionary selection pressures acting on the two loci. The phylogenetic tree showed that DQA1 alleles from two species of Cervidae and two of Bovidae grouped together, as did the DQA2 alleles. However, different genes from the four species were located in separate branches. These results lead us to suggest that these DQA alleles are derived from primordial DQA genes from a common ancestor and are maintained in Cervidae and Bovidae since their divergence around 25.5-27.8 million years ago.

Identification of MhcMafa-DRB alleles in a cohort of cynomolgus macaques of Vietnamese origin

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

Character of chicken polymorphic major histocompatibility complex class II alleles of 3 Chinese local breeds

To better understand the major histocompatibility complex (MHC) genetic character of domestic birds, we sequenced and analyzed chicken MHC II (B-L) genes of 3 local chicken breeds, derived from 3 separate areas in China. We amplified cDNA sequences from 105 individuals, accounting for 35 alleles. Some of the same B-LB alleles with a high frequency were found in all samples. The putative B-L α-chain had few polymorphic sites, whereas the B-L β-chain had several polymorphic sites. Most of the mutation positions were located in the B-LB β1 domain encoded by exon 2, especially in the peptide-binding region. This indicated that the highly polymorphic peptide-binding region could potentiate binding diverse antigen epitopes. The comparison of 3-D molecule structures of chicken B-L and human HLA-DR1 revealed a distinctly structural similarity, but the chicken B-L molecule had more polymorphic sites than the human HLA-DR1 molecule, which presumably might be a mechanism to compensate for responding to a wider array of pathogens due to fewer loci for chicken. Moreover, some conserved sites in human and chicken MHC class II molecules reflected their common ancestry and similar functions. These results suggest that the chicken B-L gene showed more polymorphic sites and distinctly dominant trans-breed alleles, potentially to adapt to pathogens.

The Père David's deer MHC class I genes show unexpected diversity patterns, with monomorphic classical genes but polymorphic nonclassical genes and pseudogenes

Père David's deer (Elaphurus davidianus) is a highly inbred species that arose from 11 founders but now comprises a population of about 3,000 individuals, making it interesting to investigate the adaptive variation of this species from the major histocompatibility complex (MHC) perspective. In this study, we isolated Elda-MHC class I loci using magnetic bead-based cDNA hybridization, and examined the molecular variations of these loci using single-strand conformation polymorphism (SSCP) and sequence analysis. We obtained seven MHC class I genes, which we designated F1, F12, G2, I7, AF, I8, and C1. Our analyses of stop codons, phylogenetic trees, amino acid conservation, and G+C content revealed that F1, F12, G2, and I7 were classical genes, AF was a nonclassical gene, and I8 and C1 were pseudogenes. Our subsequent molecular examinations showed that the diversity pattern in the Père David's deer was unusual. Most mammals have more polymorphic classical class I loci vs. the nonclassical and neutral genes. In contrast, the Père David's deer was found to be monomorphic at classical genes F1, F12, G2, and I7, dimorphic at the nonclassical AF gene, dimorphic at pseudogene I8, and tetramorphic at pseudogene C1. The adverse polymorphism patterns of Elda-I genes might provide evidence for selection too faster deplete MHC variation than drift in the bottlenecked populations, while the postbottleneck tetramorphism of the C1 pseudogene appears to be evidence of strong historical balancing selection.

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.

Study of cynomolgus monkey (Macaca fascicularis) Mhc DRB gene polymorphism in four populations

The cynomolgus macaque (Macaca fascicularis) is currently used as an animal model in various fields of immunology especially in the development of innovative vaccines for the prevention and treatment of infectious diseases. The polymorphism of the major histocompatibility complex (MHC) influences the development of adaptive immune responses, and it is crucial to characterize the polymorphism of cynomolgus MHC genes. Among all macaque species, the cynomolgus macaque has the most diversified geographical area encompassing continental and insular populations. By the study of a large sample of animals from the Philippines (N = 359), we have characterized 20 DRB haplotypes. The DRB genotyping was performed by denaturing gradient gel electrophoresis (DGGE) sequencing of exon 2 and was confirmed by polymerase chain reaction-sequence-specific oligonucleotide. The DRB and DRA cDNA of 126 animals were characterized by cloning and sequencing. By means of DGGE sequencing, we characterized the polymorphism of genomic DRB exon 2 in three other cynomolgus macaque population samples (Java, Vietnam, and Mauritius), and we discuss about the origin of the founders of the Mauritian and the Filipino cynomolgus macaque populations.

Psoriasis patients are enriched for genetic variants that protect against HIV-1 disease

An important paradigm in evolutionary genetics is that of a delicate balance between genetic variants that favorably boost host control of infection but which may unfavorably increase susceptibility to autoimmune disease. Here, we investigated whether patients with psoriasis, a common immune-mediated disease of the skin, are enriched for genetic variants that limit the ability of HIV-1 virus to replicate after infection. We analyzed the HLA class I and class II alleles of 1,727 Caucasian psoriasis cases and 3,581 controls and found that psoriasis patients are significantly more likely than controls to have gene variants that are protective against HIV-1 disease. This includes several HLA class I alleles associated with HIV-1 control; amino acid residues at HLA-B positions 67, 70, and 97 that mediate HIV-1 peptide binding; and the deletion polymorphism rs67384697 associated with high surface expression of HLA-C. We also found that the compound genotype KIR3DS1 plus HLA-B Bw4-80I, which respectively encode a natural killer cell activating receptor and its putative ligand, significantly increased psoriasis susceptibility. This compound genotype has also been associated with delay of progression to AIDS. Together, our results suggest that genetic variants that contribute to anti-viral immunity may predispose to the development of psoriasis.

The evolution of MHC diversity: evidence of intralocus gene conversion and recombination in a single-locus system

Gene conversion, the unidirectional exchange of genetic material between homologous sequences, is thought to strongly influence patterns of genetic diversity. The high diversity of major histocompatibility complex (MHC) genes in many species is thought to reflect a long history of gene conversion events both within and among loci. Theoretical work suggests that intra- and interlocus gene conversion leave characteristic signatures of nucleotide diversity, but empirical studies of MHC variation have rarely been able to analyze the effects of conversion events in isolation, due to the presence of multiple gene copies in most species. The potbellied seahorse (Hippocampus abdominalis), a species with a single copy of the MH class II beta-chain gene (MHIIb), provides an ideal system in which to explore predictions on the effects of intralocus gene conversion on patterns of genetic diversity. The genetic diversity of the MHIIb peptide binding region (PBR) is high in the seahorse, similar to other vertebrate species. In contrast, the remainder of the gene shows a total absence of synonymous variation and low levels of intronic sequence diversity, concentrated in 3 short repetitive regions and 1-12 SNPs per intron. The distribution of substitutions across the gene results in a patchwork pattern of shared polymorphism between otherwise divergent sequences. The pattern of nucleotide diversity observed in the seahorse MHIIb gene is congruent with theoretical expectations for intralocus gene conversion, indicating that this evolutionary mechanism has played an important role in MHC gene evolution, contributing to both the high diversity in the PBR and the low diversity outside this region. Neutral variation at this locus may be further reduced due to biases in nucleotide composition and functional constraints.

Sequencing of a new HLA-DRB1*04:98 allele in a Spanish donor

DRB1*04:98 shows four nucleotide changes regarding DRB1*04:01:01, resulting in two amino acids replacement at positions A73G and A74R.

Antigen-specific T cell phenotyping microarrays using grating coupled surface plasmon resonance imaging and surface plasmon coupled emission

The circulating population of peripheral T lymphocytes obtained from a blood sample can provide a large amount of information about an individual's medical status and history. Recent evidence indicates that the detection and functional characterization of antigen-specific T cell subsets within the circulating population may provide a diagnostic indicator of disease and has the potential to predict an individual's response to therapy. In this report, a microarray detection platform that combines grating-coupled surface plasmon resonance imaging (GCSPRI) and grating-coupled surface plasmon coupled emission (SPCE) fluorescence detection modalities were used to detect and characterize CD4(+) T cells. The microspot regions of interest (ROIs) printed on the array consisted of immobilized antibodies or peptide loaded MHC monomers (p/MHC) as T cell capture ligands mixed with additional antibodies as cytokine capture ligands covalently bound to the surface of a corrugated gold sensor chip. Using optimized parameters, an unlabeled influenza peptide reactive T cell clone could be detected at a frequency of 0.1% in a mixed T cell sample using GCSPRI. Additionally, after cell binding was quantified, differential TH1 cytokine secretion patterns from a T cell clone cultured under TH1 or TH2 inducing conditions was detected using an SPCE fluorescence based assay. Differences in the secretion patterns of 3 cytokines, characteristic of the inducing conditions, indicated that differences were a consequence of the functional status of the captured cells. A dual mode GCSPRI/SPCE assay can provide a rapid, high content T cell screening/characterization tool that is useful for diagnosing disease, evaluating vaccination efficacy, or assessing responses to immunotherapeutics.

Drug hypersensitivity and human leukocyte antigens of the major histocompatibility complex

The human leukocyte antigen (HLA) genes are the most polymorphic in the human genome and are critical in regulating specific immunity, hence their historical discovery as "immune response" genes. HLA allotypes are also implicated in unwanted immune reactions, including drug hypersensitivity syndrome, in which small therapeutic drugs interact with antigenic peptides to drive T cell responses restricted by host HLA. Abacavir, allo-purinol, and carbamazepine are three commonly used drugs that cause a T cell-mediated hypersensitivity that is HLA linked, with each drug exhibiting striking specificity for presentation by defined HLA allotypes. Recent findings have begun to unearth the mechanistic basis for these HLA associations, and here we review recent advances in the field of HLA-associated drug hypersensitivities.

Population coverage analysis of T-Cell epitopes of Neisseria meningitidis serogroup B from Iron acquisition proteins for vaccine design

Although the concept of Reverse Vaccinology was first pioneered for sepsis and meningococcal meningitidis causing bacterium, Neisseria meningitides, no broadly effective vaccine against serogroup B meningococcal disease is yet available. In the present investigation, HLA distribution analysis was undertaken to select three most promiscuous T-cell epitopes out of ten computationally validated epitopes of Iron acquisition proteins from Neisseria MC58 by using the population coverage tool of Immune Epitope Database (IEDB). These epitopes have been determined on the basis of their binding ability with maximum number of HLA alleles along with highest population coverage rate values for all the geographical areas studied. The comparative population coverage analysis of moderately immunogenic and high immunogenic peptides suggests that the former may activate T-cell response in a fairly large proportion of people in most geographical areas, thus indicating their potential for development of epitope-based vaccine.

HLA class I-associated diseases with a suspected autoimmune etiology: HLA-B27 subtypes as a model system

Although most autoimmune diseases are connected to major histocompatibility complex (MHC) class II alleles, a small number of these disorders exhibit a variable degree of association with selected MHC class I genes, like certain human HLA-A and HLA-B alleles. The basis for these associations, however, has so far remained elusive. An understanding might be obtained by comparing functional, biochemical, and biophysical properties of alleles that are minimally distinct from each other, but are nevertheless differentially associated to a given disease, like the HLA-B*27:05 and HLA-B*27:09 antigens, which differ only by a single amino acid residue (Asp116His) that is deeply buried within the binding groove. We have employed a number of approaches, including X-ray crystallography and isotope-edited infrared spectroscopy, to investigate biophysical characteristics of the two HLA-B27 subtypes complexed with up to ten different peptides. Our findings demonstrate that the binding of these peptides as well as the conformational flexibility of the subtypes is greatly influenced by interactions of the C-terminal peptide residue. In particular, a basic C-terminal peptide residue is favoured by the disease-associated subtype HLA-B*27:05, but not by HLA-B*27:09. This property appears also as the only common denominator of distinct HLA class I alleles, among them HLA-B*27:05, HLA-A*03:01 or HLA-A*11:01, that are associated with diseases suspected to have an autoimmune etiology. We postulate here that the products of these alleles, due to their unusual ability to bind with high affinity to a particular peptide set during positive T cell selection in the thymus, are involved in shaping an abnormal T cell repertoire which predisposes to the acquisition of autoimmune diseases.