Large-scale validation of methods for cytotoxic T-lymphocyte epitope prediction

Background  

Reliable predictions of Cytotoxic T lymphocyte (CTL) epitopes are essential for rational vaccine design. Most importantly, they can minimize the experimental effort needed to identify epitopes. NetCTL is a web-based tool designed for predicting human CTL epitopes in any given protein. It does so by integrating predictions of proteasomal cleavage, TAP transport efficiency, and MHC class I affinity. At least four other methods have been developed recently that likewise attempt to predict CTL epitopes: EpiJen, MAPPP, MHC-pathway, and WAPP. In order to compare the performance of prediction methods, objective benchmarks and standardized performance measures are needed. Here, we develop such large-scale benchmark and corresponding performance measures and report the performance of an updated version 1.2 of NetCTL in comparison with the four other methods.     

NN-align. An artificial neural network-based alignment algorithm for MHC class II peptide binding prediction

Background  

The major histocompatibility complex (MHC) molecule plays a central role in controlling the adaptive immune response to infections. MHC class I molecules present peptides derived from intracellular proteins to cytotoxic T cells, whereas MHC class II molecules stimulate cellular and humoral immunity through presentation of extracellularly derived peptides to helper T cells. Identification of which peptides will bind a given MHC molecule is thus of great importance for the understanding of host-pathogen interactions, and large efforts have been placed in developing algorithms capable of predicting this binding event.     

Prediction of MHC class I binding peptides, using SVMHC

Background  

T-cells are key players in regulating a specific immune response. Activation of cytotoxic T-cells requires recognition of specific peptides bound to Major Histocompatibility Complex (MHC) class I molecules. MHC-peptide complexes are potential tools for diagnosis and treatment of pathogens and cancer, as well as for the development of peptide vaccines. Only one in 100 to 200 potential binders actually binds to a certain MHC molecule, therefore a good prediction method for MHC class I binding peptides can reduce the number of candidate binders that need to be synthesized and tested.     

Natural Killer Cell Tolerance Persists Despite Significant Reduction of Self MHC Class I on Normal Target Cells in Mice

Background

A major group of murine inhibitory receptors on Natural Killer (NK) cells belong to the Ly49 receptor family and recognize MHC class I molecules. Infected or transformed target cells frequently downmodulate MHC class I molecules and can thus avoid CD8⁺ T cell attack, but may at the same time develop NK cell sensitivity, due to failure to express inhibitory ligands for Ly49 receptors. The extent of MHC class I downregulation needed on normal cells to trigger NK cell effector functions is not known.

Methodology/Principal Findings

In this study, we show that cells expressing MHC class I to levels well below half of the host level are tolerated in an in vivo assay in mice. Hemizygous expression (expression from only one allele) of MHC class I was sufficient to induce Ly49 receptor downmodulation on NK cells to a similar degree as homozygous expression, despite a strongly reduced cell surface level of MHC class I. Co-expression of weaker MHC class I ligands in the host did not have any further effect on the degree of Ly49 downmodulation. Furthermore, a single MHC class I allele could downmodulate up to three Ly49 receptors on individual NK cells. Only when NK cells simultaneously expressed several Ly49 receptors and hemizygous MHC class I levels, a putative threshold for Ly49 downmodulation was reached.

Conclusion

Collectively, our findings suggest that in interactions between NK cells and normal untransformed cells, MHC class I molecules are in most cases expressed in excess compared to what is functionally needed to ensure self tolerance and to induce maximal Ly49 downmodulation. We speculate that the reason for this is to maintain a safety margin for otherwise normal, autologous cells over a range of MHC class I expression levels, in order to ensure robustness in NK cell tolerance.     

Various hypotheses on MHC evolution suggested by the concerted evolution of CD94L and MHC class Ia molecules

Background  

In the accompanying paper by Virginie Rouillon and myself, our demonstration that homogenisation by gene conversion occurs readily among MHC class I genes was made possible because of the exceptional conservation of the CD94L locus between divergent species of separate taxa, suggesting that the molecules of this family are endowed with very important and well preserved biological functions. These results lead me to elaborate various hypotheses on several aspects of MHC evolution.     

Systematic Characterisation of Cellular Localisation and Expression Profiles of Proteins Containing MHC Ligands

Background

Presentation of peptides on Major Histocompatibility Complex (MHC) molecules is the cornerstone in immune system activation and increased knowledge of the characteristics of MHC ligands and their source proteins is highly desirable.

Methodology/Principal Finding

In the present large-scale study, we used a large data set of proteins containing experimentally identified MHC class I or II ligands and examined the proteins according to their expression profiles at the mRNA level and their Gene Ontology (GO) classification within the cellular component ontology. Proteins encoded by highly abundant mRNA were found to be much more likely to be the source of MHC ligands. Of the 2.5% most abundant mRNAs as much as 41% of the proteins encoded by these mRNAs contained MHC class I ligands. For proteins containing MHC class II ligands, the corresponding percentage was 11%. Furthermore, we found that most proteins containing MHC class I ligands were localised to the intracellular parts of the cell including the cytoplasm and nucleus. MHC class II ligand donors were, on the other hand, mostly membrane proteins.

Conclusions/Significance

The results contribute to the ongoing debate concerning the nature of MHC ligand-containing proteins and can be used to extend the existing methods for MHC ligand predictions by including the source protein''s localisation and expression profile. Improving the current methods is important in the growing quest for epitopes that can be used for vaccine or diagnostic purposes, especially when it comes to large DNA viruses and cancer.     

A mitochondrial mutator plasmid that causes senescence under dietary restricted conditions

Background

The definition of human MHC class I haplotypes through association of HLA-A, HLA-Cw and HLA-B has been used to analyze ethnicity, population migrations and disease association.

Results

Here, we present HLA-E allele haplotype association and population linkage disequilibrium (LD) analysis within the ~1.3 Mb bounded by HLA-B/Cw and HLA-A to increase the resolution of identified class I haplotypes. Through local breakdown of LD, we inferred ancestral recombination points both upstream and downstream of HLA-E contributing to alternative block structures within previously identified haplotypes. Through single nucleotide polymorphism (SNP) analysis of the MHC region, we also confirmed the essential genetic fixity, previously inferred by MHC allele analysis, of three conserved extended haplotypes (CEHs), and we demonstrated that commercially-available SNP analysis can be used in the MHC to help define CEHs and CEH fragments.

Conclusion

We conclude that to generate high-resolution maps for relating MHC haplotypes to disease susceptibility, both SNP and MHC allele analysis must be conducted as complementary techniques.     

Gene duplication and fragmentation in the zebra finch major histocompatibility complex

Background

Due to its high polymorphism and importance for disease resistance, the major histocompatibility complex (MHC) has been an important focus of many vertebrate genome projects. Avian MHC organization is of particular interest because the chicken Gallus gallus, the avian species with the best characterized MHC, possesses a highly streamlined minimal essential MHC, which is linked to resistance against specific pathogens. It remains unclear the extent to which this organization describes the situation in other birds and whether it represents a derived or ancestral condition. The sequencing of the zebra finch Taeniopygia guttata genome, in combination with targeted bacterial artificial chromosome (BAC) sequencing, has allowed us to characterize an MHC from a highly divergent and diverse avian lineage, the passerines.

Results

The zebra finch MHC exhibits a complex structure and history involving gene duplication and fragmentation. The zebra finch MHC includes multiple Class I and Class II genes, some of which appear to be pseudogenes, and spans a much more extensive genomic region than the chicken MHC, as evidenced by the presence of MHC genes on each of seven BACs spanning 739 kb. Cytogenetic (FISH) evidence and the genome assembly itself place core MHC genes on as many as four chromosomes with TAP and Class I genes mapping to different chromosomes. MHC Class II regions are further characterized by high endogenous retroviral content. Lastly, we find strong evidence of selection acting on sites within passerine MHC Class I and Class II genes.

Conclusion

The zebra finch MHC differs markedly from that of the chicken, the only other bird species with a complete genome sequence. The apparent lack of synteny between TAP and the expressed MHC Class I locus is in fact reminiscent of a pattern seen in some mammalian lineages and may represent convergent evolution. Our analyses of the zebra finch MHC suggest a complex history involving chromosomal fission, gene duplication and translocation in the history of the MHC in birds, and highlight striking differences in MHC structure and organization among avian lineages.     

HIV-1 Envelope Resistance to Proteasomal Cleavage: Implications for Vaccine Induced Immune Responses

Background

Antigen processing involves many proteolytic enzymes such as proteasomes and cathepsins. The processed antigen is then presented on the cell surface bound to either MHC class I or class II molecules and induces/interacts with antigen-specific CD8+ and CD4+ T-cells, respectively. Preliminary immunological data from the RV144 phase III trial indicated that the immune responses were biased towards the Env antigen with a dominant CD4+ T-cell response.

Methods

In this study, we examined the susceptibility of HIV-1 Env-A244 gp120 protein, one of the protein boost subunits of the RV144 Phase III vaccine trial, to proteasomes and cathepsins and identified the generated peptide epitope repertoire by mass spectrometry. The peptide fragments were tested for cytokine production in CD4⁺ T-cell lines derived from RV144 volunteers.

Results

Env-A244 was resistant to proteasomes, thus diminishing the possibility of the generation of class I epitopes by the classical MHC class I pathway. However, Env-A244 was efficiently cleaved by cathepsins generating peptide arrays identified by mass spectrometry that contained both MHC class I and class II epitopes as reported in the Los Alamos database. Each of the cathepsins generated distinct degradation patterns containing regions of light and dense epitope clusters. The sequence DKKQKVHALF that is part of the V2 loop of gp120 produced by cathepsins induced a polyfunctional cytokine response including the generation of IFN-γ from CD4⁺ T-cell lines-derived from RV144 vaccinees. This sequence is significant since antibodies to the V1/V2-loop region correlated inversely with HIV-1 infection in the RV144 trial.

Conclusions

Based on our results, the susceptibility of Env-A244 to cathepsins and not to proteasomes suggests a possible mechanism for the generation of Env-specific CD4⁺T cell and antibody responses in the RV144 vaccinees.     

Immunogenetic Variation and Differential Pathogen Exposure in Free-Ranging Cheetahs across Namibian Farmlands

Background

Genes under selection provide ecologically important information useful for conservation issues. Major histocompatibility complex (MHC) class I and II genes are essential for the immune defence against pathogens from intracellular (e.g. viruses) and extracellular (e.g. helminths) origins, respectively. Serosurvey studies in Namibian cheetahs (Acinonyx juabuts) revealed higher exposure to viral pathogens in individuals from north-central than east-central regions. Here we examined whether the observed differences in exposure to viruses influence the patterns of genetic variation and differentiation at MHC loci in 88 free-ranging Namibian cheetahs.

Methodology/Principal Findings

Genetic variation at MHC I and II loci was assessed through single-stranded conformation polymorphism (SSCP) analysis and sequencing. While the overall allelic diversity did not differ, we observed a high genetic differentiation at MHC class I loci between cheetahs from north-central and east-central Namibia. No such differentiation in MHC class II and neutral markers were found.

Conclusions/Significance

Our results suggest that MHC class I variation mirrors the variation in selection pressure imposed by viruses in free-ranging cheetahs across Namibian farmland. This is of high significance for future management and conservation programs of this species.     

Disruption and pseudoautosomal localization of the major histocompatibility complex in monotremes

Background

The monotremes, represented by the duck-billed platypus and the echidnas, are the most divergent species within mammals, featuring a flamboyant mix of reptilian, mammalian and specialized characteristics. To understand the evolution of the mammalian major histocompatibility complex (MHC), the analysis of the monotreme genome is vital.

Results

We characterized several MHC containing bacterial artificial chromosome clones from platypus (Ornithorhynchus anatinus) and the short-beaked echidna (Tachyglossus aculeatus) and mapped them onto chromosomes. We discovered that the MHC of monotremes is not contiguous and locates within pseudoautosomal regions of two pairs of their sex chromosomes. The analysis revealed an MHC core region with class I and class II genes on platypus and echidna X3/Y3. Echidna X4/Y4 and platypus Y4/X5 showed synteny to the human distal class III region and beyond. We discovered an intron-containing class I pseudogene on platypus Y4/X5 at a genomic location equivalent to the human HLA-B,C region, suggesting ancestral synteny of the monotreme MHC. Analysis of male meioses from platypus and echidna showed that MHC chromosomes occupy different positions in the meiotic chains of either species.

Conclusion

Molecular and cytogenetic analyses reveal new insights into the evolution of the mammalian MHC and the multiple sex chromosome system of monotremes. In addition, our data establish the first homology link between chicken microchromosomes and the smallest chromosomes in the monotreme karyotype. Our results further suggest that segments of the monotreme MHC that now reside on separate chromosomes must once have been syntenic and that the complex sex chromosome system of monotremes is dynamic and still evolving.     

Epigenetic Mechanisms Regulate MHC and Antigen Processing Molecules in Human Embryonic and Induced Pluripotent Stem Cells

Background

Human embryonic stem cells (hESCs) are an attractive resource for new therapeutic approaches that involve tissue regeneration. hESCs have exhibited low immunogenicity due to low levels of Mayor Histocompatibility Complex (MHC) class-I and absence of MHC class-II expression. Nevertheless, the mechanisms regulating MHC expression in hESCs had not been explored.

Methodology/Principal Findings

We analyzed the expression levels of classical and non-classical MHC class-I, MHC class-II molecules, antigen-processing machinery (APM) components and NKG2D ligands (NKG2D-L) in hESCs, induced pluripotent stem cells (iPSCs) and NTera2 (NT2) teratocarcinoma cell line. Epigenetic mechanisms involved in the regulation of these genes were investigated by bisulfite sequencing and chromatin immunoprecipitation (ChIP) assays. We showed that low levels of MHC class-I molecules were associated with absent or reduced expression of the transporter associated with antigen processing 1 (TAP-1) and tapasin (TPN) components in hESCs and iPSCs, which are involved in the transport and load of peptides. Furthermore, lack of β2-microglobulin (β2m) light chain in these cells limited the expression of MHC class I trimeric molecule on the cell surface. NKG2D ligands (MICA, MICB) were observed in all pluripotent stem cells lines. Epigenetic analysis showed that H3K9me3 repressed the TPN gene in undifferentiated cells whilst HLA-B and β2m acquired the H3K4me3 modification during the differentiation to embryoid bodies (EBs). Absence of HLA-DR and HLA-G expression was regulated by DNA methylation.

Conclusions/Significance

Our data provide fundamental evidence for the epigenetic control of MHC in hESCs and iPSCs. Reduced MHC class I and class II expression in hESCs and iPSCs can limit their recognition by the immune response against these cells. The knowledge of these mechanisms will further allow the development of strategies to induce tolerance and improve stem cell allograft acceptance.     

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background  

MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions.     

Female house sparrows "count on" male genes: experimental evidence for MHC-dependent mate preference in birds

Background  

Females can potentially assess the quality of potential mates using their secondary sexual traits, and obtain "good genes" that increase offspring fitness. Another potential indirect benefit from mating preferences is genetic compatibility, which does not require extravagant or viability indicator traits. Several studies with mammals and fish indicate that the genes of the major histocompatibility complex (MHC) influence olfactory cues and mating preferences, and such preferences confer genetic benefits to offspring. We investigated whether individual MHC diversity (class I) influences mating preferences in house sparrows (Passer domesticus).     

Prediction of MHC class II binding affinity using SMM-align,a novel stabilization matrix alignment method

Background  

Antigen presenting cells (APCs) sample the extra cellular space and present peptides from here to T helper cells, which can be activated if the peptides are of foreign origin. The peptides are presented on the surface of the cells in complex with major histocompatibility class II (MHC II) molecules. Identification of peptides that bind MHC II molecules is thus a key step in rational vaccine design and developing methods for accurate prediction of the peptide:MHC interactions play a central role in epitope discovery. The MHC class II binding groove is open at both ends making the correct alignment of a peptide in the binding groove a crucial part of identifying the core of an MHC class II binding motif. Here, we present a novel stabilization matrix alignment method, SMM-align, that allows for direct prediction of peptide:MHC binding affinities. The predictive performance of the method is validated on a large MHC class II benchmark data set covering 14 HLA-DR (human MHC) and three mouse H2-IA alleles.     

Marsupials and monotremes possess a novel family of MHC class I genes that is lost from the eutherian lineage

Background

Major histocompatibility complex (MHC) class I genes are found in the genomes of all jawed vertebrates. The evolution of this gene family is closely tied to the evolution of the vertebrate genome. Family members are frequently found in four paralogous regions, which were formed in two rounds of genome duplication in the early vertebrates, but in some species class Is have been subject to additional duplication or translocation, creating additional clusters. The gene family is traditionally grouped into two subtypes: classical MHC class I genes that are usually MHC-linked, highly polymorphic, expressed in a broad range of tissues and present endogenously-derived peptides to cytotoxic T-cells; and non-classical MHC class I genes generally have lower polymorphism, may have tissue-specific expression and have evolved to perform immune-related or non-immune functions. As immune genes can evolve rapidly and are subject to different selection pressure, we hypothesised that there may be divergent, as yet unannotated or uncharacterised class I genes.

Results

Application of a novel method of sensitive genome searching of available vertebrate genome sequences revealed a new, extensive sub-family of divergent MHC class I genes, denoted as UT, which has not previously been characterized. These class I genes are found in both American and Australian marsupials, and in monotremes, at an evolutionary chromosomal breakpoint, but are not present in non-mammalian genomes and have been lost from the eutherian lineage. We show that UT family members are expressed in the thymus of the gray short-tailed opossum and in other immune tissues of several Australian marsupials. Structural homology modelling shows that the proteins encoded by this family are predicted to have an open, though short, antigen-binding groove.

Conclusions

We have identified a novel sub-family of putatively non-classical MHC class I genes that are specific to marsupials and monotremes. This family was present in the ancestral mammal and is found in extant marsupials and monotremes, but has been lost from the eutherian lineage. The function of this family is as yet unknown, however, their predicted structure may be consistent with presentation of antigens to T-cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1745-4) contains supplementary material, which is available to authorized users.     

Predicting peptides binding to MHC class II molecules using multi-objective evolutionary algorithms

Background  

Peptides binding to Major Histocompatibility Complex (MHC) class II molecules are crucial for initiation and regulation of immune responses. Predicting peptides that bind to a specific MHC molecule plays an important role in determining potential candidates for vaccines. The binding groove in class II MHC is open at both ends, allowing peptides longer than 9-mer to bind. Finding the consensus motif facilitating the binding of peptides to a MHC class II molecule is difficult because of different lengths of binding peptides and varying location of 9-mer binding core. The level of difficulty increases when the molecule is promiscuous and binds to a large number of low affinity peptides.