T-cell epitope prediction: rescaling can mask biological variation between MHC molecules

Novel epitopes identified from Tembusu virus NS3 protein induce cytotoxic T lymphocyte response

Since 2010, Tembusu virus (TMUV) has spread widely in China, causing huge economic losses to the poultry industry. Due to the infectious and zoonotic nature of flaviviruses, their potential threat to public health is of great concern. Cellular immune responses usually play a critical role in combating viral infections. To study the molecular basis of cell immunity induced by TMUV, 14 cytotoxic T lymphocyte (CTL) epitope peptides of TMUV antigen E, NS1 and NS3 were predicted by bioinformatics tools. Their abilities to induce cellular immune responses were determined by IFN-γ ELISpot assay, and 4 peptides were found to exhibit highly significant responses upon stimulation. In addition, the cytotoxic activity induced by the epitope peptides was assessed by lactate dehydrogenase (LDH) release assay. Finally, among these peptides, we identified two murine TMUV NS3-derived H-2d-restricted CTL epitopes in BALB/c mice, which could be used to further study of epitope vaccines against TMUV infection.

Host genotype and time dependent antigen presentation of viral peptides: predictions from theory

The rate of progression of HIV infected individuals to AIDS is known to vary with the genotype of the host, and is linked to their allele of human leukocyte antigen (HLA) proteins, which present protein degradation products at the cell surface to circulating T-cells. HLA alleles are associated with Gag-specific T-cell responses that are protective against progression of the disease. While Pol is the most conserved HIV sequence, its association with immune control is not as strong. To gain a more thorough quantitative understanding of the factors that contribute to immunodominance, we have constructed a model of the recognition of HIV infection by the MHC class I pathway. Our model predicts surface presentation of HIV peptides over time, demonstrates the importance of viral protein kinetics, and provides evidence of the importance of Gag peptides in the long-term control of HIV infection. Furthermore, short-term dynamics are also predicted, with simulation of virion-derived peptides suggesting that efficient processing of Gag can lead to a 50% probability of presentation within 3 hours post-infection, as observed experimentally. In conjunction with epitope prediction algorithms, this modelling approach could be used to refine experimental targets for potential T-cell vaccines, both for HIV and other viruses.

Genetic Markers of the Host in Persons Living with HTLV-1, HIV and HCV Infections

Human T-cell leukemia virus type 1 (HTLV-1), hepatitis C virus (HCV) and human immunodeficiency virus type 1 (HIV-1) are prevalent worldwide, and share similar means of transmission. These infections may influence each other in evolution and outcome, including cancer or immunodeficiency. Many studies have reported the influence of genetic markers on the host immune response against different persistent viral infections, such as HTLV-1 infection, pointing to the importance of the individual genetic background on their outcomes. However, despite recent advances on the knowledge of the pathogenesis of HTLV-1 infection, gaps in the understanding of the role of the individual genetic background on the progress to disease clinically manifested still remain. In this scenario, much less is known regarding the influence of genetic factors in the context of dual or triple infections or their influence on the underlying mechanisms that lead to outcomes that differ from those observed in monoinfection. This review describes the main factors involved in the virus–host balance, especially for some particular human leukocyte antigen (HLA) haplotypes, and other important genetic markers in the development of HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) and other persistent viruses, such as HIV and HCV.

Complementarity of Binding Motifs is a General Property of HLA-A and HLA-B Molecules and Does Not Seem to Effect HLA Haplotype Composition

Different human leukocyte antigen (HLA) haplotypes (i.e., the specific combinations of HLA-A, -B, -DR alleles inherited together from one parent) are observed in different frequencies in human populations. Some haplotypes, like HLA-A1-B8, are very frequent, reaching up to 10% in the Caucasian population, while others are very rare. Numerous studies have identified associations between HLA haplotypes and diseases, and differences in haplotype frequencies can in part be explained by these associations: the stronger the association with a severe (autoimmune) disease, the lower the expected HLA haplotype frequency. The peptide repertoires of the HLA molecules composing a haplotype can also influence the frequency of a haplotype. For example, it would seem advantageous to have HLA molecules with non-overlapping binding specificities within a haplotype, as individuals expressing such an haplotype would present a diverse set of peptides from viruses and pathogenic bacteria on the cell surface. To test this hypothesis, we collect the proteome data from a set of common viruses, and estimate the total ligand repertoire of HLA class I haplotypes (HLA-A-B) using in silico predictions. We compare the size of these repertoires to the HLA haplotype frequencies reported in the National Marrow Donor Program (NMDP). We find that in most HLA-A and HLA-B pairs have fairly distinct binding motifs, and that the observed haplotypes do not contain HLA-A and -B molecules with more distinct binding motifs than random HLA-A and HLA-B pairs. In addition, the population frequency of a haplotype is not correlated to the distinctness of its HLA-A and HLA-B peptide binding motifs. These results suggest that there is a not a strong selection pressure on the haplotype level favoring haplotypes having HLA molecules with distinct binding motifs, which would result the largest possible presented peptide repertoires in the context of infectious diseases.

Deep sequencing of viral genomes provides insight into the evolution and pathogenesis of varicella zoster virus and its vaccine in humans

Immunization with the vOka vaccine prevents varicella (chickenpox) in children and susceptible adults. The vOka vaccine strain comprises a mixture of genotypes and, despite attenuation, causes rashes in small numbers of recipients. Like wild-type virus, the vaccine establishes latency in neuronal tissue and can later reactivate to cause Herpes zoster (shingles). Using hybridization-based methodologies, we have purified and sequenced vOka directly from skin lesions. We show that alleles present in the vaccine can be recovered from the lesions and demonstrate the presence of a severe bottleneck between inoculation and lesion formation. Genotypes in any one lesion appear to be descended from one to three vaccine-genotypes with a low frequency of novel mutations. No single vOka haplotype and no novel mutations are consistently present in rashes, indicating that neither new mutations nor recombination with wild type are critical to the evolution of vOka rashes. Instead, alleles arising from attenuation (i.e., not derived from free-living virus) are present at lower frequencies in rash genotypes. We identify 11 loci at which the ancestral allele is selected for in vOka rash formation and show genotypes in rashes that have reactivated from latency cannot be distinguished from rashes occurring immediately after inoculation. We conclude that the vOka vaccine, although heterogeneous, has not evolved to form rashes through positive selection in the mode of a quasispecies, but rather alleles that were essentially neutral during the vaccine production have been selected against in the human subjects, allowing us to identify key loci for rash formation.

Sequence-based in silico analysis of well studied hepatitis C virus epitopes and their variants in other genotypes (particularly genotype 5a) against South African human leukocyte antigen backgrounds

Background

Host genetics influence the outcome of HCV disease. HCV is also highly mutable and escapes host immunity. HCV genotypes are geographically distributed and HCV subtypes have been shown to have distinct repertoires of HLA-restricted viral epitopes which explains the lack of cross protection across genotypes observed in some studies. Despite this, immune databases and putative epitope vaccines concentrate almost exclusively on HCV genotype 1 class I-epitopes restricted by the HLA-A*02 allele. While both genotype and allele predominate in developed countries, we hypothesise that HCV variation and population genetics will affect the efficacy of proposed epitope vaccines in South Africa. This in silico study investigates HCV viral variability within well-studied epitopes identified in genotype 1 and uses algorithms to predict the immunogenicity of their variants from other less studied genotypes and thus rate the most promising vaccine candidates for the South African population. Six class I- and seven class II- restricted epitope sequences within the core, NS3, NS4B and NS5B regions were compared across the six HCV genotypes using local genotype 5a sequence data together with global data. Common HLA alleles in the South African population are A30:01, A02:01, B58:02, B07:02; DRB1*13:01 and DRB1*03:01. Epitope binding to 13 class I- and 8 class –II alleles were described using web-based prediction servers, Immune Epitope Database, (IEDB) and Propred. Online population coverage tools were used to assess vaccine efficacy.

Results

Despite the homogeneity of genotype 1 and genotype 5 over the epitopes, there was limited promiscuity to local HLA-alleles.Host differences will make a putative vaccine less effective in South Africa. Of the 6 well-characterized class I- epitopes, only 2 class I- epitopes were promiscuous and 3 of the 7 class-II epitopes were better conserved and promiscuous. By fine tuning the putative vaccine using an optimal cocktail of genotype 1 and 5a epitopes and local HLA data, the coverage was raised from 65.85% to 91.87% in South African Blacks.

Conclusion

While in vivo and in vitro studies are needed to confirm immunogenic epitopes, in silico HCV epitope vaccine design which takes into account HCV variation and host allele frequency will maximize population coverage in different ethnic groups.

Estimating the fitness cost of escape from HLA presentation in HIV-1 protease and reverse transcriptase

Human immunodeficiency virus (HIV-1) is, like most pathogens, under selective pressure to escape the immune system of its host. In particular, HIV-1 can avoid recognition by cytotoxic T lymphocytes (CTLs) by altering the binding affinity of viral peptides to human leukocyte antigen (HLA) molecules, the role of which is to present those peptides to the immune system. It is generally assumed that HLA escape mutations carry a replicative fitness cost, but these costs have not been quantified. In this study, we assess the replicative cost of mutations which are likely to escape presentation by HLA molecules in the region of HIV-1 protease and reverse transcriptase. Specifically, we combine computational approaches for prediction of in vitro replicative fitness and peptide binding affinity to HLA molecules. We find that mutations which impair binding to HLA-A molecules tend to have lower in vitro replicative fitness than mutations which do not impair binding to HLA-A molecules, suggesting that HLA-A escape mutations carry higher fitness costs than non-escape mutations. We argue that the association between fitness and HLA-A binding impairment is probably due to an intrinsic cost of escape from HLA-A molecules, and these costs are particularly strong for HLA-A alleles associated with efficient virus control. Counter-intuitively, we do not observe a significant effect in the case of HLA-B, but, as discussed, this does not argue against the relevance of HLA-B in virus control. Overall, this article points to the intriguing possibility that HLA-A molecules preferentially target more conserved regions of HIV-1, emphasizing the importance of HLA-A genes in the evolution of HIV-1 and RNA viruses in general.

Our immune system can recognize and kill virus-infected cells by distinguishing between self and virus-derived protein fragments, called peptides, displayed on the surface of each cell. One requirement for a successful recognition is that those peptides bind to the human leukocyte antigen (HLA) class I molecules, which present them to the immune system. As a counter-strategy, human immunodeficiency virus type 1 (HIV-1) can acquire mutations that prevent this binding, thereby helping the virus to escape the surveillance of T-lymphocytes. It is likely that the virus pays a replicative cost for such escape mutations, but the magnitude of this cost has remained elusive. Here, we quantified this fitness cost in HIV-1 protease and reverse transcriptase by combining two computational systems biology approaches: one for prediction of in vitro replicative fitness, and one for the prediction of the efficiency of peptide binding to HLA. We found that in viral proteins targeted by HLA-A molecules, mutations which disrupt binding to those molecules carry a lower replicative fitness than mutations which do not have such an effect. We argue that these results are consistent with the hypothesis that our immune systems might have evolved to target genetic regions of RNA viruses which are costly for the pathogen to alter.

Prediction of MHC class I binding peptides with a new feature encoding technique

The recognition of specific peptides, bound to major histocompatibility complex (MHC) class I molecules, is of particular importance to the robust identification of T-cell epitopes and thus the successful design of protein-based vaccines. Here, we present a new feature amino acid encoding technique termed OEDICHO to predict MHC class I/peptide complexes. In the proposed method, we have combined orthonormal encoding (OE) and the binary representation of selected 10 best physicochemical properties of amino acids derived from Amino Acid Index Database (AAindex). We also have compared our method to current feature encoding techniques. The tests have been carried out on comparatively large Human Leukocyte Antigen (HLA)-A and HLA-B allele peptide binding datasets. Empirical results show that our amino acid encoding scheme leads to better classification performance on a standalone classifier.

Immunoinformatic evaluation of multiple epitope ensembles as vaccine candidates: E coli 536

Epitope prediction is becoming a key tool for vaccine discovery. Prospective analysis of bacterial and viral genomes can identify antigenic epitopes encoded within individual genes that may act as effective vaccines against specific pathogens. Since B-cell epitope prediction remains unreliable, we concentrate on T-cell epitopes, peptides which bind with high affinity to Major Histacompatibility Complexes (MHC). In this report, we evaluate the veracity of identified T-cell epitope ensembles, as generated by a cascade of predictive algorithms (SignalP, Vaxijen, MHCPred, IDEB, EpiJen), as a candidate vaccine against the model pathogen uropathogenic gram negative bacteria Escherichia coli (E-coli) strain 536 (O6:K15:H31). An immunoinformatic approach was used to identify 23 epitopes within the E-coli proteome. These epitopes constitute the most promiscuous antigenic sequences that bind across more than one HLA allele with high affinity (IC50 < 50nM). The reliability of software programmes used, polymorphic nature of genes encoding MHC and what this means for population coverage of this potential vaccine are discussed.

The efficiency of the human CD8+ T cell response: how should we quantify it, what determines it, and does it matter?

Multidisciplinary techniques, in particular the combination of theoretical and experimental immunology, can address questions about human immunity that cannot be answered by other means. From the turnover of virus-infected cells in vivo, to rates of thymic production and HLA class I epitope prediction, theoretical techniques provide a unique insight to supplement experimental approaches. Here we present our opinion, with examples, of some of the ways in which mathematics has contributed in our field of interest: the efficiency of the human CD8+ T cell response to persistent viruses.

Is There a Role for HTLV-1-Specific CTL in Adult T-Cell Leukemia/Lymphoma?

ATLL is an aggressive malignancy of T cells that affects about 5% of individuals infected with HTLV-1. The precise mechanism of oncogenesis is not known, but there is evidence that two regulatory viral proteins, Tax and HBZ, are involved. A high set point proviral load is associated with development of ATLL or a chronic inflammatory condition, HAM/TSP. Several lines of evidence, including HLA class 1 association studies and in vitro killing assays, indicate that cytotoxic T lymphocytes are instrumental in determining this proviral load set point. Prior studies have focused chiefly on the CTL response to the immunodominant Tax protein: efficient lysis of Tax-expressing cells inversely correlates with proviral load in nonmalignant infection. However, a recent study showed that strong binding of peptides from HBZ, but not Tax, to HLA class 1 molecules was associated with a low proviral load and a reduced risk of developing HAM/TSP, indicating an important role for HBZ-specific CTL in determining infection outcome. In comparison with nonmalignant infection, HTLV-1-specific CTLs in ATLL patients are reduced in frequency and functionally deficient. Here we discuss the nature of protective CTL responses in nonmalignant HTLV-1 infection and explore the potential of CTLs to protect against ATLL.

Short communication an interferon-γ ELISPOT assay with two cytotoxic T cell epitopes derived from HTLV-1 tax region 161-233 discriminates HTLV-1-associated myelopathy/tropical spastic paraparesis patients from asymptomatic HTLV-1 carriers in a Peruvian population

HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) is a chronic and progressive disorder caused by the human T-lymphotropic virus type 1 (HTLV-1). In HTLV-1 infection, a strong cytotoxic T cell (CTL) response is mounted against the immunodominant protein Tax. Previous studies carried out by our group reported that increased IFN-γ enzyme-linked immunospot (ELISPOT) responses against the region spanning amino acids 161 to 233 of the Tax protein were associated with HAM/TSP and increased HTLV-1 proviral load (PVL). An exploratory study was conducted on 16 subjects with HAM/TSP, 13 asymptomatic carriers (AC), and 10 HTLV-1-seronegative controls (SC) to map the HAM/TSP-associated CTL epitopes within Tax region 161-233. The PVL of the infected subjects was determined and the specific CTL response was evaluated with a 6-h incubation IFN-γ ELISPOT assay using peripheral blood mononuclear cells (PBMCs) stimulated with 16 individual overlapping peptides covering the Tax region 161-233. Other proinflammatory and Th1/Th2 cytokines were also quantified in the supernatants by a flow cytometry multiplex assay. In addition, a set of human leukocyte antigen (HLA) class I alleles that bind with high affinity to the CTL epitopes of interest was determined using computational tools. Univariate analyses identified an association between ELISPOT responses to two new CTL epitopes, Tax 173-185 and Tax 181-193, and the presence of HAM/TSP as well as an increased PVL. The HLA-A*6801 allele, which is predicted to bind to the Tax 181-193 peptide, was overpresented in the HAM/TSP patients tested.

KIR2DL2 enhances protective and detrimental HLA class I-mediated immunity in chronic viral infection

Killer cell immunoglobulin-like receptors (KIRs) influence both innate and adaptive immunity. But while the role of KIRs in NK-mediated innate immunity is well-documented, the impact of KIRs on the T cell response in human disease is not known. Here we test the hypothesis that an individual's KIR genotype affects the efficiency of their HLA class I-mediated antiviral immune response and the outcome of viral infection. We show that, in two unrelated viral infections, hepatitis C virus and human T lymphotropic virus type 1, possession of the KIR2DL2 gene enhanced both protective and detrimental HLA class I-restricted anti-viral immunity. These results reveal a novel role for inhibitory KIRs. We conclude that inhibitory KIRs, in synergy with T cells, are a major determinant of the outcome of persistent viral infection.

HLA class I allele promiscuity revisited

The peptide repertoire presented on human leukocyte antigen (HLA) class I molecules is largely determined by the structure of the peptide binding groove. It is expected that the molecules having similar grooves (i.e., belonging to the same supertype) might present similar/overlapping peptides. However, the extent of promiscuity among HLA class I ligands remains controversial: while in many studies T cell responses are detected against epitopes presented by alternative molecules across HLA class I supertypes and loci, peptide elution studies report minute overlaps between the peptide repertoires of even related HLA molecules. To get more insight into the promiscuous peptide binding by HLA molecules, we analyzed the HLA peptide binding data from the large epitope repository, Immune Epitope Database (IEDB), and further performed in silico analysis to estimate the promiscuity at the population level. Both analyses suggest that an unexpectedly large fraction of HLA ligands (>50%) bind two or more HLA molecules, often across supertype or even loci. These results suggest that different HLA class I molecules can nevertheless present largely overlapping peptide sets, and that “functional” HLA polymorphism on individual and population level is probably much lower than previously anticipated.

HLA class I binding of HBZ determines outcome in HTLV-1 infection

CD8(+) T cells can exert both protective and harmful effects on the virus-infected host. However, there is no systematic method to identify the attributes of a protective CD8(+) T cell response. Here, we combine theory and experiment to identify and quantify the contribution of all HLA class I alleles to host protection against infection with a given pathogen. In 432 HTLV-1-infected individuals we show that individuals with HLA class I alleles that strongly bind the HTLV-1 protein HBZ had a lower proviral load and were more likely to be asymptomatic. We also show that in general, across all HTLV-1 proteins, CD8(+) T cell effectiveness is strongly determined by protein specificity and produce a ranked list of the proteins targeted by the most effective CD8(+) T cell response through to the least effective CD8(+) T cell response. We conclude that CD8(+) T cells play an important role in the control of HTLV-1 and that CD8(+) cells specific to HBZ, not the immunodominant protein Tax, are the most effective. We suggest that HBZ plays a central role in HTLV-1 persistence. This approach is applicable to all pathogens, even where data are sparse, to identify simultaneously the HLA Class I alleles and the epitopes responsible for a protective CD8(+) T cell response.

Frequent associations between CTL and T-Helper epitopes in HIV-1 genomes and implications for multi-epitope vaccine designs

BACKGROUND

Epitope vaccines have been suggested as a strategy to counteract viral escape and development of drug resistance. Multiple studies have shown that Cytotoxic T-Lymphocyte (CTL) and T-Helper (Th) epitopes can generate strong immune responses in Human Immunodeficiency Virus (HIV-1). However, not much is known about the relationship among different types of HIV epitopes, particularly those epitopes that can be considered potential candidates for inclusion in the multi-epitope vaccines.

RESULTS

In this study we used association rule mining to examine relationship between different types of epitopes (CTL, Th and antibody epitopes) from nine protein-coding HIV-1 genes to identify strong associations as potent multi-epitope vaccine candidates. Our results revealed 137 association rules that were consistently present in the majority of reference and non-reference HIV-1 genomes and included epitopes of two different types (CTL and Th) from three different genes (Gag, Pol and Nef). These rules involved 14 non-overlapping epitope regions that frequently co-occurred despite high mutation and recombination rates, including in genomes of circulating recombinant forms. These epitope regions were also highly conserved at both the amino acid and nucleotide levels indicating strong purifying selection driven by functional and/or structural constraints and hence, the diminished likelihood of successful escape mutations.

CONCLUSIONS

Our results provide a comprehensive systematic survey of CTL, Th and Ab epitopes that are both highly conserved and co-occur together among all subtypes of HIV-1, including circulating recombinant forms. Several co-occurring epitope combinations were identified as potent candidates for inclusion in multi-epitope vaccines, including epitopes that are immuno-responsive to different arms of the host immune machinery and can enable stronger and more efficient immune responses, similar to responses achieved with adjuvant therapies. Signature of strong purifying selection acting at the nucleotide level of the associated epitopes indicates that these regions are functionally critical, although the exact reasons behind such sequence conservation remain to be elucidated.

NetCTLpan: pan-specific MHC class I pathway epitope predictions

Reliable predictions of immunogenic peptides are essential in rational vaccine design and can minimize the experimental effort needed to identify epitopes. In this work, we describe a pan-specific major histocompatibility complex (MHC) class I epitope predictor, NetCTLpan. The method integrates predictions of proteasomal cleavage, transporter associated with antigen processing (TAP) transport efficiency, and MHC class I binding affinity into a MHC class I pathway likelihood score and is an improved and extended version of NetCTL. The NetCTLpan method performs predictions for all MHC class I molecules with known protein sequence and allows predictions for 8-, 9-, 10-, and 11-mer peptides. In order to meet the need for a low false positive rate, the method is optimized to achieve high specificity. The method was trained and validated on large datasets of experimentally identified MHC class I ligands and cytotoxic T lymphocyte (CTL) epitopes. It has been reported that MHC molecules are differentially dependent on TAP transport and proteasomal cleavage. Here, we did not find any consistent signs of such MHC dependencies, and the NetCTLpan method is implemented with fixed weights for proteasomal cleavage and TAP transport for all MHC molecules. The predictive performance of the NetCTLpan method was shown to outperform other state-of-the-art CTL epitope prediction methods. Our results further confirm the importance of using full-type human leukocyte antigen restriction information when identifying MHC class I epitopes. Using the NetCTLpan method, the experimental effort to identify 90% of new epitopes can be reduced by 15% and 40%, respectively, when compared to the NetMHCpan and NetCTL methods. The method and benchmark datasets are available at http://www.cbs.dtu.dk/services/NetCTLpan/.