Identifiying human MHC supertypes using bioinformatic methods

Classification of MHC molecules into supertypes in terms of peptide-binding specificities is an important issue, with direct implications for the development of epitope-based vaccines with wide population coverage. In view of extremely high MHC polymorphism (948 class I and 633 class II HLA alleles) the experimental solution of this task is presently impossible. In this study, we describe a bioinformatics strategy for classifying MHC molecules into supertypes using information drawn solely from three-dimensional protein structure. Two chemometric techniques–hierarchical clustering and principal component analysis–were used independently on a set of 783 HLA class I molecules to identify supertypes based on structural similarities and molecular interaction fields calculated for the peptide binding site. Eight supertypes were defined: A2, A3, A24, B7, B27, B44, C1, and C4. The two techniques gave 77% consensus, i.e., 605 HLA class I alleles were classified in the same supertype by both methods. The proposed strategy allowed “supertype fingerprints” to be identified. Thus, the A2 supertype fingerprint is Tyr9/Phe9, Arg97, and His114 or Tyr116; the A3-Tyr9/Phe9/Ser9, Ile97/Met97 and Glu114 or Asp116; the A24-Ser9 and Met97; the B7-Asn63 and Leu81; the B27-Glu63 and Leu81; for B44-Ala81; the C1-Ser77; and the C4-Asn77. action fields calculated for the peptide binding site. Eight supertypes were defined: A2, A3, A24, B7, B27, B44, C1, and C4. The two techniques gave 77% consensus, i.e., 605 HLA class I alleles were classified in the same supertype by both methods. The proposed strategy allowed “supertype fingerprints” to be identified. Thus, the A2 supertype fingerprint is Tyr9/Phe9, Arg97, and His114 or Tyr116; the A3-Tyr9/Phe9/Ser9, Ile97/Met97 and Glu114 or Asp116; the A24-Ser9 and Met97; the B7-Asn63 and Leu81; the B27-Glu63 and Leu81; for B44-Ala81; the C1-Ser77; and the C4-Asn77.

A naturally selected dimorphism within the HLA-B44 supertype alters class I structure, peptide repertoire, and T cell recognition

HLA-B*4402 and B*4403 are naturally occurring MHC class I alleles that are both found at a hi,,h frequency in all human populations, and vet they only differ by one residue on the alpha2 helix (B*4402 Aspl56-->B*4403 Leu156) CTLs discriminate between HLA-B*4402 and B*4403, and these allotypes stimulate strong mutual allogeneic responses reflecting their known barrier to hemopoeitic stem cell transplantation. Although HLA-B*4402 and B*4403 share >95% of their peptide repertoire, B*4403 presents more unique peptides than B*4402, consistent with the stronger T cell alloreactivity observed toward B*4403 compared with B*4402. Crystal structures of B*4402 and B*4403 show how the polymorphisin at position 156 is completely buried and yet alters both the peptide and the heavy chain conformation, relaxing ligand selection by B*4403 compared with B*4402. Thus, the polymorphism between HLA-B*4402 and B 4403 modifies both peptide repertoire and T cell recognition, and is reflected lit the paradoxically powerful alloreactivity that occurs across this minimal mismatch. The findings suggest that these closely related class I genes are maintained lit diverse human populations through their differential impact on the selection of peptide ligands and the T cell repertoire.

Variety matters: adaptive genetic diversity and parasite load in two mouse opossums from the Brazilian Atlantic forest

The adaptive potential of a species to a changing environment and in disease defence is primarily based on genetic variation. Immune genes, such as genes of the major histocompatibility complex (MHC), may thereby be of particular importance. In marsupials, however, there is very little knowledge about natural levels and functional importance of MHC polymorphism, despite their key role in the mammalian evolution. In a previous study, we discovered remarkable differences in the MHC class II diversity between two species of mouse opossums (Gracilinanus microtarsus, Marmosops incanus) from the Brazilian Atlantic forest, which is one of the most endangered hotspots for biodiversity conservation. Since the main forces in generating MHC diversity are assumed to be pathogens, we investigated in this study gastrointestinal parasite burden and functional associations between the individual MHC constitution and parasite load. We tested two contrasting scenarios, which might explain differences in MHC diversity between species. We predicted that a species with low MHC diversity would either be under relaxed selection pressure by low parasite diversity (`Evolutionary equilibrium` scenario), or there was a recent loss in MHC diversity leading to a lack of resistance alleles and increased parasite burden (`Unbalanced situation` scenario). In both species it became apparent that the MHC class II is functionally important in defence against gastrointestinal helminths, which was shown here for the first time in marsupials. On the population level, parasite diversity did not markedly differ between the two host species. However, we did observe considerable differences in the individual parasite load (parasite prevalence and infection intensity): while M. incanus revealed low MHC DAB diversity and high parasite load, G. microtarsus showed a tenfold higher population wide MHC DAB diversity and lower parasite burden. These results support the second scenario of an unbalanced situation.

Survival strategy of Echinococcus multilocularis in the human host

As exemplified by aborted calcified liver lesions commonly found in patients from endemic areas, Echinococcus multilocularis metacestodes develop only in a minority of individuals exposed to infection with the papasite. Clinical research has disclosed some aspects of the survival strategy of E. multilocularis in human hosts. Clinical observations in liver transplantation and AIDS suggest that suppression of cellular/Th1related immunity increases disease severity. Most of the studies have stressed a role for CD8+ T cells and for Interleukin-10 in the development of tolerance. A spontaneous secretion of IL-10 by the PBMC seems to be the immunological hallmark of patients with progressive forms of alveolar echinococcosis (AE). IL-10-induced inhibition of effector macrophages, but also of antigen-presenting dendritic cells, may be operating and allowing parasite growth and survival. The genetic correlates of susceptibility to infection with E. multilocularis are clearer in humans than in the mouse model. A significant link between MHC polymorphism and clinical presentation of AE has been shown, and the spontaneous secretion of IL-10 in patients with a progressive AE is higher in patients with the HLA DR3+, DQ2+ haplotype. Clustering of cases in certain families, in communities otherwise exposed to similar risk factors, also points to immuno-genetic predisposition factors that may allow the larva to escape host immunity more easily. The first stage of larval development may be crucial in producing danger signals stimulating the initial production of cytokines. Therapeutic use of Interferon alpha is an attempt to foil the survival strategy of E. multilocularis. (C) 2005 Elsevier Ireland Ltd. All rights reserved.

Visualisation of bioinformatics datasets

Analysing the molecular polymorphism and interactions of DNA, RNA and proteins is of fundamental importance in biology. Predicting functions of polymorphic molecules is important in order to design more effective medicines. Analysing major histocompatibility complex (MHC) polymorphism is important for mate choice, epitope-based vaccine design and transplantation rejection etc. Most of the existing exploratory approaches cannot analyse these datasets because of the large number of molecules with a high number of descriptors per molecule. This thesis develops novel methods for data projection in order to explore high dimensional biological dataset by visualising them in a low-dimensional space. With increasing dimensionality, some existing data visualisation methods such as generative topographic mapping (GTM) become computationally intractable. We propose variants of these methods, where we use log-transformations at certain steps of expectation maximisation (EM) based parameter learning process, to make them tractable for high-dimensional datasets. We demonstrate these proposed variants both for synthetic and electrostatic potential dataset of MHC class-I. We also propose to extend a latent trait model (LTM), suitable for visualising high dimensional discrete data, to simultaneously estimate feature saliency as an integrated part of the parameter learning process of a visualisation model. This LTM variant not only gives better visualisation by modifying the project map based on feature relevance, but also helps users to assess the significance of each feature. Another problem which is not addressed much in the literature is the visualisation of mixed-type data. We propose to combine GTM and LTM in a principled way where appropriate noise models are used for each type of data in order to visualise mixed-type data in a single plot. We call this model a generalised GTM (GGTM). We also propose to extend GGTM model to estimate feature saliencies while training a visualisation model and this is called GGTM with feature saliency (GGTM-FS). We demonstrate effectiveness of these proposed models both for synthetic and real datasets. We evaluate visualisation quality using quality metrics such as distance distortion measure and rank based measures: trustworthiness, continuity, mean relative rank errors with respect to data space and latent space. In cases where the labels are known we also use quality metrics of KL divergence and nearest neighbour classifications error in order to determine the separation between classes. We demonstrate the efficacy of these proposed models both for synthetic and real biological datasets with a main focus on the MHC class-I dataset.

MHC promoter polymorphism in grey wolves and domestic dogs

A functional immune system requires a tight control over major histocompatibility complex (MHC) gene transcription, as the abnormal MHC expression patterns of severe immunodeficiency and autoimmune diseases demonstrate. Although the regulation of MHC expression has been well documented in humans and mice, little is known in other species. In this study, we detail the level of polymorphism in wolf and dog MHC gene promoters. The promoter regions of the DRB, DQA and DQB locus were sequenced in 90 wolves and 90 dogs. The level of polymorphism was high in the DQB promoters, with variation found within functionally relevant regions, including binding sites for transcription factors. Clear associations between DQB promoters and exon 2 alleles were noted in wolves, indicating strong linkage disequilibrium in this region. Low levels of polymorphism were found within the DRB and DQA promoter regions. However, a variable site was identified within the T box, a TNF-alpha response element, of the DQA promoter. Furthermore, we identified a previously unrecognised 18-base-pair deletion within exon 1 of the DQB locus.

Major Histocompatibility Complex (MHC) Class II sequence polymorphism in long-finned pilot whale (Globicephala melas) from the North Atlantic

Funding Silvia S. Monteiro and Marisa Ferreira were supported by a Ph.D. grant from Fundação para a Ciência e Tecnologia (ref SFRH/BD/38735/2007 and SFRH/BD/30240/2006, respectively). Alfredo López was supported by a postdoctoral grant from Fundação para a Ciência e Tecnologia (ref SFRH/BPD/82407/2011). Catarina Eira is supported by CESAM (UID/AMB/50017), from FCT/MEC through national funds and FEDER (PT2020, Compete 2020). The work related with strandings and tissue collection in Portugal was partially supported by the SafeSea Project EEAGrants PT 0039 (supported by Iceland, Liechtenstein and Norway through the EEA Financial Mechanism), by the Project MarPro–Life09 NAT/PT/000038 (funded by the European Union–Program Life+) and by the project CetSenti FCT RECI/AAG-GLO/0470/2012; FCOMP-01-0124-FEDER-027472 (Funded by the Program COMPETE and Fundação para a Ciência e Tecnologia).

HLA-G 14-bp Insertion/Deletion Polymorphism Is a Risk Factor for HTLV-1 Infection

About 95% of HTLV-1 infected patients remain asymptomatic throughout life, and the risk factors associated with the development of related diseases, such as HAM/TSP and ATL, are not fully understood. The human leukocyte antigen-G molecule (HLA-G), a nonclassical HLA class I molecule encoded by MHC, is expressed in several pathological conditions, including viral infection, and is related to immunosuppressive effects that allow the virus-infected cells to escape the antiviral defense of the host. The 14-bp insertion/deletion polymorphism of exon 8 HLA-G gene influences the stability of the transcripts and could be related to HTLV-1-infected cell protection and to the increase of proviral load. The present study analyzed by conventional PCR the 14-bp insertion/deletion polymorphism of exon 8 HLA-G gene in 150 unrelated healthy subjects, 82 HTLV-1 infected patients with symptoms (33 ATL and 49 HAM), and 56 asymptomatic HTLV-1 infected patients (HAC). In addition, the proviral load was determined by quantitative real-time PCR in all infected groups and correlated with 14-bp insertion/deletion genotypes. The heterozygote genotype frequencies were significantly higher in HAM, in the symptomatic group, and in infected patients compared to control (p < 0.05). The proviral load was higher in the symptomatic group than the HAC group (p < 0.0005). The comparison of proviral load and genotypes showed that -14-bp/-14-bp genotype had a higher proviral load than +14-bp/-14-bp and +14-bp/+14-bp genotypes. Although HLA-G 14-bp polymorphism does not appear to be associated

Frequency of insertion/deletion polymorphism in exon 8 of HLA-G and kidney allograft outcome

Human leukocyte antigen-G (HLA-G) is a non-classical major histocompatibility complex (MHC) class Ib molecule predominantly expressed in cytotrophoblasts, where it acts as a specific immunosuppressor. Literature data have shown that grafts in some settings, such as cardiac and liver/kidney-associated transplantations, express HLA-G and this expression is associated with less severe rejection and also reduces the incidence of rejection. Fourteen-base pair deletion/insertion polymorphism has been reported in exon 8 of the 3`-untranslated region of HLA-G. This polymorphism within exon 8 of the HLA-G gene might influence transcription activity, which in turn may influence the stability of HLA-G transcripts. This influences the stability of the HLA-G protein and therefore is of potential functional relevance. In order to determine a possible correlation between the 14-bp insertion/deletion polymorphism and kidney allograft outcome, we isolated genomic DNA from 83 patients who had received isolated kidney allografts, and we classified the 83 specimens into two groups, grafts presenting Banff features of rejection group and a non-rejection group, and compared them with a control group of 97 healthy subjects. The 14-bp polymorphism at exon 8 was genotyped in all groups. There was no significant difference in allelic frequencies of 14-bp insertion/deletion polymorphism between normal controls and kidney transplant patients. In the RG, the homozygous genotype +14/+14 bp (P = 0.0238) was significantly increased in the group with acute rejection compared with the healthy control group. Analysis of other HLA-G polymorphisms and functional studies on immune regulation are essential to elucidate the role of HLA-G in kidney allografts.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection.

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.

Monitoring of NY-ESO-1 specific CD4+ T cells using molecularly defined MHC class II/His-tag-peptide tetramers.

MHC-peptide tetramers have become essential tools for T-cell analysis, but few MHC class II tetramers incorporating peptides from human tumor and self-antigens have been developed. Among limiting factors are the high polymorphism of class II molecules and the low binding capacity of the peptides. Here, we report the generation of molecularly defined tetramers using His-tagged peptides and isolation of folded MHC/peptide monomers by affinity purification. Using this strategy we generated tetramers of DR52b (DRB3*0202), an allele expressed by approximately half of Caucasians, incorporating an epitope from the tumor antigen NY-ESO-1. Molecularly defined tetramers avidly and stably bound to specific CD4(+) T cells with negligible background on nonspecific cells. Using molecularly defined DR52b/NY-ESO-1 tetramers, we could demonstrate that in DR52b(+) cancer patients immunized with a recombinant NY-ESO-1 vaccine, vaccine-induced tetramer-positive cells represent ex vivo in average 1:5,000 circulating CD4(+) T cells, include central and transitional memory polyfunctional populations, and do not include CD4(+)CD25(+)CD127(-) regulatory T cells. This approach may significantly accelerate the development of reliable MHC class II tetramers to monitor immune responses to tumor and self-antigens.

Dramatic influence of V beta gene polymorphism on an antigen-specific CD8+ T cell response in vivo.

According to recent crystallographic studies, the TCR-alpha beta contacts MHC class I-bound antigenic peptides via the polymorphic V gene-encoded complementarity-determining region 1 beta (CDR1 beta) and the hypervariable (D)J-encoded CDR3 beta and CDR3 alpha domains. To evaluate directly the relative importance of CDR1 beta polymorphism on the fine specificity of T cell responses in vivo, we have taken advantage of congenic V beta a and V beta b mouse strains that differ by a CDR1 polymorphism in the V beta 10 gene segment. The V beta 10-restricted CD8+ T cell response to a defined immunodominant epitope was dramatically reduced in V beta a compared with V beta b mice, as measured either by the expansion of V beta 10+ cells or by the binding of MHC-peptide tetramers. These data indicate that V beta polymorphism has an important impact on TCR-ligand binding in vivo, presumably by modifying the affinity of CDR1 beta-peptide interactions.

MHC-restricted antigen presentation and recognition: constraints on gene, recombinant and peptide vaccines in humans

The target of any immunization is to activate and expand lymphocyte clones with the desired recognition specificity and the necessary effector functions. In gene, recombinant and peptide vaccines, the immunogen is a single protein or a small assembly of epitopes from antigenic proteins. Since most immune responses against protein and peptide antigens are T-cell dependent, the molecular target of such vaccines is to generate at least 50-100 complexes between MHC molecule and the antigenic peptide per antigen-presenting cell, sensitizing a T cell population of appropriate clonal size and effector characteristics. Thus, the immunobiology of antigen recognition by T cells must be taken into account when designing new generation peptide- or gene-based vaccines. Since T cell recognition is MHC-restricted, and given the wide polymorphism of the different MHC molecules, distinct epitopes may be recognized by different individuals in the population. Therefore, the issue of whether immunization will be effective in inducing a protective immune response, covering the entire target population, becomes an important question. Many pathogens have evolved molecular mechanisms to escape recognition by the immune system by variation of antigenic protein sequences. In this short review, we will discuss the several concepts related to selection of amino acid sequences to be included in DNA and peptide vaccines.