Impact Factor 6.429

The 5th most cited journal in Immunology

This article is part of the Research Topic Role of HLA and KIR in viral infections

SUPPLEMENTAL DATA

Original Research ARTICLE

Front. Immunol., 29 October 2015 | https://doi.org/10.3389/fimmu.2015.00552

HLA preferences for conserved epitopes: a potential mechanism for hepatitis C clearance

imageXiangyu Rao1†, imageIlka Hoof1†, imageDebbie van Baarle2‡, imageCan Keşmir1 and imageJohannes Textor1*
  • 1Theoretical Biology and Bioinformatics, Utrecht University, Utrecht, Netherlands
  • 2Laboratory of Translational Immunology, Department of Immunology, University Medical Center Utrecht, Utrecht, Netherlands

Hepatitis C virus (HCV) infections affect more than 170 million people worldwide. Most of these individuals are chronically infected, but some clear the infection rapidly. Host factors seem to play a key role in HCV clearance, among them are the human leukocyte antigen (HLA) class I molecules. Certain HLA molecules, e.g., B*27 and B*57, are associated with viral clearance. To identify potential mechanisms for these associations, we assess epitope distribution differences between HLA molecules using experimentally verified and in silico predicted HCV epitopes. Specifically, we show that the NS5B protein harbors the largest fraction of conserved regions among all HCV proteins. Such conserved regions could be good targets for cytotoxic T-cell (CTL) responses. We find that the protective HLA-B*27 molecule preferentially presents cytotoxic T-cell (CTL) epitopes from NS5B and, in general, presents the most strongly conserved epitopes among the 23 HLA molecules analyzed. In contrast, HLA molecules known to be associated with HCV persistence do not have similar preferences and appear to target the variable P7 protein. Overall, our analysis suggests that by targeting highly constrained – and thereby conserved – regions of HCV, the protective HLA molecule HLA-B*27 reduces the ability of HCV to escape the cytotoxic T-cell response of the host. For visualizing the distribution of both experimentally verified and predicted epitopes across the HCV genome, we created the HCV epitope browser, which is available at theory.bio.uu.nl/ucqi/hcv.

1. Introduction

Hepatitis C virus (HCV) is a severe infectious agent and affects millions of people per year (1). Approximately 25% of the infected individuals spontaneously resolve the infection, whereas most develop a persistent viral infection, which may lead to liver cirrhosis, liver failure, or hepatocellular carcinoma (2). Although the determinants of HCV clearance, persistence, and disease outcome are not well known, it has been shown in numerous studies that the magnitude, diversity, and quality of T-cell responses play an essential role in controlling the infection and determining the outcome of HCV infection (36). Especially, efficient (CD4+) T-cell responses against non-structural proteins might be important for viral control (7, 8). The importance of T-cell responses is also supported by the finding that certain human leukocyte antigen (HLA) molecules are associated with differential HCV infection outcome. Numerous HCV-HLA association studies have been performed covering different ethnic groups and different HCV genotypes [reviewed in Ref. (9, 10)]. Interestingly, two allele groups, B*27 and B*57, which have previously been associated with slow disease progression in HIV-1-positive individuals [reviewed by Bashirova et al. (11)], are also associated with viral clearance in HCV infection. For HLA molecules in these groups, which we will call protective alleles, several dominantly targeted CTL epitopes have been identified. These epitopes were shown to contain mutations in individuals unable to clear the HCV infection (12, 13) and lead to fitness loss (14, 15) unless restored by compensatory mutations.

The HCV genome encodes a polyprotein precursor of about 3000 amino acid residues, which is processed into several mature structural proteins (Core, E1, E2) and non-structural proteins (NS2, NS3, NS4A, NS4B, NS5A, NS5B, P7) (16). A possible mechanism behind the observed associations between certain HLA class I molecules and HCV infection outcome could be the differential presentation of HCV proteins by certain HLA molecules, as was shown to be the case for HLA associations with HIV-1 (17) and HTLV-1 (18). Such protein targeting could be beneficial for various reasons. For instance, the targeted proteins could be highly expressed, or the timing of their expression could exhibit a critical period which is essential for viral replication. Alternatively, a targeted protein could harbor many conserved epitopes in which mutations would come at a high fitness cost, making it more difficult to escape from the CTL response. Here, we therefore address two related questions: (1) Is there evidence for differential presentation of HCV proteins by protective alleles? and (2) Do protective alleles target conserved regions of the HCV genome?

To address these questions, we analyzed the HCV epitope repertoires of several HLA molecules that have been associated with viral clearance or persistence, employing experimental data from public databases and in silico predictions of HLA-peptide binding to define these HLA epitope repertoires. We found that the protective HLA molecule B*27 shows a preference to present epitopes from the HCV protein NS5B, whereas other HLA molecules show no preferential targeting or target other HCV proteins such as P7. Analyzing the sequence variability of HCV proteins, we found that NS5B harbors the largest fraction of strongly conserved regions among all HCV proteins and that the predicted B*27 epitope repertoire contains the largest amount of strongly conserved epitopes of all alleles that were analyzed. Taken together, our analysis suggests a relationship between the protective potential of an HLA molecule and the degree of sequence conservation of the HCV epitopes targeted by that HLA molecule.

2. Materials and Methods

2.1. Experimentally Verified HCV T-Cell Epitopes

All experimentally verified HCV CD8+ T-cell epitopes restricted by HLA class I molecules were downloaded (October 2014) from two public databases: (1) the Los Alamos HCV immunology database1 [Ref. (19); note that maintenance of this database stopped in 2007] and (2) the Immune Epitope Database and Analysis Resource2 (IEDB (20)). HCV is classified into 7 phylogenetically distinct genotypes (21). As HCV genotype 1 is the most dominant strain worldwide and is well studied (22), we aligned each known epitope to the HCV reference strain H77 (accession number AF009606) using the blastp program (23). Only those epitopes for which an alignment with 100% source coverage could be found were included in the analysis. To make sure that all epitopes we considered are from the HCV strains infecting humans, epitopes identified using HLA transgenic mammalian cells were excluded. This procedure resulted in 398 experimentally verified combinations of an HCV genotype 1 CTL epitope and its known HLA restriction (26 peptides appear in more than one combination). Only 7 epitopes were restricted by HLA-C molecules. To calculate the distribution of CTL epitopes over the HCV proteins for each HLA allele, we limited our analysis to 263 non-redundant experimental epitopes: whenever multiple epitopes aligned to the same positions in the H77 reference strain, we only included the epitope with the highest alignment score. The final set of curated epitopes, which forms the basis for Figure 1, is provided as supplementary data (Data Sheet S1 in Supplementary Material).

FIGURE 1
www.frontiersin.org

Figure 1. NS5B is enriched in experimentally verified epitopes restricted by protective HLA allele groups (HLA-B*27, HLA-B*57). (A) The distribution of (non-redundant) experimentally verified epitopes restricted by the protective alleles. (B) The distribution of all non-redundant HCV epitopes reported in the IEDB and LANL databases. (C) The normalized length of each protein in the complete HCV proteome. The abundance of NS5B in (A) significantly differs from what would be expected based on the distributions in (B) (p = 0.0006) and (C) (p = 0.006; Fisher’s exact test).

2.2. HLA-Peptide-Binding Predictions

We used the artificial neural network-based MHC binding predictor NetMHCpan 2.8 (24) to predict MHC binding affinities for 9-mer peptides from the HCV reference strain H77 (accession number AF009606). For every allele group, 2-digit resolution, we have performed the predictions for the most dominant allele, 4-digit resolution. Alternatively, we predicted MHC binding for common HLA-A and HLA-B alleles by using the Stabilized Matrix Method (25), but we found that this did not affect the conclusions drawn from our analysis so we omit these data from the present article for the sake of simplicity. We also performed a version of our analysis in which we additionally used the algorithm NetChop (26) to predict proteasomal cleavage and transportation by TAP, which are important selection steps during peptide generation for antigen presentation. Also this analysis gave results similar to those reported here.

2.3. Sequence Conservation

Prealigned HCV protein sequences (Core, E1, E2, NS2, NS3, NS4A, NS4B, NS5A, NS5B, and P7 sequences) dated 2008 or older were downloaded (October 2014) from the Los Alamos HCV sequence database. The alignments contain only the “oldest” sequence per patient based on their GenBank submission date and cover sequences from all HCV genotypes. To estimate sequence conservation, the Shannon entropy

E(i)=LpLlog2(pL)

was calculated for each position i of each protein, where pL is the observed frequency of a certain amino acid L at position i in the alignment. The Shannon entropy for position i reaches its maximum value of log2(20) (≈4.3) if all amino acids are equally frequent at that position and is 0 if the position is fully conserved. Different numbers of sequences were available for different HCV proteins, ranging from 466 sequences for NS5B to 918 sequences for NS3. To prevent a bias due to different numbers of sequences in the alignments, we used the Miller-Madow entropy estimator which corrects for such bias (27). We repeated this analysis separately for sequences of HCV genotype 1.

To map sequence conservation of the HCV genome at the spatial resolution that is most relevant for epitope presentation, we applied a running median filter with window size nine to obtain an estimate of the overall conservation of a nine-mer centered at the respective position. The window size nine was chosen because it matches the typical length of peptides presented by HLA class I molecules. The first and last four positions of each protein were not assigned a median entropy, as no nine-mer is centered at these positions. We chose the median filter instead of other alternatives such as a running average because the distribution of entropies is strongly skewed. In order to separate strongly conserved regions from more variable regions, we applied a threshold of 0.037 bits – the lower quartile of the overall entropy distribution in the HCV proteome. Changing this threshold did not affect the qualitative results.

2.4. Statistical Analysis

All statistical analysis was performed using the R environment (28). To analyze the targeting of proteins by HLA alleles, we ranked the 9-mer peptides from each protein by their predicted binding affinities and considered the n best-ranking peptides as epitopes, where n was varied as described in the results. We then compared the resulting distribution of epitopes across proteins to the distribution expected from a random ordering of the peptides. We used a chi-square test to investigate the significance of the difference between the expected and observed distributions and computed standardized residuals for each HLA-protein combination to assess the sources of those differences.

3. Results

3.1. HCV-HLA Associations and Analysis of Known HCV Epitopes

We studied the HLA restriction of all known CTL epitopes from the HCV LANL and IEDB databases (see Methods). Kuniholm et al. (9) conducted a literature review of HCV-HLA associations. Using stringent selection criteria, the authors identified five HLA class I allele groups (B*18, B*27, B*57, Cw*01, and Cw*04) that are strongly associated with either HCV clearance or persistence (Table 1). Interestingly, the two allele groups B*27 and B*57, which are strongly associated with slow disease progression in HIV-1-positive individuals, are also associated with a favorable infection outcome in HCV-infected patient cohorts. Next, Kuniholm et al. (9) tested the association of HLA alleles with HCV viremia in a large cohort of HCV-seropositive women and confirmed all the associations, except for B*18, though the effects were not always significant. Of the list of protective HLA allele groups, multiple CTL epitopes have only been reported for B*27 and B*57 (Table 2). We observed that 6 out of the 10 experimentally verified HCV CD8+ T-cell epitopes restricted by the protective HLA allele groups are located in the HCV protein NS5B, whereas the remaining 4 epitopes are distributed over 4 other HCV proteins (Figure 1A). The share of NS5B epitopes from B*57 and B*27 is larger than would be expected from the distribution of all known epitopes across HCV proteins (Figure 1B) or the protein length (Figure 1C). Thus, this small set of experimentally verified epitopes suggests that the epitopes of protective alleles are concentrated in the NS5B protein.

TABLE 1
www.frontiersin.org

Table 1. HLA class I allele groups with significant associations with HCV viremia as reported by Kuniholm et al. (9).

TABLE 2
www.frontiersin.org

Table 2. Experimentally verified CTL epitopes restricted by protective HLA alleles (used in Figure 1A).

3.2. Preferential Targeting Patterns of HCV Proteins by Class I HLA Molecules

The power of our analysis of experimental epitopes is limited because for many alleles, no CTL epitopes have been reported yet. We therefore used in silico binding prediction as an alternative approach to define HLA epitope repertoires (see Materials and Methods for details). The cohort used by Kuniholm et al. (9) to define protective allele groups (Table 1) was largely infected with HCV genotype 1 (55% with genotype 1a, 29% with genotype 1b, and 4% with genotype 1 of undetermined subtype). Therefore, we assumed that the associations identified by Kuniholm et al. (9) should be valid in the context of HCV genotype 1 infection and used the genotype 1 reference strain H77 from the LANL HCV database (19) for our predictions. The HLA-binding affinities of all HCV peptides of length nine were estimated for several HLA alleles using the prediction tool NetMHCpan (see Materials and Methods). This method predicts an IC50 value of <500 nM for 8 out of the 10 known HLA-B*27 and HLA-B*57 epitopes that are part of Figure 1A (Table 2). Therefore, the in silico predictions should be a valid approach to overcome the limited number of experimentally identified epitopes.

In a first attempt to study preferential targeting of HCV proteins by protective and susceptible HLA alleles, we followed the strategy taken in previous studies on HIV-1 (17) and HTLV-1 infection (18). We ranked the predicted HCV epitopes by their predicted HLA-binding affinity within all HCV proteins and then determined the ranks of the two best-binding peptides derived from each HCV protein among all other HCV epitopes (Figure 2). These ranks are expected to depend on the protein length (Figure 1C): shorter proteins are less likely to contain many top-ranking peptides than long proteins. This intuition was confirmed for the NS4A protein, which is the shortest of the HCV proteins (56 amino acids). Interestingly, the ranks for the second shortest protein P7 (64 amino acids) were much higher and comparable to the other proteins. However, this analysis did not show any convincing differences between protective and susceptible alleles (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Binding affinity rankings of the best epitopes do not differ between susceptible and protective alleles. The relative ranks of epitopes from each HCV protein are plotted for the protective (B*5701, B*2705, and Cw*0101; filled circles) and susceptible (B*1801, Cw*0401; open circles) HLA alleles (see Table 1). The median of each group is indicated by a horizontal line. The relative rank is calculated with respect to the binding scores of all peptides from the HCV proteome and separately for each HLA allele. None of the differences between susceptible and protective alleles are significant at α = 0.05.

We reasoned that focusing only on the top two binders might not be an optimal approach because of the low number of peptides involved, which is expected to produce significant uncertainty with respect to the ranking. Thus, we decided to alternatively look at the whole set of peptides that are predicted to bind each HLA class I molecule. The fraction of peptides bound by each HLA is normally considered to be around 2%. We therefore looked at the top 30–120 binders (corresponding to 1–4% of the HCV proteome, respectively) and compared the observed distribution of proteins among those binders to the expected distribution based on the protein length. Overall, HLAs seem to deviate significantly from uniform protein targeting (p = 0.007, chi-square test, for the 120 binders cutoff). To investigate which HLA-protein combinations gave rise to the difference between the observed protein targeting and a uniform targeting based only on protein length, we investigated the standardized residuals of the chi-square test. Standardized residuals measure the difference between an observed frequency and the expected frequency; a positive residual means that a combination occurs more often than expected, whereas a negative residual means that it occurs less often than expected. The magnitude of the residual gives an indication of statistical significance. For instance, a standardized residual of 2 or −2 means that the difference is twice as large as the standard error in measuring the expected frequency and roughly corresponds to a significance level of α = 0.05. In the following, we will call the residuals enrichment scores.

In contrast to the previous analysis (Figure 2), the enrichment scores revealed differential targeting for many HLA molecules. Specifically, plotting the scores as a function of the cutoff shows that HLA-B*2705, but none of the other protective or susceptible HLA alleles, targets NS5B (Figure 3A). We also performed a two-dimensional analysis of the targeting patterns of all HLA-protein combinations by averaging the standardized residuals across the cutoffs from 60 to 120 (Figure 3B), which showed a rich pattern of protein targeting by many HLA alleles. Including a filter of the predicted HLA ligands for their probability to be generated by the proteasome (26) does not appreciably change these results (data not shown).

FIGURE 3
www.frontiersin.org

Figure 3. HLA-B*27 targets the NS5B protein. (A) Enrichment score (standardized residual; explained in main text and Materials and Methods) as a function of the rank cutoff between epitopes and non-epitopes. For example, at a cutoff of 60, the epitope repertoire consists of the 60 best binders and the proportion of predicted epitopes from a given protein is compared with the proportion that would be expected based on the protein length. An enrichment score of 2 roughly corresponds to significance at α = 0.05 and a score of 2.58 corresponds to α = 0.01 (gray horizontal line). The top 60 binders correspond to the top 2% binders in HCV. (B) Matrix plot of enrichment scores per protein-HLA combination. Enrichment scores in each cell are averaged over cutoffs between 60 and 120 (or the top 2–4%). For clarity, all cells that would not be individually significant at the α = 0.05 are muted (white cells). *These proteins are shorter than 100 amino acids (aa). Therefore, the enrichment scores are based on only a few peptides. The protein lengths are CORE 191aa; E1 192aa; E2 363aa; P7 63aa; NS2 217aa; NS3 631aa; NS4A 54aa; NS4B 261aa; NS5A 448aa; and NS5B 591aa.

Taken together, both the experimentally determined and the predicted epitope repertoires indicate that the protective HLA-B*27:05 allele targets the NS5B protein of HCV.

3.3. Identifying Strongly Conserved Regions of the HCV Genome

Targeting NS5B could be important for several reasons. For example, NS5B could be expressed in infected cells in large amounts or early during the infection, resulting in both cases in an efficient antigen presentation and a good CTL response. Unfortunately, little is known about the expression of HCV proteins in infected cells. Alternatively, since HCV is known to escape CTL pressure by continuous evolution (33), CTL responses to structurally and functionally constrained parts of the HCV genome are likely to be most beneficial. Thus, a possible mechanism how targeting of NS5B by HLA-B*27 would increase the odds of viral clearance could be that NS5B harbors many constrained epitopes. Therefore, we sought to identify conserved regions in the HCV genome and to compare their locations with those of HLA-B*27 epitopes.

The relative variability of HCV proteins was assessed by estimating the average Shannon entropy of all positions for each protein using prealigned amino acid sequences from all genotypes available from the Los Alamos HCV database. Each of these alignments includes only one sequence per patient. In order to allow for direct comparison between the proteins, we corrected for the number of sequences per protein alignment (see Materials and Methods). We then identified strongly conserved regions in the HCV genome as nine-mers in which the median entropy was lower than the lower quartile of the overall entropy distribution; this was the case for about 10% of all nine-mers in the genome. Not surprisingly, strongly conserved regions are not uniformly distributed across HCV proteins (Figure 4). Almost half of the sequence of NS5B and Core are covered by strongly conserved regions, whereas such regions are not found at all in the NS2, NS4A, and P7 proteins. The degree of conservation probably reflects the structural flexibility of these proteins: the structural protein Core, for example, constitutes the capsid of the virus particle and affects several cellular processes (34). For most of the HCV proteins, however, the functional and structural constraints are not yet known. The number of sequences available is not evenly distributed over the different HCV genotypes. Therefore, we could only repeat the conservation analysis on a genotype basis for genotype 1, as the number of available sequences per protein is less than 50 for the other genotypes. The sequence conservation within genotype 1 (plotted downward in Figure 4) largely mirrors the conservation across all genotypes (plotted upward), a notable exception being the NS4A protein which is quite conserved within genotype 1 but highly variable across all genotypes.

FIGURE 4
www.frontiersin.org

Figure 4. NS5B is the most conserved HCV protein. Protein variability was estimated using bias-corrected Shannon entropy scores calculated for each position in HCV protein alignments (see Materials and Methods). For each protein, the plot shows raw entropies (gray bars) and median entropies in nine-mers (black lines) estimated from all sequences (plotted upward) or genotype 1 sequences only (plotted downward; see inset for the scale). Blue arrows point to the centers of experimentally verified HLA-B*27 epitopes (see Figure 1). Strongly conserved nine-mers with low median entropy (see Materials and Methods and main text for details) are highlighted in red. According to the coverage by strongly conserved nine-mers (shown to the right of each plot), NS5B is the most conserved HCV protein (49%). The number of sequences per protein varies between 466 (NS5B) and 918 (NS3).

3.4. Presentation of Strongly Conserved Epitopes by HLA-B*27

We mapped the six experimentally identified HLA-B*27 epitopes to the multiple sequence alignment used for assessing variability of the HCV genome (blue arrows in Figure 4). Four of the epitopes lie within or border of strongly conserved regions of NS5B and NS3, one lies within a less conserved region of NS5B, and one is found in the short highly variable P7 protein. While these sequence conservation of these epitopes thus largely reflects the average sequence conservation of the proteins themselves, the small number of known HLA-B*27 epitopes does not allow meaningful statistical assessment.

Therefore, we also analyzed the localization of in silico predicted epitopes for all HLA alleles used in Figure 3 using the top 30, 60, and 120 binders, which correspond to the top 1%, 2%, and 4%, as explained above. For all three cutoffs, HLA-B*27 is predicted to present the largest number of strongly conserved epitopes of all the HLA alleles that were analyzed (Figure 5), whereas such a preference was not found for HLA-B*57. As was the case for the few experimentally verified epitopes, the epitope conservation level of the predicted epitopes mirrors the protein targeting pattern (Figure 3): indeed, 11 of the 16 strongly conserved predicted epitopes for HLA-B27 lie within NS5B (60-binders cutoff; Table 3).

FIGURE 5
www.frontiersin.org

Figure 5. HLA-B*27 presents the most conserved HCV epitopes. For each HLA supertype, the number of strongly conserved epitopes among the top 30 (black), top 60 (gray), and top 120 (white) predicted binders is shown. HLA-B*27 shows the largest fraction of highly conserved epitopes for all three thresholds.

TABLE 3
www.frontiersin.org

Table 3. Predicted HLA-B*27 epitopes that have both low entropy and high binding affinity.

Taken together, the in silico analysis as well as our analysis based on a limited number of known HCV CD8+ T-cell epitopes both suggest that the protective HLA allele HLA-B*2705 associated with clearance of HCV viremia tends to target conserved regions within the NS5B protein, whereas many other HLA molecules target different (and less conserved) regions of the HCV genome. However, the in silico analysis did not confirm such an effect for the other protective HLA alleles.

4. Discussion

Several studies have shown that both CD4+ and CD8+ T-cell responses play important roles in controlling HCV replication during infection (36), and certain HLA alleles have been significantly associated with differential infection outcome (9). The underlying mechanism behind these associations remains unclear. In case of HIV-1 and HTLV-1, it was shown that the immune system can effectively control infectious diseases by targeting certain pathogen-derived proteins (17, 18). Similarly, other studies found that T-cell epitopes are not equally distributed over the genome of HIV-1 and M. tuberculosis but that they cluster together in conserved regions (3537). We hypothesized that targeting of conserved regions within the HCV proteome by HLA alleles could be a potential mechanism to control viral replication during the acute phase of HCV infection, which could then further determine disease outcome.

To test this hypothesis, we performed a literature study to compile a list of protective and susceptible HLA alleles. This turned out to be rather difficult due to limited data and differences between ethnic groups. Using the five strong HLA class I/HCV associations defined by Kuniholm et al. (9) as a basis for our analysis, we found that experimentally verified CD8+ T-cell epitopes restricted by the protective HLA molecules HLA-B*27 and HLA-B*57 are enriched within the NS5B protein. When carrying out a similar analysis using in silico predicted epitopes, we could only see a preference for presentation of NS5B by HLA-B*2705, but not for HLA-B*57:01. Interestingly, however, for precisely this HLA-protein combination, several epitopes have been identified for which mutations appear to have a high fitness cost (15, 38). In line with previous studies (39), we found that NS5B is the most conserved protein, followed by Core and NS3, while the other proteins harbor more highly variable regions. Combined with these results, our study thus suggests that the protective effect of HLA-B*2705 might be mediated by preferential targeting of conserved epitopes.

Unexpectedly, we found that quite many HLA class I alleles seem to preferentially present epitopes from the P7 protein (Figure 3B), which according to our analysis does not harbor any conserved regions. It is interesting to observe that the HLAs targeting P7 include the two detrimental alleles B*18 and Cw*04. Given that P7 is also the protein whose functional role for HCV infection is most unclear (40), it is tempting to speculate that the P7 targeting prevents those HLAs from inducing more protective responses to other proteins. However, two caveats are that (1) the protective allele Cw*01 also targets P7 and (2) P7 is the second shortest HCV protein, hence the uncertainty involved in ranking its best epitopes is high compared to longer proteins like NS5B.

HCV genotype differences also play a role in potential protective effects. It was shown, for example, that the protective effect of B*27 is dependent on the genotype, with genotype 1-infected patients having T-cell responses against the protective B*27 epitope, and genotype 3-infected patients showing no evidence of T-cell responses to any epitope from the same region in the HCV proteome, which differed in its sequence by three amino acids (13). Our analysis focused on genotype 1 sequences (i.e., we use genotype 1a and 1b reference strains in our in silico analysis and analyze only experimentally verified CTL epitopes from genotype 1 sequences), as HLA associations taken into account here were defined largely in the context of HCV genotype 1 infections. When more sequence and epitope data become available from the other genotypes, it will be possible to repeat our analysis on a genotype basis.

So far, there is no direct evidence that targeting more conserved HCV proteins by CD8+ T cells helps controlling HCV replication. However, it has been widely argued that CD8+ T-cell responses specific for conserved regions of HIV-1 (especially Gag p24) should slow down disease progression (17, 38, 41). As in the case of HIV-1, the most conserved HCV proteins are also the most functionally and structurally constrained ones, which makes it more difficult for the pathogen to escape from host immune control in these proteins than in other regions. In line with this reasoning, it has been shown that a point mutation within a B*27 restricted dominant T-cell epitope in NS5B (pos. 421-429, ARMILMTHF) entails a large fitness cost (15). A study of viral evolution in non-structural proteins during the acute phase of HCV infection showed that mutations in this NS5B epitope are also found and maintained in subjects that are not B*27 positive (42). This lack of reversion over time suggests that additional (compensatory) mutations must have occurred, which do not allow simple reversion. Moreover, a rare mutation in the anchor residue of a subdominant B*27 restricted NS5B epitope (pos. 516-524, GRAAICGKY) was shown to nearly abrogate viral replication, which could be restored by two upstream mutations (38). Interestingly, also for B*57, it has been shown that mutations in an NS5B epitope (pos. 209-217, KSKKTPMGF) substantially impaired viral replication, which was restored by compensatory mutations further upstream (14). Another recent study found that B*57-positive HCV-infected chronic carriers harbor HCV strains with a high frequency of mutations in key residues of this NS5B epitope as well as in an E2 epitope [pos. 158-166, NTRPPLGNW; Ref. (12)]. These epitope variants were also shown to be recognized to a lesser extent in chronic carriers, which indicates that presentation of these epitopes is important in viral clearance. Unfortunately, to our knowledge, none of the other known B*57 or B*27 epitopes were studied in such detail. Still, these data indicate that NS5B is a highly constrained region of the HCV genome and shed some light on why immune responses targeting this conserved protein of HCV might result in a favorable infection outcome.

We wish to stress that the research presented in this article is exploratory and that more clinical data are required to test the mechanistic hypothesis suggested by our results. Specifically, the sixteen predicted epitopes shown in Table 3 should be good CTL targets, as they have low entropy and high predicted binding affinity. This prediction is testable by measuring CTL responses against these epitopes in vivo and correlating the strength of immune escapes against these epitopes to the ability to clear the HCV infection.

Going beyond hepatitis C infections, we believe that the methodology we developed here for measuring HLA protein targeting and epitope distributions may be useful for studying class I responses to other infections as well. To explore the epitope distributions identified in this study for both experimentally verified and in silico predicted epitopes, we have generated the “HCV epitope browser” – an extended, interactive version of Figure 4. This tool is available at theory.bio.uu.nl/ucqi/hcv.

Author Contributions

CK, DB, and JT designed the research. XR, IH, and JT performed the research and analyzed the data. XR, IH, DB, JT, and CK wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

This work was supported by a High Potential Grant from Utrecht University. JT was supported by NWO Earth and Life Sciences grant number 823.02.014.

Supplementary Material

The Supplementary Material for this article can be found online at http://journal.frontiersin.org/article/10.3389/fimmu.2015.00552

Footnotes

References

1. Lavanchy D. The global burden of hepatitis C. Liver Int (2009) 29(Suppl 1):74–81. doi: 10.1111/j.1478-3231.2008.01934.x

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Lauer GM, Walker BD. Hepatitis C virus infection. N Engl J Med (2001) 345(1):41–52. doi:10.1056/NEJM200107053450107

CrossRef Full Text | Google Scholar

3. Cooper S, Erickson AL, Adams EJ, Kansopon J, Weiner AJ, Chien DY, et al. Analysis of a successful immune response against hepatitis C virus. Immunity (1999) 10(4):439–49. doi:10.1016/S1074-7613(00)80044-8

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Grüner NH, Gerlach TJ, Jung MC, Diepolder HM, Schirren CA, Schraut WW, et al. Association of hepatitis C virus-specific CD8+ T cells with viral clearance in acute hepatitis C. J Infect Dis (2000) 181(5):1528–36. doi:10.1086/315450

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lechner F, Wong DK, Dunbar PR, Chapman R, Chung RT, Dohrenwend P, et al. Analysis of successful immune responses in persons infected with hepatitis C virus. J Exp Med (2000) 191(9):1499–512. doi:10.1084/jem.191.9.1499

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Thimme R, Oldach D, Chang KM, Steiger C, Ray SC, Chisari FV. Determinants of viral clearance and persistence during acute hepatitis C virus infection. J Exp Med (2001) 194(10):1395–406. doi:10.1084/jem.194.10.1395

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Ruys TA, Nanlohy NM, van den Berg CH, Hassink E, Beld M, van de Laar T, et al. HCV-specific T-cell responses in injecting drug users: evidence for previous exposure to HCV and a role for CD4+ T cells focussing on nonstructural proteins in viral clearance. J Viral Hepat (2008) 15(6):409–20. doi:10.1111/j.1365-2893.2007.00963.x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. van den Berg CH, Ruys TA, Nanlohy NM, Geerlings SE, van der Meer JT, Mulder JW, et al. Comprehensive longitudinal analysis of hepatitis C virus (HCV)-specific T cell responses during acute HCV infection in the presence of existing HIV-1 infection. J Viral Hepat (2009) 16(4):239–48. doi:10.1111/j.1365-2893.2009.01076.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Kuniholm MH, Kovacs A, Gao X, Xue X, Marti D, Thio CL, et al. Specific human leukocyte antigen class I and II alleles associated with hepatitis C virus viremia. Hepatology (2010) 51(5):1514–22. doi:10.1002/hep.23515

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Mina MM, Luciani F, Cameron B, Bull RA, Beard MR, Booth D, et al. Resistance to hepatitis C virus: potential genetic and immunological determinants. Lancet Infect Dis (2015) 15(4):451–60. doi:10.1016/S1473-3099(14)70965-X

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Bashirova AA, Thomas R, Carrington M. HLA/KIR restraint of HIV: surviving the fittest. Annu Rev Immunol (2011) 29:295–317. doi:10.1146/annurev-immunol-031210-101332

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Kim AY, Kuntzen T, Timm J, Nolan BE, Baca MA, Reyor LL, et al. Spontaneous control of HCV is associated with expression of HLA-B 57 and preservation of targeted epitopes. Gastroenterology (2011) 140(2):686–96. doi:10.1053/j.gastro.2010.09.042

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Neumann-Haefelin C, Timm J, Schmidt J, Kersting N, Fitzmaurice K, Oniangue-Ndza C, et al. Protective effect of human leukocyte antigen B27 in hepatitis C virus infection requires the presence of a genotype-specific immunodominant CD8+ T-cell epitope. Hepatology (2010) 51(1):54–62. doi:10.1002/hep.23275

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Oniangue-Ndza C, Kuntzen T, Kemper M, Berical A, Wang YE, Neumann-Haefelin C, et al. Compensatory mutations restore the replication defects caused by cytotoxic T lymphocyte escape mutations in hepatitis C virus polymerase. J Virol (2011) 85(22):11883–90. doi:10.1128/JVI.00779-11

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Dazert E, Neumann-Haefelin C, Bressanelli S, Fitzmaurice K, Kort J, Timm J, et al. Loss of viral fitness and cross-recognition by CD8+ T cells limit HCV escape from a protective HLA-B27-restricted human immune response. J Clin Invest (2009) 119(2):376–86. doi:10.1172/JCI36587

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Shimotohno K, Tanji Y, Hirowatari Y, Komoda Y, Kato N, Hijikata M. Processing of the hepatitis C virus precursor protein. J Hepatol (1995) 22(1 Suppl):87–92.

Google Scholar

17. Borghans JA, Mølgaard A, de Boer RJ, Kesmir C. HLA alleles associated with slow progression to AIDS truly prefer to present HIV-1 p24. PLoS One (2007) 2(9):e920. doi:10.1371/journal.pone.0000920

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Macnamara A, Rowan A, Hilburn S, Kadolsky U, Fujiwara H, Suemori K, et al. HLA class I binding of HBZ determines outcome in HTLV-1 infection. PLoS Pathog (2010) 6(9):e1001117. doi:10.1371/journal.ppat.1001117

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yusim K, Richardson R, Tao N, Dalwani A, Agrawal A, Szinger J, et al. Los Alamos hepatitis C immunology database. Appl Bioinformatics (2005) 4(4):217–25. doi:10.2165/00822942-200504040-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Vita R, Zarebski L, Greenbaum JA, Emami H, Hoof I, Salimi N, et al. The immune epitope database 2.0. Nucleic Acids Res (2010) 38(Database issue):D854–62. doi:10.1093/nar/gkp1004

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Smith DB, Bukh J, Kuiken C, Muerhoff AS, Rice CM, Stapleton JT, et al. Expanded classification of hepatitis C virus into 7 genotypes and 67 subtypes: updated criteria and genotype assignment web resource. Hepatology (2014) 59(1):318–27. doi:10.1002/hep.26744

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kuiken C, Yusim K, Boykin L, Richardson R. The Los Alamos hepatitis C sequence database. Bioinformatics (2005) 21(3):379–84. doi:10.1093/bioinformatics/bth485

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics (2009) 10:421. doi:10.1186/1471-2105-10-421

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Hoof I, Peters B, Sidney J, Pedersen LE, Sette A, Lund O, et al. NetMHCpan, a method for MHC class I binding prediction beyond humans. Immunogenetics (2009) 61:1–13. doi:10.1007/s00251-008-0341-z

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kim Y, Sidney J, Pinilla C, Sette A, Peters B. Derivation of an amino acid similarity matrix for peptide: MHC binding and its application as a Bayesian prior. BMC Bioinformatics (2009) 10:394. doi:10.1186/1471-2105-10-394

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Nielsen M, Lundegaard C, Lund O, Kesmir C. The role of the proteasome in generating cytotoxic T-cell epitopes: insights obtained from improved predictions of proteasomal cleavage. Immunogenetics (2005) 57:33–41. doi:10.1007/s00251-005-0781-7

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Paninski L. Estimation of entropy and mutual information. Neural Comput (2003) 15:1191–252. doi:10.1162/089976603321780272

CrossRef Full Text | Google Scholar

28. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing (2014).

Google Scholar

29. McKiernan SM, Hagan R, Curry M, McDonald GS, Kelly A, Nolan N, et al. Distinct MHC class I and II alleles are associated with hepatitis C viral clearance, originating from a single source. Hepatology (2004) 40(1):108–14. doi:10.1002/hep.20261

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Thio CL, Gao X, Goedert JJ, Vlahov D, Nelson KE, Hilgartner MW, et al. HLA-Cw*04 and hepatitis C virus persistence. J Virol (2002) 76(10):4792–7. doi:10.1128/JVI.76.10.4792-4797.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Chuang WC, Sarkodie F, Brown CJ, Owusu-Ofori S, Brown J, Li C, et al. Protective effect of HLA-B57 on HCV genotype 2 infection in a West African population. J Med Virol (2007) 79(6):724–33. doi:10.1002/jmv.20848

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Fanning LJ, Kenny-Walsh E, Shanahan F. Persistence of hepatitis C virus in a white population: associations with human leukocyte antigen class 1. Hum Immunol (2004) 65(7):745–51. doi:10.1016/j.humimm.2004.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Guglietta S, Garbuglia AR, Pacciani V, Scotta C, Perrone MP, Laurenti L, et al. Positive selection of cytotoxic T lymphocyte escape variants during acute hepatitis C virus infection. Eur J Immunol (2005) 35(9):2627–37. doi:10.1002/eji.200526067

PubMed Abstract | CrossRef Full Text | Google Scholar

34. McLauchlan J. Properties of the hepatitis C virus core protein: a structural protein that modulates cellular processes. J Viral Hepat (2000) 7(1):2–14. doi:10.1046/j.1365-2893.2000.00201.x

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Fontaine Costa A, Rao X, Lechenadec E, van Baarle D, Kesmir C. HLA-B molecules target more conserved regions of the HIV-1 proteome. AIDS (2010) 24(2):211–5. doi:10.1097/QAD.0b013e328334442e

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Yusim K, Kesmir C, Gaschen B, Addo MM, Altfeld M, Brunak S, et al. Clustering patterns of cytotoxic T-lymphocyte epitopes in human immunodeficiency virus type 1 (HIV-1) proteins reveal imprints of immune evasion on HIV-1 global variation. J Virol (2002) 76(17):8757–68. doi:10.1128/JVI.76.17.8757-8768.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Comas I, Chakravartti J, Small PM, Galagan J, Niemann S, Kremer K, et al. Human T cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconserved. Nat Genet (2010) 42(6):498–503. doi:10.1038/ng.590

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Goulder PJ, Watkins DI. Impact of MHC class I diversity on immune control of immunodeficiency virus replication. Nat Rev Immunol (2008) 8(8):619–30. doi:10.1038/nri2357

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Li H, Hughes AL, Bano N, McArdle S, Livingston S, Deubner H, et al. Genetic diversity of near genome-wide hepatitis C virus sequences during chronic infection: evidence for protein structural conservation over time. PLoS One (2011) 6(5):e19562. doi:10.1371/journal.pone.0019562

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Atoom AM, Taylor NG, Russell RS. The elusive function of the hepatitis C virus p7 protein. Virology (2014) 46(2–463):377–87. doi:10.1016/j.virol.2014.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Kiepiela P, Ngumbela K, Thobakgale C, Ramduth D, Honeyborne I, Moodley E, et al. CD8+ T-cell responses to different HIV proteins have discordant associations with viral load. Nat Med (2007) 13(1):46–53. doi:10.1038/nm1520

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Pfafferott K, Gaudieri S, Ulsenheimer A, James I, Heeg M, Nolan D, et al. Constrained pattern of viral evolution in acute and early HCV infection limits viral plasticity. PLoS One (2011) 6(2):e16797. doi:10.1371/journal.pone.0016797

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hepatitis C virus, sequence conservation, HLA disease association, peptide-MHC binding predictions, CTL responses

Citation: Rao X, Hoof I, van Baarle D, Keşmir C and Textor J (2015) HLA preferences for conserved epitopes: a potential mechanism for hepatitis C clearance. Front. Immunol. 6:552. doi: 10.3389/fimmu.2015.00552

Received: 16 July 2015; Accepted: 15 October 2015;
Published: 29 October 2015

Edited by:

Linda F. Van Dyk, University of Colorado Denver School of Medicine, USA

Reviewed by:

Luka Cicin-Sain, Helmholtz Zentrum für Infektionsforschung GmbH, Germany
Diego A. Vargas-Inchaustegui, National Cancer Institute, USA

Copyright: © 2015 Rao, Hoof, van Baarle, Keşmir and Textor. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Johannes Textor, johannes.textor@gmx.de

Xiangyu Rao and Ilka Hoof have contributed equally to this work.

Present address: Debbie van Baarle, Department Immune Mechanisms, Center for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), Bilthoven, Netherlands